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PREFACE 


Since the mid 1970's, NASA, industry, and universities have worked together to 
conduct important research focused at developing laminar-flow technology that could 
reduce fuel consumption for general aviation, commuter, and transport aircraft by as 
much as 40 to 50 percent. This research, which was initiated under the NASA Aircraft 
Energy Efficiency Program and continued through the Research and Technology Base 
Program, has proved very successful with many significant and impressive results 
having been obtained. 

This symposium was planned in view of the recent accomplishments within the 
areas of laminar-flow control and natural laminar flow and the potential benefits of 
laminar-flow tephnology to the civil and military aircraft communities in the United 
States. The symposium included technical sessions on advanced theory and design tool 
development, wind tunnel and flight research, transition measurement and detection 

techniques, low and high Reynolds number research, and subsonic and supersonic 
research. 
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TRANSITION PHYSICS RESEARCH 


In modern laminar flow flight and wind tunnel research, it is important to 
understand the specific cause(s) of laminar to turbulent boundary layer 
transition. Such information is crucial to the exploration of the limits 
of practical application of laminar flow for drag reduction on aircraft. 

The process of transition involves both the possible modes of disturbance 
growth, and the environmental conditioning of the instabilities by free- 
stream or surface conditions (see fig. 1). The possible modes of 
disturbance growth include viscous (e.g. Tollmien-Schlichting), inviscid 
(e.g. crossflow), and modes which may bypass these natural ones. Theory 
provides information on the possible modes of disturbance amplification, 
but experiment must be relied upon to determine which of those modes 
actually dominates the transition process in a given environment. This 
presentation covers the results to date of research on advanced devices 
and methods used for the study of transition phenomena in the subsonic and 
transonic flight and wind tunnel environments. 


Cause! s) of 

Mode(s) of 


Environmentally conditioned 

transition 

disturbance growth 

X 

instabilities 


Viscous 

Inviscid 

Bypasses 


Turbulence 

Vorticity 

Thermal 

Acoustics 

Vibration 


Figure 1 
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TRANSITION IN THE FLIGHT ENVIRONMENT 


Different experimental environments can produce different transition 
modes. In the subsonic flight environment, the instabilities of practical 
Interest include the Tollmien-Schlichting instability in two-dimensional 
boundary layers, laminar separation-induced inflectional instability, 
crossflow instability in three-dimensional boundary layers, and 
instabilities which may bypass these "natural" modes, caused by surface 
imperfections or acoustics for example {fig. 2). Very few flight 
transition experiments have attempted to document the mode of disturbance 
growth or cause of transition (ref. 1). The data published in nearly all 
past experiments focus on the transition locations without exploring the 
dominant instability(ies) responsible for Initiating the transition 
process. Ultimately it is very important that transition studies be 
conducted at the speeds and altitudes at which the laminar-flow technology 
is intended to be applied. This is important because on laminar airframe 
surfaces which are practical to build and maintain, the transition process 
will be influenced by unit Reynolds number as well as length Reynolds 
number. The simultaneous scaling of both of these parameters is not 
possible; large-scale flight conditions are required for complete 
understanding of transition behavior. 


• Tollmien-Schlichting instability 

• 2-D boundary layers/unswept wings 

• Laminar separation 

• Crossflow instability 

• 3-D boundary layers/ swept wings 

• Bypasses 

• Manufacturing tolerances 

• Engine/airframe noise 

• Insect contamination 


Figure 2 



TRANSITION MEASUREMENT TECHNIQUES 

Laminar, iransitional, and turbulent flows exhibit a variety of changes in 
fluid properties which provide the means to detect the state of the flow. 
Such property changes include skin friction, boundary— layer thickness, 
local turbulence, heat transfer, mass transfer, local flow direction, and 
others. Past research has made extensive use of a wide variety of 
transition measurement methods which rely on the flow properties 
mentioned. Such methods include sublimating chemical, oil flow, China 
clay, acoustic sensors, total surface tubes, liquid films, hot films 
sensors, and liquid crystals. Each technique has advantages and 
disadvantges. Hot films, for example, are extremely useful for continuous 
recording of point-measurements of transition during continuously changing 
test conditions. In the past, however, hot-film transition sensors were 
not widely used to detect the spatial or temporal mode behavior of the 
transition process (e.g. for laminar separation-induced transition 
behavior). Some of the available transition visualization methods such as 
sublimating chemicals have the advantage that they can provide for post- 
flight observations of transition details (ref 2). However, flow 
visualization of transition using sublimating chemicals was not practical 
at flight altitudes much above 20,000 feet. This presentation 
concentrates on recent research on arrayed hot-film sensors for use in 
transition mode measurements and on the development of the liquid-crystal 
method for flow visualization (fig. 3). 


• Sublimating chemicals 

• Oil flow 

• China clay 

• Hot films 

• Acoustic detection 

• Total pressure surface tubes 

• Liquid film 

• Liquid crystals 


Figure 3 
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ARRAYED HOT-FILM DEVICES FOR MEASUREMENT OF TRANSITION LOCATION AND MODE 


Applications of thin hot-fiim devices to measure transition location began 
in the early 1930's, and their use since then has been principally limited 
to single-point measurements of transition location. Thus, except for a 
few experiments where Tollmien-Schiichting wave frequencies were measured 
with hot— films, these devices have not been used to detect transition 
mode. Three arrayed hot-film device concepts have recently been developed 
for detection of transition modes, the multi-element hot-film transition 
sensor, the laminar separation bubble sensor, and the laminar crossflow 
vorticity sensor (ref 3). While these arrayed hot-film concepts can apply 
to either wind tunnel or flight research, the present discussion focuses 
on flight applications. For many flight research applications, through- 
the-surface types of sensor installation are not practical. For such 
situations, the mounting of the sensors and associated signal leads must 
be done on the exterior surface. The problems associated with surface- 
mounted sensors are addressed by these arrayed hot-film devices. The 
multi-element transition sensor and the laminar separation sensor were 
evaluated in flight. The laminar separation sensor has been evaluated in 
wind tunnel tests, and the prototype crossflow sensor has not been tested 
(fig. 4). 


• Multi-element hot-film transition sensor 

• Laminar separation bubble sensor 

• Laminar crossflow vorticity sensor 

Figure 4 
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GATES LEARJET MODEL 28/29 AIRCRAFT 


The flight evaluations of the advanced measurement techniques discussed in 
this presentation were conducted on a NASA-operated Lear Model 28/29 
business jet airplane used at Langley for viscous drag reduction flight 
research (see figure 5). The airplane provides a flight envelope 
including Mach numbers up to 0.805, maximum unit Reynolds numbers up to 
2.65 million per foot, and maximum altitudes up to 51,000 feet. The wing 
on this airplane incorporates a non-production, modified airfoil section. 

Both the wing and the winglet on this airplane have very smooth surfaces 
on which to conduct laminar-flow research. 



Figure 6 
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IS 

MULTI-ELEMENT HOT-FILM TRANSITION SENSOR 

For many research efforts, it is important to study the growth of 
disturbances as neariy along streamlines as possible. Behavior along 
streamlines is of particular concern in three-dimensional flows, for 
example on swept wings and on three-dimensional bodies such as fuselages 
and engine nacelles. For these kinds of investigations, the multi-element 
hot-film sensor overcomes the disadvantages of the individual ("postage 
stamp") sensors often used in the past. The disadvantages of the 
individual surface-mounted hot-film sensors concern installation and 
signal interpretation. With individual hot-film sensors, downstream 
contamination of a sensor by an upstream neighboring sensor creates the 
requirement for staggering of the devices. This staggering of the sensors 
provides information on transition behavior along different streamwise 
stations for each sensor. In three-dimensional flows, staggering of the 
individual sensors may not provide useful transition behavior information. 
Other examples where the multi-element arrayed sensor may be useful 
include investigations of laminar-flow behavior in propeller slipstreams, 
during flight through clouds, or in the presence of acoustic disturbances 
where noise radiation patterns may not be uniform. The sensor is 
fabricated by imbedding the required number and spatial distribution of 
hot-film elements and signal leads between two layers of polyimide film. 
Figure 6 presents the geometric detail of a 25-element sensor design. 
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LEAR 28/29 HOT-FILM ANEMOMETRY INSTALLATION 


Figure 7 depicts the arrangement of sensors on the left wing upper surface 
on the Lear 28/29. The figure shows the location of the multi element 
sensor and the laminar separation sensors which were flown. Transition 
data acquisition with the multi-element sensor was accomplished using an 
electronic switching system which allows rapid switching of all sensor 
elements through six anemometers on board the airplane. Thus at any 
instant, six of the 25 sensor elements could be recorded simultaneously. 
Prior to operation of the multi-element sensor for transition 
investigation, the sensor and surrounding wing surface regions were 
sprayed with sublimating chemicals to ensure that the imperfect sur ace o 
the sensor did not cause transition. The results of these tests showe 
that for altitudes above about 29.000 feet where the cruise unit Reynolds 
number was below about 2.1 million per foot, the sensor roughness did not 
affect transition behavior. 



Figure 7 
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LEAR 28/29 FLIGHT-MEASURED TRANSITION AND PRESSURE DISTRIBUTION 

Figure 8a) illustrates typical transition data which can be acquired with 
the multi-element transition sensor; the pressure distribution (figure 8b) 
was measured with the pressure belts previously illustrated in figure 7. 
Transition onset is located at 41% chord where the first turbulent spikes 
are observed in the hot-film traces. The transitional boundary layer 
produces larger amplitude fluctuations in the signal traces, and at the 
end of the transitional region, the fully developed turbulent boundary 
layer produces the large amplitude uniformly fluctuating hot-film signal 
seen at the 52% chord element. The total length of the transition region 
extends over about eight percent of the chord. This transition region 
length and the character of the hot-film signal are indicative of Tollraien- 
Schlichting (T-S) initiated transition. To assure that this transition 
mechanism is active, analytical predictions of disturbance growth rates In- 
factors) and most amplified T-S frequencies are compared to amplified hot- 
film signal data. 


FLIGHT-MEASURED TRANSITION USING MULTI-ELEMENT SENSOR 

M = 0.77, hp = 39,000 ft, R^. = 8.8 x 106 



Time (sec) 


Figure 8a 
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LEAR 28/29 PRESSURE DISTRIBUTION 

M = 0.80, h = 40427 ft, R' = 1.51 x 10^ ft'\ n = .362 



Chord station, x/c 


Figure 8b 



TOLLMIEN-SCHLICHTING INSTABILITY DETECTION FROM HOT-FILMS 

The identification of the Tollmien-Schiichting disturbance frequencies 
associated with transition is important for validation of laminar 
stability theory, and for ruling out possible bypass-induced initiators of 
transition not related to "natural" T-S amplification. Figure 9 shows a 
hot-film signal trace from an individual hot-film sensor on the Lear 28/29 
wing. With typical signal conditioning, hot-film records such as 
previously shown in figure 8 cannot be used to measure T-S frequencies: 
the hot-film signal must be amplified in order to detect these 
disturbances. The full-scale voltage of the laminar hot-film signal in 
figure 9 is of the same amplitude as the turbulent signals shown in figure 
8. To provide the necessary sensitivity to detect T-S frequencies, the 
hot-film anemometry was operated with as large an overheat ratio as could 
be achieved without melting the sensor element, and the signal was given a 
high external gain prior to recording. Figure 9 shows 0.005 seconds of 
the data on an expanded time scale, revealing the periodic T-S 
oscillations. Spectral analysis of these signals can provide information 
on the most amplified frequencies for correlation with theory. 
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Figure 9 
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PREDICTED AND MEASURED TOLLMIEN-SCHLICHTING WAVE CHARACTERISTICS 

The identification of the Tollmien-Schlichting frequencies associated with 
transition is important for validation of laminar stability theory and for 
ruling out possible bypass-induced initiators of transition not related to 
"natural" T-S amplification. Figure 10 illustrates the power spectral 
density (PSD) analysis of the T-S frequencies near transition measured in 
flight on the Lear 28/29 wing. The spectral analysis of the hot-film 
signal at transition onset shows that the most amplified frequencies occur 
near 4800 Hz. This frequency is in close agreement with theory. The hot- 
film spectrum can also show subharmonic or superharmonic frequencies, 
whereas the linear theory used only accounts for the growth of the 
fundamental frequency. 


Lear 28/29 IV\ = . 79 R = 9. 3 x 10^ A = 17° 




Figure 10 
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CONCEPTUAL OPERATING PRINCIPLE OF THE LAMINAR SEPARATION SENSOR 


For many practical flight conditions of interest, boundary-layer 
transition may occur at the location of laminar separation, in a pressure 
recovery far downstream from the airfoil leading edge. This is most 
likely to occur in dominantly two-dimensional flows at moderate Reynolds 
numbers, but can occur in three-dimensional flows as well, for example on 
swept wings at lower Reynolds numbers. It is particularly important to 
identify the presence of laminar separation and the associated transition 
for calibration of laminar stability theory. Without such knowledge, it 
can be misleading to correlate predicted linear amplification ratios In- 
factors) with experimental transition locations for the purpose of theory 
calibration for use in transition prediction. The laminar separation 
sensor conceptually illustrated in figure 11 provides a simple means for 
the detection of the presence of laminar separation and the associated 
transition. The sensor consists of an array of hot-films (three in the 
illustration) aligned in the direction of the local freestream. The 
center film is electronically heated by a constant temperature anemometer 
(CTA). The remaining two films, upstream and downstream of the center 
element, are incorporated into two legs of a bridge circuit to operate as 
resistance thermometers in a flow reversal meter (FRM) circuit. Thus, 
flow direction is detected when the heat from the center element is 
transported upstream or downstream. 



Figure 11 













MULTI-ELEMENT LAMINAR SEPARATION ARRAYED HOT-FILM SENSOR 
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For detailed studies of laminar separation and transition behavior it is 
desirable to simultaneously measure the locations of separation of the 
laminar boundary-layer, transition on the bubble, and turbulent 
reattachment. These measurements have been accomplished using the array 
of hot-film elements shown in figure 12. With this configuration, the 
elements were operated in sets of three (sub-arrays) at a time. By 
switching the sub-arrays sequentially from the front of the sensor to the 
rear, the details of the laminar separation bubble can be measured. 
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Figure 12 
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LAMINAR SEPARATION SENSOR WIND-TUNNEL EXPERIMENT 

In figure 13, the experimental setup is shown of the laminar separation 
sensor installed on an airfoil model in the Langley 14x22-Foot Wind 
Tunnel. The arrangement includes the multi-element array with 48 elements 
at 1/16 inches spacing in the streamwise direction. The sensor is mounted 
on the upper surface of a 6-inch chord Wortraann airfoil. The model was 
designed and built at Notre Dame University for research on laminar bubble 
behavior planned by Dr. T. Mueller. This cooperative research effort also 
includes planned flight measurements of laminar bubble behavior on the 
same model to be conducted on a sailplane at Texas A&M University by Dr. 

S. Miley. The purpose of the Langley wind tunnel experiment shown in the 
figure was to evaluate the operation of the multi-element laminar 
separation bubble sensor prior to the planned flight experiments. 







Figure 13 
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LAMINAR SEPARATION SENSOR DATA SAMPLE 

An example of transition and flow direction data from the laminar 
separation sensor is shown in figure 14. The figure presents oscillograph 
signals, both from the center elements (CTA's) of each sensor, and from 
the flow reversal meters (FRM's) monitoring the temperature differences 
between the outer two elements of each sensor. The data run began at a 
high angle of attack with transition occurring upstream of the forward- 
most sensor (no. 1). As shown in the figure, the center element in sensor 
no. 1 (CTA-1) indicated a turbulent boundary-layer state at the higher 
angle of attack, with transition moving across the element as angle of 
attack was reduced, until the element indicated a laminar state. The 
corresponding data from FRM— 1 for this sensor showed that no flow reversal 
occurs. Thus, transition initiated by T-S amplification had occurred at 
location no. 1 for these conditions. The same data for sensor no. 3 
(located in the middle of the bubble at zero degrees angle of attack) 
showed that this sensor experienced laminar separation. As angle of 
attack decreases, the passage of transition is indicated by CTA-2, and 
flow reversal is shown by the FRM-2 signals. The traces from CTA-2 behave 
differently from those on CTA-1. As transition moved across sensor no. 2, 
the CTA-2 signal shown in the figure changed from laminar to turbulent 
with a different behavior than was observed from CTA-1. The hot-film 
signal behavior on CTA-2 is typical of laminar-separation-induced 
transition signals. 


CIA Sensor no. 4 
% full scale volts 

FRM Sensor no, 4 
% full scale volts 

CTA Sensor no. 1 
% full scale volts 

FRM Sensor no. 1 
% full scale volts 

Angle of attack 
a . deg 
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Figure 14 




CONCEPTUAL OPERATING PRINCIPLE OF THE CROSSFLOW VORTICITY SENSOR 

Near the leading edge of a swept wing, or on a body with non-axisyinmetric 
pressure gradients, longitudinal vortices can develop in the laminar 
boundary layer, and can cause transition. Locally, these vortices cause 
spatial variations in heat transfer. As illustrated in figure 15, these 
variations in heat transfer can provide the source for detecting the 
wavelength (vortex spacing) and, if nonstationary, the most amplified 
frequencies of crossflow vorticity. A prototype crossflow vorticity 
sensor is illustrated in the figure. The spacing of the individual hot- 
film elements in the sensor is selected to provide accurate wavelength 
determination. Wind tunnel and flight experiments are planned to evaluate 
the performance of this sensor concept. 
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Figure 15 
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SUBLIMATING CHEMICALS FOR BOUNDARY-LAYER TRANSITION VISUALIZATION 


In the past, several methods for in-flight and for wind tunnel 
visualization of transition have been productively used. One most useful 
in recent work has been the sublimating chemical method, illustrated in 
figure 16. It has the advantage of providing, in many situations, data 
which can be recorded on the ground after the flight or after the wind 
tunnel is shut down. Unfortunately, sublimating chemicals have two 
significant disadvantages also. First, only one data point (transition 
measurement) can be obtained with a coating of chemicals; second, for 
altitudes above about 20,000 feet, sublimating chemicals or oil flows are 
Impractical to use. For practical purposes, no means have been available 
in the past for transition visualization at the cruise altitudes of modern 
aircraft. Results will be presented here to illustrate the use of liquid 
crystals to overcome some of the shortcomings of prior means for 
transition visualization. (Fig 17.) 



Figure 16 



ORIGINAL’ PAGE 

BLACK AND WHITE PHOTOGRAPH 







HO\V LIQUID CRYSTALS INDICATE BOUNDARY-LAYER TRANSITION 

crystals are a peculiar state of matter between solid and liouid 
Although they appear as oily liquids, they have certain mechanical 
properties which are siialliar to solid crystals. irparUcIl^lar niuid 
crystals scatter light very selectively. Within a liquid crystal film 

jHchVeUh“o?Z'"her' “ f<=lical fashion. The 

^ helix IS within the range of visible light and wh^n 

subjected to certain physical Influences, thi helix pitch cSnpes and the 

wavelength of the reflected light changes according?r In thl| fasEion 

stress "7 *" to changes m shear 

vapors Heferenc'^^' ferromagnetism, and certain chemical 

P . References 4, 5, and 6 describe recent development and 
applications of the liquid crystal flow visualization method, (pfg 17 .) 



Figure 17 
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TRANSITION VISUALIZATION IN THE LEAR 28/29 WING USING LIQUID CRYSTALS 

Exploratory flight experiments were conducted using many liquid crystal 
formulations on the wing and winglet of the Lear 28/29. Figure 18 shows 
the right wing of the airplane. The surface was prepared for testing by 
spraying on a deep matte flat black paint. A liquid crystal mixture is 
then formulated for the range of surface temperatures anticipated during 
the experiment. The mixture is formulated to have the color-play 
bandwidth (temperature range over which colors will appear) which will be 
most useful for the particular test needs. Narrower bandwidths (about 5 
to 10 degrees Fahrenheit) have proven very useful for providing detailed 
resolution of transition location. Wider bandwidths (as much as 50 or 
more degrees Fahrenheit) have been useful for tests which experience wider 
temperature fluctuations. The liquid crystal material can be either 
brushed on neat, or thinned about 8 to 1 with solvent and sprayed on. By 
either technique, the objective is to obtain a very thin liquid crystal 
film. When applied with the proper film thickness, the liquid crystal 
coating will neither run, nor collect lint and other debris. The figure 
shows transition near the 30- to 35-percent chord location at Mach number 
of .8 at 48,000 feet. The regions of laminar flow appear blue in color, 
and no color appears in the turbulent areas; a few turbulent wedges appear 
emmanating from tape edges near the leading edge. 



Figure 18 
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CORRELATION BETWEEN LIQUID CRYSTALS AND HOT-FILM TRANSITION INDICATIONS 


In order to correlate the transition indications of the liquid crystais 
with another transition method, a liquid crystal coating was applied to 
the region on the left wing of the Learjet, where the multi-element hot- 
film sensor was installed. Figure 19 shows the transition as indicated by 
the liquid crystal coating at Mach = 0.8 at 48,000 feet. On-board real 
time observations of the hot-film signals on the oscilloscopes showed that 
the color change in the liquid crystal coating correlated with the first 
appearance of turbulent spots in the hot-films. Thus, with transition by 
Tollmien-Schlichting amplification, the liquid crystals show the onset of 
the transition region. This differs from the response of sublimating 
chemicals, which show the end of the transition region. 



M = .8, h = 48,000 ft 


Figure 19 
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DYNAMIC RESPONSE OF LIQUID CRYSTAL COLOR CHANGES TO TRANSITION MOTION 

An additional objective of the liquid crystal flight experiments on the 
Lear 28/29 was to qualitatively evaluate the time response of liquid 
crystal coatings to changes in transition location. Figure 20 shows 
elapsed time photographs taken from video record of the transition front 
motions on the winglet. Liquid crystals were brushed on the left winglet 
and flown at an altitude of 17,150 feet at Mach = 0.525. During this 
flight the airplane was oscillated to sideslip angles of +3.9 and -3.5 
degrees with a period of 0.56 seconds. The resulting manuever provided 
variations in local angles of attack and pressure distributions on the 
winglet, causing transition to move from near the leading edge to near the 
70 percent chord location. On this particular airfoil, pressure recovery 
at small angles of attack begins at about 65 percent chord. The observed 
movement of transition was in phase with the sideslip oscillation, 
indicating very rapid time response of the coating. 
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LIQUID CRYSTAL TRANSITION VISUALIZATION 


The use of liquid crystal coatings in future transition research can 
provide data which is very difficuit to acquire by other means. Such 
appiications include the study of transition mode phenomena which are 
affected by changes in freestram disturbance states (e.g. acoustics) or by 
changes in geometry (e.g. sweep). The method overcomes one of the 
iimitations of sublimating chemicals by providing continuous (reversible) 
indications of transition locations for flight conditions which are 
difficult to stabilize, such as at high angles of attack or sideslip. 

Liquid crystals have been used for visualization of shock and shock- 
induced separation in wind-tunnel testing, and for vizualization of 
transition in underwater testing (ref 7). Further evaluations are planned 
of encapsulated liquid crystal coatings which may be useful in 
applications where the thermal response to transition at compressible 
speeds may be sufficiently strong to provide needed information. Examples 
include tests on rotating surfaces (high-speed propellers) where the use 
of neat materials would not be practical. (Fig. 21.) 


• A(j vantages 

• Reversible 

• Rapid time response 

• Non -toxic 

• Low cost 

• Applications 

• Transient test conditions 

• Variable wing sweep 

• Variable acoustic disturbances 

• High altitudes 


Figure 21 


CONCLUDING REMARKS 


Historically, as applications for laminar technology have been pushed to 
higher speeds, the necessary research tools have developed in parallel. 
With the progress in use of arrayed hot-film transition mode sensor 
technology and high— altitude liquid— crystal transition visualization 
methods, we see this trend in tool development continuing. Applications 
for laminar flow throughout future aircraft flight envelopes will require 
detailed knowledge of transition behavior. Transition mode behavior must 
be understood at the flight conditions of interest (i.e. real Mach numbers 
and altitudes) if passive and active laminar flow control applications are 
to be optimized. 
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In recent years a number of measurement techniques have been developed and 
refined for use in fluid mechanic flow diagnostics. These diagnostic tools have been 
primarily used in wind tunnels and, in several instances, have been extended for use 
in flight. The purpose of this paper is to present a brief review of some of the 
advanced techniques that have been developed by eight groups of researchers. The 
eight specific advanced techniques encompass five basic areas of measurements (fig. 
1 ). 

The first section describes the basic measurement of flow quality by a recently 
developed technique which uses a three -wire hot-wire probe to obtain velocity, 
density, and total temperature fluctuation in subsonic compressible flow. The three- 
wire -hot wire probe technique was primarily developed for use in the determination of 
wind tunnel disturbance levels, but, as indicated in this paper, it has also been 
used in flight (refs. 1, 2, 3, and 4). 

The second section describes the basic type of measurement of boundary- layer 
transition detection by two somewhat different hot-film methods and by an optical 
method. The segment including the first hot- film method discusses the recent 
developments with a specialized hot- film system for on-line transition detection in 
cryogenic wind tunnels such as the U.S. National Transonic Facility (refs, 4, 5, 6, 

7, 8, and 9). The segment describing the second hot-film technique discusses the 
measurement and theoretical validation of the existence of Tollmien-Schlichting waves 
that were obtained with a different hot-film system (ref. 10). The third transition 
detection technique discussed involves a recently developed nonintrusive optical 
system which uses a high-resolution interferometer coupled with electronic data 
acquisition system that senses changes in density fluctuation in the boundary layer 
to indicate the transition region (ref. 11). 

The third type of basic measurement encompasses the general area of flow 
diagnostics (i.e., flow quality) made with two types of three component laser 
velocimeter (LV) systems. The first LV system discussed is a single-axis five-beam 
system which indicates excellent resolution with the longitudinal and vertical (u and 
v) components of velocity but poor resolution with the spanwise component (w) of 
velocity due to a high degree of nonorthogonality between the laser beam and w 
component of velocity (refs. 12, 13, and 14). The second LV system discussed 
utilizes three orthogonal beams (and fringe patterns) for use in a "Basic Aerodynamic 
Research Facility" (ref. 15) which has test section side walls and ceiling made of 
glass to accommodate the "orthogonal LV" technique. 

The fourth type of basic measurement is skin friction. This section describes 
a unique moving belt balance which was developed for the direct measurement of the 
wall shear stress (refs. 16, 17, and 18). Measurements with the moving belt balance 
made on the test section side walls of the two Langley facilities are presented. 

The fifth basic type of measurement is flow visualization as applied to 
cryogenic wind tunnels . This section describes various background shadowgraph 
measurements made at different pressures and temperatures in a cryogenic tunnel (ref, 
19) . The modifications to the "test section optical access" that were required for 
an interference free background shadowgraph are discussed. 
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BASIC AREAS OK Kc,ASUREf<ENTS 


Flow quality: three-wire hot-wire probe 

Transition detection: two hot-film techniques and an optical technique 

Flow diagnostics: two LV systems 

Skin friction measurements*, moving belt balance 

Flow visualization: shadowgraphs in cryogenic tunnels 


Figure 1 
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MEASUREMENTS USING THREE-WIRE PROBE TECHNIQUE 


Until recently, very few hot-wire measurements were reported for subsonic 
compressible, transonic, and low supersonic flows. There are several reasons for 
this. First, for a given probe, the heat transfer from a heated wire oriented normal 
to the flow is a function of the three fluid variables: velocity, density, and total 
temperature. This makes the calibration of the hot-wire probe very time consuming. 
Also, since the heat transfer from the hot wire is a function of the combined effect 
of these three fluid variables, it is difficult to separate the individual effects of 
these variables. Finally, the dynamic pressures in these speed ranges are high and 
wire breakage can be excessive due to aerodynamic loads and, of course, any particles 
in the flow can cause wire failure. 

Before the development of the three-wire hot-wire probe (fig. 1, refs. 1 and 
2) , the application of hot-wire anemometry to compressible subsonic and transonic 
flow was similar to that existing in high supersonic flow. In this earlier 
application, the heat transfer from hot wires was considered to be a function of mass 
flow and total temperature (i.e., two-parameter method) rather than the individual 
quantities of velocity, density, and total temperature. The results using a two- 
parameter method made the calibration of hot-wire probes and data reduction somewhat 
less complex, and fluctuating quantities using this technique have been reported for 
compressible subsonic, transonic, and low supersonic flows. The two-parameter method 
represented an advance in hot-wire anemometry and renewed interest in the problems 
encountered in making such measurements in these flow regimes. However, the two- 
parameter method requires several restrictive assumptions and required additional 
measurements to obtain quantitative flow quality results. In contrast, the three- 
wire hot-wire probe technique requires only the assumption usually made for hot-wires 
anemometry with small perturbations (refs. 1, 2, 20, and 21). 

The purpose of this paper (fig. 1) is to: (1) re-examine the heat transfer 

from a hot-wire probe in the compressible subsonic flow regime, (2) describe the 
three-wire hot-wire probe calibration and data reduction technique used to measure 
the velocity, density, and total temperature fluctuations, and (3) present flow 
quality results obtained in the Langley 0.3-Meter Transonic Cryogenic Tunnel (0.3-m 
TCT) and in flight with the NASA JETSTAR from the same three-wire hot-wire probe. 

• Compare wind tunnel and flight measurements with same probe and 
data reduction technique 

• Langley 0,3-meter transonic cryogenic tunnel flow quality (8 by 24-in. 
test section) 

• Jet Star flights at Dryden to obtain atmospheric turbulence (over 
Mojave Desert) 

• Three-wire hot-wire probe method: 

• No restrictive assumptions — such as used with the single 
wire technique 

• Data reduction technique separates velocity, density and 
total temperature fluctuations 


GENERAL HOT-WIRE EQUATIONS 


The mean voltage measured across a given hot wire, oriented normal to the flow, 
is a function of velocity, density, total temperature, and wire temperature. The 
total change in the mean voltage due to a change in the independent variables can be 
obtained from the calculus. If a constant temperature hot-wire anemometer is used, 

Q within the limits of the ability of the feedback amplifier to hold the 
temperature of the wire constant. The equation for the fluctuating quantities and 
the definitions of the three sensitivity coefficients are given in figure 2, 

Attempts usually are made to solve this equation by using the mean square values of 
e*. However, this method results in an equation with six unknowns. The squared 
equation can, in principle, be solved by operating a single wire at six overheats and 
solving a system of six equations for the six unknowns. However, when this was 
tried, difficulties were encountered inverting the 6x6 matrix. 
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EQUATIONS FOR THREE-UIRE HOT-WIRE PROBE 

The form of the general hot-wire equation in figure 2 suggests that the 
instantaneous values of u* , p\ and T^ can be obtained from the solution of an 
equation with three unknowns. The means used for obtaining a solution for the three 
instantaneous quantities was to fabricate a single probe with three hot wires normal 
to the flow, with each wire operated at a slightly different overheat. Then an 
equation can be written for each of the three wires as shown in figure 3 and the set 
of three equations can then be solved for the three unknowns, u* , p* , and T^* , as a 
function of time. Prior to the solution of the sets of equations, the sensitivity 
coefficients must be determined from a calibration which uses the measured values of 
E, u, p, and (ref. 2). Once the fluctuating quantities are obtained conventional 
statistical techniques can be used to obtain mean and root mean square (RMS) values, 
auto correlations, cross products, spectra, etc. 

In principle, the ability to obtain fluctuating quantities using the three -wire 
hot-wire probe technique can be extended into the transonic and low supersonic flow 
regimes without difficulty except those typically associated with lengthy calibration 
and possible wire breakage. 



u', p', T'^, m' { instantaneous values) 


u'P' uT P'T' 

t t (mean of cross product terms) 

From instantaneous values statistical values are determined 
u, P, Tj, m, etc 

Figure 3 
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THREE -WIRE HOT-WIRE PROBE 



The hot-wire probe which was used to obtain flow quality data in the 0.3-m TCT 
and in flight with NASA JETSTAR is shown in figure 4. The probe was fabricated with 
six "needles” made from piano wires which were used to mount the three hot wires 
normal to the flow. The needles, which were about 1 inch long, were glued into a 
section of six-hole ceramic tubing. The ceramic tubing, with the needles and lead 
wires, was glued into a section of 0.250 -inch stainless -steel tubing. Platinum- 
coated tungsten wire having a diameter of 0.0004 inch was mounted between the tips of 
the needles with about 1/4 circle slack in each wire. The wires were about 0.10 inch 
long and about 0.040 inch apart. Each of the three wires on the probe was operated 
with its own constant temperature hot-wire anemometer and signal condition system. 

For the wind tunnel test, the hot-wire probe was mounted in the 8- by 24- inch 
test section of the 0.3-m TCT, (both floor and ceiling 5 percent open), in a mounting 
position normally used for the drag rake. The sensing elements were located midway 
between the centerline of the tunnel and the ceiling at an axial position 
approximately at the center of turntables which were usually used to hold two- 
dimensional models. A shorting probe with three low resistance wires attached to 
pairs of needles was positioned below and rearward of the active probe. The shorting 
probe with leads going to the anemometer was used to balance out any lead resistance 
changes which occurred due to changes in ambient temperature. The resistances of the 
shorting probe and leads were matched to those of the active probe and both sets of 
leads followed the same path from the tunnel test section to the anemometers in order 
to ensure that each set of leads was subjected to the same temperature environment. 
The shorting probe was used to balance out any lead resistance changes which occurred 
due to changes in temperature during the test program. This procedure was followed 
during the wind tunnel (ref. 2) and flight tests. 



Figure 4 
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VELOCITY, DENSITY, AND TOTAL TEMPERATURE FLUCTUATIONS FROM 
0.3-M TCT USING THREE-WIRE HOT-WIRE PROBE 


The RMS values of the velocity, density and total temperature, normalized by 
their respective mean values (fig. 5) were calculated from the large amount of 
detailed data taken during the calibration of the three-wire hot-wire probe. The 
calibration of the three hot wires (ref. 2) (to obtain wire sensitivity coefficient 
indicated in figs. 2 and 3) was started at a condition of maximum dynamic pressure 
(V - 747 ft/sec, p. = 48 psia, and T^ - 504°R) . The three wires were operated at 
ov^heat ratios of 1 6 , 1 . 8 , and 2 . 0) . A square wave signal was used to set I 

approximately equal time constants for the three wires. The upper frequency response 
of the hot-wire anemometry system was about 20 kHz. During the calibration, the 
velocity was varied from 747 to 220 ft/sec in five even increments with the mean and 
fluctuating voltages recorded on line at each velocity. In addition, the fluctuating 
voltages were recorded on odd-numbered tracks of a seven- track tape recorder. As the 
velocity was reduced, the total pressure was decreased to maintain a constant density 
at this fixed total temperature condition. The density was reduced in three steps 
with a similar velocity variation for each of the density steps. The same procedure 
with varying velocities at constant density was repeated for total temperatures of 

522°R and 540°R. _ 

This constant parameter method of calibration, which is a cumbersome procedure, 
is not required if the multiple linear regression technique is used to obtain the , 

sensitivity coefficients. However, the constant parameter method is useful for 
determining the form of the correlation equation to be used for the multiple 

regression method (refs. 2, 20, and 21). ; 

The three sets of fluctuating quantities (fig. 5) were obtained from the j 

simultaneous solution of the hot-wire equation i^n figure 3 for the three hot wires. . 

The data indicate a significant increase in u, p, and with increasing unit j 

Reynolds number for a given Mach number and also indicate an increase with increasing 
Mach number. In general, the velocity and density fluctuations are quite high and * 

have maximum values of about 7.5 and 8.4 percent, respectively. The total | 

temperature fluctuations, which have maximum value of about 0.57 percent, are i 

approximately an order of magnitude lower than velocity and density fluctuations. | 
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MASS FLOW FLUCTUATIONS FROM THREE-WIRE HOT-WIRE PROBE DATA 

In the past, mass flow fluctuations were often used to infer wind tunnel 
disturbance levels in a manner similar to the way velocity fluctuations (as shown in 
fig' 5) were used to quantify wind tunnel disturbance levels. However, the mass flow 
fluctuations are about an order of magnitude lower than velocity fluctuations at the 
higher Mach numbers and Reynolds numbers . An example of this can be seen from the 
calculation of the normalized mass flow fluctuations (fig. 6) from the same three- 
wire hot-wire probe data that were used to determine the velocity and density 
fluctuations in figure 5. The mass flow RMS fluctuations were calculated from the 
sum of the instantaneous velocity and density values. The maximum value of the mass 
flow fluctuations is about 0.78 percent compared to a maximum value of 7.5 percent 
for the velocity fluctuations at the same conditions. The reduction in the mass flow 
fluctuations is predominately due to the fact that at high subsonic Mach numbers the 
major disturbances in the test section were due to upstream propagation of sound from 
the diffuser. For these conditions, the velocity and density fluctuations are anti- 
correlated, with a correlation coefficient nearly equal to -1.0, thereby resulting in 
a lower value for the mass flow fluctuations (refs. 2 and 4). 
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Figure 6 
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NASA JETSTAR WITH THREE-WIRE PROBES MOUNTED ON WING 


Flicht tests which were conducted at the Dryden Flight Research Facility, 

using the Lockheed’ JETSTAR were made with the same NASa'^JETStS (f^r 7) 

0 1-m TCT tests (figs. 4, 5. and 6). A photograph of the NASA JETSTAR (tig. n 

indicates the location of the hot-wire probes on the leading edge on the lef - an 
portion of the wing. 
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CLOSE VIEW OF THREE -WIRE PROBES MOUNTED ON LEADING EDGE OF JETSTAR WING 

A closeup view of the probes (fig. 8) indicates the t 3 rpe of installation used 
to mount the three -wire hot-wire probe holders on the leading edge of the wing. The 
dummy probe was mounted nearer to the fuel tanks where possible flow interferences 
could occur. The active probe was mounted in a more outboard position, in what was 
considered to be a more interference -free flow field. Note that when this photograph 
was taken, the leading- edge flap was lowered to facilitate the installation of the 
two probes on the leading edge of the wing. 

The dummy probe, or shorting probe, contained three low resistance wires across 
the three pairs of needles. This probe was used to balance out any change in lead 
resistance due to ambient temperature changes, in the aircraft. A similar procedure 
with a dummy probe was used during the flow quality test in the 0.3-m TCT (fig. 5 and 
ref. 2). 



THREE -WIRE PROBE DATA ACQUISITION EQUIPMENT IN JETSTAR 


A photograph of the three -wire hot-wire probe anemometry data acquisition 
system mounted in the NASA JETSTAR (fig. 9) was esentially the same acquisition 
system used in the flow quality test in the 0.3-m TCT (ref • 2) . ^ 

available, constant temperature anemometers (CTAs) (DISA 55M system wi^ 55M10CTA 
standard bridge) were used in both the wind tunnel and flight tests pie mean and 

RMS voltages were recorded by hand from the respective instruments during t e i.g . 

The signal conditioner was used to bandpass the signal (1.0 Hz lower fre^ency cuto 

and 20 kHz upper frequency cutoff) and amplify the signal as required. The 

fluctuating and mean voltages were also recorded on a 28 -track tap pcordp upd 
specifically for the flight test. The anemometers were operated by batteries in an 
effort to eliminate 60 Hz electrical noise present in the JETSTAR. 


Oscilloscopes 
_ 

RMS voltage meters 
Mean voltage meters i 


CTAs with standard bridgej] 


Signal conditioners 
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ATMOSPHERIC TURBULENCE FLIGHT DATA FROM THREE-WIRE PROBE 

The RMS values of velocity, density, and total temperature, normalized by their 
respective free-stream quantities (fig. 10), were obtained with the three-wire hot- 
wire probe in a flight with the NASA JETSTAR in January 1982. The data were obtained 
at an altitude of 5000 feet over the Mojave Desert at Mach numbers of 0.2, 0.3, 0.4, 

desert, relative to sea level, ranges from about 2200 to’ 
2400 feet, which meant the flight altitude was about 2700 feet relative to the 
ground. The proximity of the desert floor, coupled with the high winds that 
generally prevail during the time of year of the flight, accounted for the rather 
Mgh level of atmospheric turbulence (i.e., normalized RMS velocity fluctuation). 

The turbulence level and the extremely rough "ride quality" would put these 
measurements in the classification of ”clear air turbulence.” 

Previous measurements of atmospheric turbulence (refs. 22 and 23), in terras of 
normalized RMS mass flow fluctuations, which were made at altitudes of 12 29 and 

39 thousand feet and Mach numbers from 0.28 to 0.80, indicated levels an order of 
magnitude less than those obtained with the three-wire probe (fig. 10). The density 
fluctuations (fig. 10) were somewhat lower than the velocity fluctuations, and the 
total temperature fluctuations were even lower. A comparison of the flight data in 
figure 10 with the wind tunnel data in figure 5 (from the identical hot wire probe) 
indicates that all three fluctuating quantities from the flight data are of at least 
an order of magnitude lower than the wind tunnel data. 
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SPECTRA OF VELOCITY, DENSITY, AND TOTAL TEMPERATURE FLUCTUATIONS 
FROM FLIGHT MEASUREMENTS OF ATMOSPHERIC TURBULENCE 

Spectra were obtained from the flight measurements using the three -wire hot- 
wire probe at an altitude of 5000 feet for a Mach number of 0.5 and a Reynolds number 
about 3 million per foot (fig. 11). The slopes of the spectra for the velocity, 
density, and total temperature fluctuations appear to follow the expected -5/3 slope 
for a region of equilibrium turbulence. In this region, referred to as the inertial 
subrange, there is no production or dissipation and the spectra will have a -5/3 
slope. The 60-Hz noise and its harmonics are present to some extent up to the limit 

of the spectra. 
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COMPARISON OF SEVERAL MEASUREMENTS OF ATMOSPHERIC TURBULENCE 


The atmospheric turbulence measurements made with the three -wire hot-wire probe 
are compared with other atmospheric flight measurements (fig. 12) which include high 
and low altitude clear air turbulence, cumulus turbulence, and thunderstorm 
turbulence (ref. 24). The comparison, which is made in terms of the power spectral 
density versus reduced frequency, indicates good agreement between the three wire 
probe data and the other clear air turbulence data. All the data shown indicate a 
-5/3 slope over the entire range of reduced frequencies. The high and low altitude 
clear air turbulence data (ref, 24) was measured with sensors that had a low pass 
frequency cutoff of 16 Hz. In contrast, three -wire hot-wire probe data had a low 
pass frequency cutoff of 2500 Hz and thus extended out to reduced frequencies that 
were about two orders of magnitude higher than the other clear air turbulence data. 

It is interesting to note that the three wire data and the high and low altitude 
clear air turbulence data fall into a nearly linear band. The fact that the three 
sets of data "line up" is believed to be somewhat fortuitous. 



Reduced frequency, radian/m 


Figure 12 
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INTRODUCTION TO HOT-FILM SYSTEMS FOR CRYOGENIC WIND TUNNELS 


It is well known that the determination of the location of boundary- layer 
transition is necessary for the correct interpretation of aerodynamic data in 
transonic wind tunnels such as the U.S. National Transonic Facility (NTF) and the 
Langley 0.3-Meter Transonic Cryogenic Tunnel (0.3-m TCT) (refs. 5 and 22). In the 
past, hot- film gauges (hot films) have proved to be one of the best devices for 
detecting the beginning and end of boundary- layer transition on models in 
conventional (near-ambient temperature and pressure) wind tunnels (ref. 7). However, 
the types of hot films used in the conventional wind tunnels are unsuitable for use 
in cryogenic wind tunnels, such as the NTF, because of the large range of 
temperatures (cryogenic to ambient) and because of the high Reynolds numbers 
encountered in testing (refs, 8 and 9) (fig. 1). In addition, the wing of a typical 
NTF model of a transport aircraft, with a wing span of 4 to 5 feet, would require as 
many as 100 hot films to adequately locate transition over the wide range of test 
conditions encountered in a typical force and moment test. Furthermore, the hot 
films would have to be located on the surface of 3-D shapes such as leading edges of 
wings and winglets. 

In the late 1970 *s the Douglas Aircraft Company (DAC) developed a vapor 
deposition hot- film, system for transition detection in cryogenic wind tunnels (refs. 
23 and 24). The ability of the DAC hot-film system to detect transition was first 
validated in an ambient temperature test (ref. 6) and later validated in a cryogenic 
test in the 0.3-m TCT (ref. 8). Although a significant amount of transition data was 
obtained with the DAC hot -film system in the 0.3-m TCT, the original goal of on-line 
detection of transition was not achieved. Since the tests in the TCT, Langley has 
developed an improved deposition technique for cryogenic hot films. The new 
technique includes a new dielectric and a new process for the buildup of the hot -film 
substrate. Tests of the new hot films in a low- speed tunnel demonstrated the ability 
to obtain on-line transition data with an enhanced simultaneous hot- film data 
acquisition system. 

Langley is also developing a large-scale vapor deposition system which will 
include two large (5 -foot diameter) chambers that will be used to put hot films on 
large NTF models . 

• Special requirements 

• Operate from ambient to cryogenic condition 

• Nonintrustive (less than critical roughness height) 

• Large number of films and realistic models li.e., 3-0 shapes) 

• Obtain transition on line as adjunct to primary force and moment test 

• A cryo hot film system was developed by Fancher of DAC 

• Validated system in ambient temperature and cryogenic tunnels 

• Since validation tests, Langley has been developing: 

• A new hot film and substrate deposition technique 

• Two new on-line data acquisition systems 

• A system for transition detection of NTF models 

Figure 1 
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MODEL USED IN CRYOGENIC TRANSITION DETECTION STUDY: 0.3-M TCT 

The 2-D airfoil used in the boundary layer transition detection study in the 
0 3-m TCT (ref. 8), with 38 hot films on the upper surface, is shown mounted between 
the two test section turntables (fig. 2). With a technique they had developed (refs. 
6 23 and 24) , DAG vapor deposited the hot films and their respective 7 go 
on a ^electric substrate that was uniformly applied to the entire flow surface o 
the airfoil The surface installation of the vapor deposited hot films S 

uLfwari to 8 microinches thick, well below the critical roughness height required 
to trip the thinnest adiabatic laminar boundary layer that was encountered during 

*^^^*^^ The 9- inch chord model, made of beryllium copper, was an advanced- te^nolop 
superc^^tlLl designed by NASA, designated the NASA SC(3).0712. The model 

was fabricated with 16 internal liquid nitrogen (LN 2 ) cooling passages and 7 
thermocouples, 0.050 below the surface, to monitor the ®°del wall temperature Th 
ability to cool the model to wall temperatures as low as 169 R demonstrated the 
ability of the cryogenic hot-film system to operate at cryogenic conditions. 
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ENLARGED VIEW OF HOT -FILM AND GOLD LEADS 


The hot films, which are extremely small and barely visible in figure 2 can 
easily be distinguished in an enlarged photograph of the hot-film at an x/c of 0.6211 
(fig. 3). The photograph, which has a magnification of about 12, clearly indicates 
the hot-film located between the two gold leads. The gold leads are so thin that a 
portion of the metallic hot films (i.e., the end pads) can be seen through the leads. 

The DAC procedure used to deposit an array of specialized cryogenic hot films 
on the surface of a model (refs. 6, 23, and 24) consisted of first applying a coating 
of dielectric substrate on the model, in this case an epoxy paint, and then polishing 
the substrate to a smooth finish of a uniform thickness of about 0.001 in. (0.0254 
mm). Next, a pattern of hot films was applied to the surface by a vapor deposition 
technique. Each end of the hot film was then connected to a copper wire contact 
(i.e,, model surface penetration) by a set of vapor deposited gold leads. 



Figure 3 
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ENLARGED VIEW OF WIRE PENETRATION CONTACT FOR HOT -FILM LEADS 


An important component of the vapor deposited cryogenic hot- film system is the 
model wire penetrations which are flush to the model surface and serve to carry the 
hot- film signal from the vapor deposited electrical leads, through insulated wires in 
the model, to lead wires outside the model. For 2-D airfoils (refs. 6, 8, and 9), 
the wires run from the sides of the model at the point where the vapor deposited lead 
connects to the wire penetrations. This connection which is a vital link in the hot- 
film system must have a constant, low resistance electrical connection. 

For the transition detection study in the 0.3-m TCT (ref. 8), the wire 
penetration was accomplished by cutting 76 slots (38 in each side) in the upper 
surface of the model and 20 slots in the lower surface to provide lead wires for 38 
upper surface hot films and 10 lower surface hot films (fig. 4) . The wire 
penetration contact points for the hot -film leads protrude above the surface of the 
airfoil about 0.001 in (0.0254 mm) in order to be flush with the surface dielectric 
coating, which is also about 0.001 in thick. The penetration wires were bonded into 
the slots with an electrically nonconducting epoxy. After the installation, each 
wire was checked to ensure that it was isolated from the model. (The 7 tubes on the 
side of the model (fig. 4) are part of the liquid nitrogen cooling system.) 



Figure 4 
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INITIAL APPROACH TO METALIZATION FOR CRYOGENIC HOT -FILM SYSTEM 


For a number of years prior to the hot-film tests in the 0.3-m TCT (refs. 8 and 
9) , Langley had been evaluating various dielectric substates that would be 
considerably thinner than the epoxy paint used in the DAC method and that could be 
applied by a vapor deposition process that would ensure a uniformly thick substrate 
over 3-D shapes. Finally, a vapor deposited polymer called Parylene C was selected 
as the potential dielectric for use on a 3-D model to be tested in the NTF. The 
initial approach was to coat a number of polished aluminum coupons (about 4 in by 4 
in) with various thicknesses of Parylene C (about 0.001 in to 0.003 in) and then 
deposit a chrome/nickel metal coating over the Parylene. It appeared at first the 
approach of depositing Parylene C as a dielectric on a model surface and then 

depositing the hot films on top of the Parylene C in a manner similar to the DAC hot- 

film method (refs. 6, 8, 9, 26, and 27) would be a satisfactory technique for 3-D 
model hot-film application. However, when a metallic coupon with an array of hot 
films vapor deposited by this technique were checked for uniformity of resistance, it 
was discovered that there was a wide variation in the level of resistance. In 

addition, when these same films were powered with a moderately low level of current 

by an anemometer, there was a significant number of hot -film failures. When the 
metalization of a chrome/nickel strip deposited on top of a coating of Parylene C 
was examined under a microscope (ref. 9), a number of cracks were discovered in the 
metals (fig. 5). It was soon realized the reason for the cracks was due to the 
flexible nature of the Parylene C dielectric substrate. The Parylene C has the 
capability to accommodate different linear coefficients of expansion. However, when 
hot films, which have intrinsic stresses, were deposited on the Parylene surface, 
fractures resulted in the metal film (fig. 5). 



Figure 5 
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LANGLEY -DEVELOPED METALIZATION FOR CRYOGENIC HOT-FILM DEPOSITION 

The problem of the fractures in the metalization procedure, indicated in figure 
5, was solved after several months of intensive research and development. The 
solution to the fractures in the metalization was to coat the Parylene dielectric 
substrate with a thin layer of fused silica (Si02) (ref. 9). When this Langley- 
developed technique was examined in the microscope, there were no fractures (i.e., 
failures) in the metal film (fig. 6). This addition of the Si02 "buffer" layer 
enabled the Parylene C to act as a strain isolation pad between the model surface and 
the layer of Si02. With the addition of a layer of SiO^, a dielectric substrate was 
produced that was thermally, chemically, and mechanically stable enough to prevent 
the cracking failure of the hot films at both ambient and cryogenic conditions (ref. 

9). 



Figure 6 
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AIRFOIL WITH LANGLEY -DEVELOPED CRYOGENIC HOT-FILM Pj^--gQQ£. .QUAL rTV 
ON UPPER SURFACE INSERT 

Throughout the transition detection studies in the 0.3-m TCT with the DAC 
cryogenic hot-film system (ref. 8), the on-line voltage output from various hot-film 
anemometers versus time traces on oscillographs and oscilloscopes showed a clear 
distinction between laminar, transitional, or turbulent boundary -layer states. 

However, on-line plots of the RMS voltages versus chordwise location (for an array of 
hot films) could not be obtained in real time due to numerous problems with the hot 
films and the data acquisition system. The results from a posttest analysis, which 
required the use of digitized data obtained from the FM tape recordings made during 
the test, clearly validated the ability of this type of system to detect boundary- 
layer transition at cryogenic conditions. However, it was also apparent from the 
tests that if this type of cryogenic hot- film system was to be used for the on-line 
detection of boundary- layer transition in the NTF, significant advances had to be 
made: (1) in the hot-film metalization process, (2) in the on-line data acquisition 

systems, and (3) in the type of dielectric that was used on the wings of large NTF 
models . 

Since the transition detection study in the 0.3-m TCT (ref. 8), Langley has 
developed a new technique for a cryogenic hot- film metalization process, which 
includes a new type of dielectric substrate (ref. 9 and fig. 6) . In addition, the 
simultaneous hot-film data acquisition system used in the 0.3-m TCT test was 
significantly enhanced (ref. 9) , In order to evaluate both the new Langley hot -film 
deposition technique and the enhanced data acquisition system a low speed wind tunnel 
test was conducted with a model of the NACA 0012 airfoil. One surface of the 12 inch 
chord (12 inch span) had a chordwise array of 30 of the Langley hot films (fig. 7). 

The model also had a chordwise installation of 23 pressure orifices on the surface 
that also contained the 30 hot films. 



Figure 7 
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METALLIC INSERT USED IN AIRFOIL TEST OF CRYOGENIC HOT -FILM SYSTEM 

The model used for the low-speed tests of the Langley cryogenic hot- film 
system, shown in figure 7, was fabricated with an internal core of high strength, 
closed-cell foam insulation and an outer shell of reinforced fiberglass. The 
stainless steel pressure tubing (0,20 in. i.d. and 0.30 o.d.) was routed through the 
internal core and was hand sanded to be flush to the fiberglas surface of the 
airfoil. A slot was cut in the upper surface of the airfoil for a 3 in. by 10 in. 
aluminum insert with a thickness of 0.125 in. on which 30 of the Langley developed 
cryogenic hot films were deposited. 

The 3 by 10 inch aluminum insert (fig. 8) was closely fitted to the slot (i.e., 
no gaps) in the upper surface of the model and worked to conform to the design 
contour of the airfoil. Screws which fit into four countersunk holes in the corners 
of the insert held the insert in place so that the entire upper surface of the 
airfoil had the correct airfoil contour. There were 31 contact points (i.e., wire 
penetrations) on each side of the insert flush to the insert's outer surface. 



Figure 8 
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BACKSIDE OF METALLIC INSERT USED FOR CRYOGENIC HOT-FILM INSTALLATION 


The surface electrical wire penetrations (sometimes referred to as pin 
connections in 3-D model installations), shown in figures 4 and 7, must be installed 
in a manner that does not introduce a surface roughness which could trip the laminar 
boundary layer in a localized area. The backside of the aluminum insert (fig. 9), 
used for the 30 hot films , shows the epoxy compound used to bond the individual 
copper wire to the backside of the insert to ensure the wire penetrations remain 
flush to the surface of the insert. For the low speed tests (ref. 9), which were 
made at a constant ambient temperature, dissimilar metal for the model insert and 
wire penetrations (i.e., aluminum and copper) were used and presented no problem 
because of the near isothermal condition of the model during the test. However, for 
the wire surface penetrations used for cryogenic testing, the wires (or pins) must be 
of the same material to avoid a differential growth between the model and the pins 
when the model temperature varies from ambient condition to cryogenic condition. 





Contact point wires — 
penetration points 


Lead wires to 
anemometers 


Figure 9 
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METALLIC INSERT WITH CRYOGENIC SENSORS INSTALLED 
FOR LOW SPEED WIND TUNNEL TEST 

Once the aluminum insert shown in figures 8 and 9 had been coated with Parylene 
C and Si02, the contact points were cleared of any material that would impede the 
necessary secure electrical connection without exposing any of the surface of the 

aluminum insert to the vapor deposited aluminum leads. Next, the array of 30 hot 

films was vapor deposited on the insert by using a mask fabricated by a 

photolithographic technique. The array of 30 hot films and the 60 surface contact 

points were then connected by vapor deposited aluminum electrical leads (fig. 10). 

The hot films were deposited in two swept arrays (skewed 15° from the model 
chord line) , with 23 films in the array near the leading edge and the remaining 7 hot 
films in the downstream array. The 15° skew reduces the concentration of energy into 
the boundary layer along a given chord line that an unskewed array of films would 
produce. In addition, the skew reduces the downstream influence that a spurious 
turbulent wedge from a given hot film might have on the array of hot films. There is 
a shorted lead at the downstream end of the insert which is used to balance out any 
variation in lead resistance due to temperature changes, which will be an important 
consideration in testing in a cryogenic wind tunnel. 



Figure 10 
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CLOSEUP OF HOT-FILK INSERT INSTALLED IN AIRFOIL 


An enlarged view of the forward portion of the NACA 0012 airfoil shows the 
upstream end of the hot- film insert installed in the upper surface (fig. 11) . The 
hot films can be seen between the aluminum leads, and the electrical contact points 
can be seen at the chordwise edges of the insert. In addition, several of the 
pressure orifices can be seen in the lower right hand portion of the photograph. 

The upstream end of the 10 in. long hot- film insert was located 0.4 in. from 
the leading edge of the airfoil, as measured along the chord line of the symmetrical 
airfoil. The chordwise centerline of the metallic insert was 5.5 in. from the right 
side of the model. The row of pressure orifices was 7.6 in. from the right side of 
the model. The model coordinates were within ±0.001 in. of the design coordinate of 
the NACA 0012 airfoil. 



Figure 1 1 




ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 


369 


INSTRUMENTED AIRFOIL MOUNTED IN LANGLEY SMALL CALIBRATION FACILITY 

The instrumented NACA 0012 airfoil was tested in Langley’s Instrument Research 
Division (IRD) small calibration facility (ref. 9). The facility is an open-circuit 
tunnel that takes in room temperature air, at ambient pressure, from within the 
building where the facility is located and exhausts the air in the same area of the 
building. The tunnel operates over a test section speed range from very low subsonic 
to a Mach number of about 0.25 and is driven by a 50 bp motor that operates a 
squirrel cage blower at either 885 revolutions per min (RPM) or 1780 RPM. At a given 
fan RPM, the tunnel speed is set by regulating 12 radial damper vanes which are 
located at the downstream end of the diffuser. 

The airfoil was tested in a vertical position in the 17 in. by 12 in. test 
section with a chord line (at a - 0°) at the center of the floor and ceiling (fig. 
12). The 12-in. span model was secured between two turntables located in the 
longitudinal center of the test section floor and ceiling. The angle of attack was 
set by the alignment of a finely marked circular scale on the top turntable to a 
finely scribed line on top of the test section that marked the longitudinal flow 
axis. The model was tested at angles of attack of -2,0, and +2 , a^ Mach numbers 
of 0.092, 0.122, and 0.247; and at chord Reynolds numbers of 0.65 x 10 , 0.86 x 10 , 
and 1.70 X 10 . 



Figure 12 


.^tlGlNAT. PAnTj; 

ORIGINAL PAGE DE POOR-^fi^TTY 

BLACK AND WHITE PHOTOGRAPW 


370 


ON-LINE CHORDWISE FLUCTUATING VOLTAGES FOR VARIOUS ANGLES OF ATTACK 


One of the goals of the low- speed wind tunnel tests, with the hot- film 
instrumented NACA 0012 airfoil described with figure 12, was to demonstrate an on- 
line capability to obtain the location of the beginning and end of transition with a 
"specialized" cryogenic hot-film system (ref. 9). On-line transition data, presented 
as the ratio of the RMS voltage to the mean voltage across the anemometer bridge, 
were obtained for the chordwise array of hot films at a Mach number of 0.122 and 
chord Reynolds number of 0.86 x 10° at angles of attack of 0^ and ±2^ (fig. 13). The 
data indicate three distinct regions where the boundary layer is first laminar, then 
transitional, and finally fully turbulent. The transitional region is between the 
beginning and end of transition. The sudden rise in the normalized voltage indicates 
the beginning of transition. In the transition region the normalized voltage reaches 
a peak value and then rapidly decreases to a point where the voltage has a nearly 
constant value which indicates the end of transition. The normalized voltage for the 
laminar boundary layer upstream of the beginning of transition has a value slightly 
less than 0.1 while the normalized voltage for the fully turbulent boundary layer 
downstream of the end transition has a value between about 0.5 and 0.8. As the angle 
of attack increases from -2° to 0^ to +2^, the location of the beginning of 
transition correspondingly moves forward from 51 to 43 to 24 percent of chord. This 
forward movement of the beginning of transition is expected based on the type of 
pressure gradients noted in the pressure distribution which were obtained during the 
test program (ref. 9). 



Figure 13 



ON-LINE CHORDWISE FLUCTUATING VOLTAGES FOR VARIOUS MACH NUMBERS 

On-line normalized fluctuating voltages, which indicated the beginning and end 
of transition, were also obtained from the chordwise films for Mach 

numbers of 0.092, 0.122, and 0.247 at an angle of attack of 0 (fig. 14). The 
fluctuating voltages for the three Mach numbers again indicate the beginning and end 
of transition based on the sharp increase and rapid decrease of normalized voltages 
in the region where the boundary layer is transitional. When the Mach number 
increased from 0.092 to 0.122, there was a slight forward movement in the voltage 
distribution that defined the transitional region, but a clear difference in the 
beginning and end of transition could not be distinguished because of variations in 
the spacing of the operational films in the array. For the two lowest Mach numbers 
the beginning and end of transition appears to be 43 and 59 percent, respectively. 
When the Mach number was increased to 0.247, there was a significant upstream 
movement in the location of the beginning of transition to about 10 percent of chord. 
In addition, the chordwise extent of the transition region more than doubled and the 
level of the normalized voltage for the fully turbulent boundary in was 60 to 70 
percent higher than the normalized voltage for the two lower Mach numbers. 

Thus fluctuating voltage obtained in a real-time mode from a closely spaced 
chordwise array of hot films (figs. 13 and 14) will give a clear indication of the 
beginning and end of transition with the new Langley hot films and an enhanced 
simultaneous data acquisition system. 
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FLUCTUATING CHORDWISE VOLTAGES AND CORRESPONDING 
VOLTAGE VERSUS TIME TRACES 


The fluctuating voltage output from the hot films, which was recorded on FM 
tape with a 28 track tape recorder (ref. 9), was examined on oscilloscopes both 
during (i.e., on-line) and after the tests. From the posttest examinations, selected 
voltage -versus -time traces, that correspond to RMS voltages in figures 13 and 14, 
were obtained with a digital storage oscilloscope and were then plotted directly on 
an analog x-y plotter (ref. 9), The voltage -versus -time traces (fig. 15), which 
correspond to specific values of RMS voltage, at an angle of attack of 0^, a Mach 
number of 0.122, and a chord Reynolds number of 0.86 x 10^, clearly indicate the 
waveform for a laminar boundary layer associated with the low RMS voltages upstream 
of the beginning of transition. Downstream of the end of transition the waveform is 
that of a fully developed turbulent boundary with the associated elevated level of 
RMS voltage. The transition region, between the beginning and end of transition, 
indicates the expected positive voltage spikes causes by turbulent bursts in the area 
where the RMS voltage is increasing and also indicates the negative voltage spikes 
associated with "laminar intervals" in the area where the RMS voltage is decreasing. 
The voltage versus time traces (fig. 15) are typical of those from tests with similar 
cryogenic hot films systems (ref. 9) and from tests with conventional hot films 
systems (ref. 3) , 

The data acquisition system which enabled the realization of on-line transition 
detection consisted of 30 channels of hot- film anemometry, a computer-controlled 
switching network, scanner, digital multimeter, filter, and true RMS meter. All 
components in the system were calibrated very carefully to ensure that each channel 
of hot- film anemometry would have a nearly uniform sensitivity and frequency 
response. Measurements were directly taken with the digital multimeter and 
assimilated in the desktop computer. They were then displayed as a plot of e/E vs 
x/c . 
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TEST OF CRYOGENIC HOT -FILM DIELECTRIC WITH AIRFOIL 
MOUNTED ON TEST SECTION SIDEWALL OF NTF 


Two test entries were made in the NTF to evaluate the durability and 
survivability of the Parylene C dielectric substrate over a wide range of Mach 
numbers, Reynolds numbers, and total temperatures. Each facility entry was about 2 
weeks in duration and remained in the test section during the entire "primary" test 
program. The two qualitative evaluation tests were made with an 8-in. span, 5.575- 
in. chord, 20-percent thick symmetrical airfoil and were mounted on the side wall of 
the test section of the NTF (fig. 16). Two separate, but identical, airfoil sections 
with a 0.0003 -in. thick coating of Parylene C were used in each of the tunnel 
entries. The airfoil sections were mounted on the right side wall of the test 
section well downstream of the primary force and moment model. During both tests, a 
minimum of 4 inches of the outboard portion of the coated airfoil was in the 
freestream flow of the test section, and at highest Reynolds numbers as much as 7 
inches of the airfoil span were in the free stream. During the first entry, the NTF 
was operated in the air mode at Mach numbers from 0.2 to 0.9, a total temperature of 
lOO^F, and a total pressure of about 1 atm. A posttest examination of the Parylene G 
coating indicated no mechanical failure and an electrical check indicated the coating 
would be an excellent dielectric substrate for the hot -film system. During the 
second tunnel entry, a second airfoil, also coated with 0.0003 in. of Parylene C, was 
tested in the cryogenic mode. During the tests, the Mach number varied from, low 
subsonic to 1.15, and the total temperature was as low as -247^F (118 K) , with a 
dynamic pressure as high as 3600 PSF. These test conditions represented a severe set 
of test conditions for the Parylene C dielectric. The results of a posttest 
examination again indicated no mechanical failure of the dielectric substrate at the 
high dynamic pressure - cryogenic condition. An electrical check again indicated the 
Parylene C would serve as an excellent dielectric substrate for the cryogenic hot- 
films system. 



Figure 16 
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DAC - NTF PATHFINDER I WING WITH PROVISIONS FOR THE 
INSTALLATION OF CRYOGENIC HOT FILMS 

The application of the cryogenic hot- film system to three-dimensional models is 
the next step in the development of this measurement technique. It is well known 
that future test programs of airplane configurations in the NTF all have the general 
user requirement of a boundary- layer transition detection capability, preferably on- 
line. At this time, the only cryogenically validated transition detection system is 
the vapor deposited hot-film system as discussed in this section of the publication. 
Currently there are NTF models, in the general class of transport aircraft, that 
contain the necessary surface pin penetrations on the surface to accommodate the 
installation of the Langley cryogenic hot films. One of these NTF models is the 
Douglas Aircraft Company's (DAC) Pathfinder I wing (fig. 17) which has 44 surface 
penetrations pins installed at 4 spanwise stations (i.e., four chordwise arrays of 
hot films) on the upper surface of the right hand span of the wing. The installation 
technique that DAC has developed for hot- film pin penetration produces an extremely 
smooth surface at the point of penetration. The pins for the hot- film leads are 
electrically isolated from the model and are attached to lead wires that connect to 
the hot-film data acquisition system. 



Figure 17 
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VACUUM CHAMBERS AT LANGLEY FOR DEPOSITION OF 
CRYOGENIC HOT FILMS ON LARGE MODELS 

Since there has been a requirement for transition detection on NASA models to 
be tested in the NTF for a number of years, coupled with the DAC requirement that a 
cryogenic hot-film system be vapor deposited on the 51-inch span wing, shown in 
figure 17, Langley has been developing a large-scale deposition system. The large 
scale system will consist of two large vacuum chambers (5 ft i.d.) and the required 
large capacity vacuum system to go with both chambers (fig. 18). The two large 
chambers were obtained from surplus equipment and are being refurbished for use as 
large vacuum chambers . 

One of the large chambers (fig. 18) will be used to deposit the Parylene C 
dielectric substrate on the large wings of NTF models and will have a Parylene C 
system similar to Langley's small scale system (ref. 9). The other large chamber 
(fig. 18), for the metalization/ deposition process, is made from stainless steel and 
will require a much lower vacuxom than the chamber used for the dielectric. The 
system for the large metalization/deposition chamber will be similar to Langley s 
small scale system (ref. 9). In addition, two large 30 horse power mechanical vacuum 
pumps with a pumping capacity of 730 CFM were also obtained from surplus. One of the 
mechanical pumps will provide a sufficiently low vacuum for the chamber used for 
dielectric deposition. The other vacuum pump will serve as the roughing pump for the 
large metalization/deposition chamber which requires a much lower pressure (i.e. , 

10 to 10*° mm of mercury) . A diffusion pump will be used in the vacuum system for 
this chamber to obtain the necessary high vacuum required for metalization/ 
deposition. 



Figure 18 
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I N90-12523 

PREDICTED AND HOT-FILM MEASURED 
TOLLMIEN-SCHLICHTING WAVE CHARACTERISTICS 
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Time Traces Obtained From a Thin Film 

The Tollmien-Schlichting (T-S) instability (ref. 28) is a time -dependent 
instability which can lead to transition of laminar boundary layers on airfoils. 

This paper presents a comparison of theoretical predictions and experimental 
observations of the T-S instability on the NLF(1)-0414F airfoil designed by Viken and 
Pfenninger (ref. 29). The theoretical predictions were obtained using the SALLY 
stability code (ref. 30). The experimental observation of the T-S instabilities was 
accomplished by enhancing the output obtained from the conventional hot film 
transition detection technique (refs. 10, 31, and 32). 

The measurements were made in the Langley Low-Turbulence Pressure Tunnel (ref. 

33) on a 3 -foot chord model of the NASA designed NLF(1)-0414 at - 0.143, R - 3 x 
10^ and a - 3.0^. Pressure orifices on both surfaces were used to determine the 
pressure distribution, and hot-film gages were used to locate the beginning of 
transition. j 

Test results, from the same hot films that were used to detect transition, | 

revealed that T-S waves could be detected by the hot film if the hot-film signal | 

(fig. 1) was adequately amplified. | 

I 

i 

I 

• Tests conducted in Langley's LTPT on 3.0 ft. chord NLF (D-0414 
airfoil at M„= 0. 143, R,. = 3.0 x 10^ 

^ I 

• The hot film system also defined transition region 


• Detection of T-S waves on laminar signal required 
high signal-to-noise ratio 



Figure 1 
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Wave Form and Spectra From a Thin Film 


The spectral analysis from the wave form (fig. 2) clearly shows a dominant 
frequency at 1.4 KHz. The measurements were made at 70 percent of chord on the 
NLF(1)-0414 at - 0.143 and - 3 x 10®. 


Wave form Spectrum 



Figure 2 
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Predicted Tollmlen-Schlichting Wave Characteristics 

The dominant frequency, shown in figure 2, can be seen to fall within the 
theoretically predicted incompressible T-S amplification range for -= 0.143, and 
R - 3 X 10^ at 70 percent of chord on the for the NLF(1)-0414 airfoil (fig. 3). 

This comparison of the predicted and experimental measurement of the T-S instability 
on an NLF airfoil demonstrates that conventional hot films, generally used for the 
detection of boundary- layer transition, can be extended for use to measure the T-S 
most amplified frequency. 


Amplified T-S frequencies at 
x/c = 0.7 and values of n 


T-S freq., Hz 

n 

o 1250 

11.1 

o 1500 

12.2'' 

^ 1750 

10.7 


^Most amplified frequency is 1500 Hz 
at transition location, with n = 12.2 

Measured freq. is 1400 Hz 



Chord station, x/c 


Figure 3 
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THE BOUNDARY- LAYER TRANSITION DETECTOR BASED ON 
DIFFERENTIAL INTERFEROMETRY 

The boundary- layer transition detector, developed by Spectron Development 
Laboratories (see reference 11), uses a differential interferometer to nonintrusively 
detect transition in compressible boundary layers. The heart of the device is a 
highly sensitive differential interferometer capable of detecting optical path length 
differences of less than one - thousandth of the wavelength of the laser light. The 
high sensitivity of the interferometer is due, in part, to a compensator loop which 
is employed to adjust the phase relationship between the two beams of the 
interferometer to optimize the performance of the interferometer and to null out any 
low frequency (less than 1 khz) signals that might result from, for example, model 
vibration. Twin photodetectors generate the output voltage signal as well as 
information for controlling the compensator loop. Electronic signal processing is 
also essential for real-time assessments of the boundary -layer activity. (Fig. 1.) 


• Nonintrusive - two beam method 

• Capable of detecting variations in optical path lengths 
(i.e., density fluctuations) on the order A/1000 

• Transition detector consists of: 

• Optical system 

• Compensator loop 

• Detection system 

• Electronic signal processing 

• The RMS voltage output is proportional to the difference 

in density fluctuations between the two beams (i.e. , initial 
rise indicates onset of transition) 

Figure 1 
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SCHEMATIC OF OPTICAL DETECTION SYSTEM 


The differential interferometer uses a 5-mW helium-neon laser for its light 
source. The beam is then passed through a polarizing plate, which aligns the beam at 
a 45 degree angle from vertical and results in one beam having components of 
polarization in both the vertical and horizontal directions. The Pockell*s cell, the 
key element for the compensator loop, provides a means of adding pathlength to one of 
the components so that the relative phase between the two components of the polarized 
beam is controllable. The beam expander and lens provide a collimated beam to the 
beam splitter and Wollaston (1) . Wollaston (1) splits the original beam into two 
beams of orthogonal polarization which leave the Wollaston prism at predetermined, 
but different angles. These two beams pass through the boundary layer on the model 
at different locations. Differences in density fluctuations in the boundary layer at 
the two model locations manifest themselves as optical path length differences 
between the two beams. Light reflecting back from the model surface returns through 
the lens and Wollaston (1) and is directed by the beam splitter to Wollaston (2) and 
the two photodetectors. The photodetectors generate the signal voltage out and 
provide information for adjusting the Pockell cell for optimum interferometer 
performance. (Fig. 2.) 

Reflecting 



Figure 2 
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TRANSITION DETECTION SYSTEM AT THE BOEING MODEL 
TRANSONIC WIND TUNNEL 


An Important evaluation of the instrument occurred during June of 1986 when it 
was tested in the Boeing Model Transonic Wind Tunnel, which is a pilot tunnel with a 
5- by 7 -inch test section. The tunnel test section does not include a plenum and its 
sidewalls are of optical quality glass. As seen in figure 3, the instrument was set 
up on an optical bench next to the tunnel and the beam entered the tunnel normal to 
the test section side wall. An airfoil was mounted between the bottom and top walls. 
The beam struck the airfoil in a direction approximately normal to its surface. The 
electronic instruments that were used to store and assess the signal data in nearly 
^gal time are seen on the workbench behind the optical bench. 



ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 


ORiOjUNAE"f»A€E IS 

O F POOR QUa utv 


I 

m 


384 



oklGTTgTtfc-PAGEjg 

OF 

NATURAL AND TRIPPED BOUNDARY -LAYER TRANSITION 
PATTERN FROM FLOW VISUALIZATION 

A 6 -inch chord NACA 66-006 airfoil was used during the Boeing tests. Both 
laminar and turbulent flow are known to occur on this airfoil. Concurrently with the 
interferometry measurements, sublimating chemicals (napthalene) were used to 
independently determine where boundary- layer transition was occurring. The tunnel 
was operated at freestream Mach numbers of 0.7 and 0.8, which resulted in values of 
Reynolds number per foot approximately equal to 4 million. Figure 4 illustrates the 
type of pattern seen with the sublimating chemicals for a Mach number of 0.7 and with 
a roughness element attached to the leading edge of the airfoil. 

Several interesting features of the flow are apparent. First, a distinct 
turbulent wedge emanates from the leading edge as a result of the roughness element, 
as marked by the absence of the sublimating chemical. (Absence of whitish chemical 
corresponds to a region of high heat transfer, or high surface shear stress, which is 
usually indicative of a turbulent boundary- layer state.) Second, at a location of 
about 5 inches from the leading edge (corresponding to the 1-inch location on the 
scale), it appears that natural transition is occurring independent of the turbulent 
wedges associated with either the tunnel side walls or the roughness element. 



Figure 4 
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COMPARISON OF INTERFEROMETER SIGNAL AND FLOW 
VISUALIZATION FOR NATURAL TRANSITION 

As an example of the type of data generated by the interferometer, figure 5 shows 
a plot of the RMS level of the signal for natural transition. The freestream Mach 
number was 0.70 and the chord Reynolds number was approximately 2 million. No 
roughness elements were on the model. 

The instrument data, shown by the circles, demonstrate an increase in signal 
fluctuations (indicative of density fluctuations in the boundary layer) from a region 
3.5 inches from the leading edge rearward. This area near 3.5 inches appears to 
coincide with the beginning of boundary- layer transition. The unsteadiness in the 
signal increases until 5.0 inches from the leading edge, where the unsteadiness 
reaches a maximum and where, in fact, the line of demarcation between absence and 
presence of chemicals was photographically determined at an earlier time. After this 
peak in activity, the signal RMS once again falls off as the fully turbulent region 
develops. However, the RMS level in the fully turbulent region appears to be 
approaching a value of about 6 millivolts as opposed to the laminar value of about 2 
millivolts . 



Figure 5 
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COMPARISON OF INTERFEROMETER SIGNAL AND FLOW 
VISUALIZATION FOR FORCED TRANSITION 

A second example of transition detection is given for the case of transition 
triggered by the roughness element. The resulting turbulent wedge is traversed by 
the instrument in such a manner that the instrument moves from outside of the 
turbulent wedge to the centerline of the turbulent wedge. During the traverse, the 
chordwise distance is kept constant at 3.5 inches from the leading edge. The 
freestream Mach number was again equal to 0.70. 

Once again, the RMS level of the signal is plotted on the vertical axis while 
the spanwise position of the beams is plotted along the horizontal axis. A spanwise 
distance of 0 inches is actually l~inch from the centerline of the turbulent wedge 
while a spanwise distance of 1-inch corresponds to the centerline of the turbulent 
wedge. Again, the unsteadiness peak in the optical signal coincides with the 
boundary between the presence and absence of sublimating chemical. It is interesting 
that the interferometer would suggest that the actual turbulent wedge generated by 
the roughness element is not as well defined as the sublimating chemical test would 
suggest. (Fig. 6.) 



Figure 6 
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SUMMARY 


This instrument development program has been funded because of the urgent need 
to measure boundary- layer transition in wind tunnels. In the course of this 
development program, .prototype Instrument was designed, built, end tested^ Recent 
transonic experiments in the Boeing Model Transonic Wind Tunnel show that the 
interferometer results correlate very well with sublimating chemical tests. 
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THREE-COMPONENT LASER DOPPLER VELOCIMETER 
MEASUREMENTS IN A JUNCTURE FLOW 
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INTRODUCTION 


In a wing- fuselage juncture the oncoming turbulent boundary layer on the 
fuselage surface experiences steep adverse pressure gradients as it approaches the 
wing leading edge . These gradients cause the boundary layer to separate ahead of the 
leading edge, resulting in a vortex which rolls up and trails downstream in the 
juncture. The use of laser velocimetry facilitated making detailed measurements near 
the non-planar surfaces forming the juncture, and within near- separating flows close 
to the trailing edge of the wing (refs. 12, 13, and 14). 

The experiment was carried out in the Low- Turbulence Pressure Tunnel (LTPT) at 
the NASA Langley Research Center. The LTPT is a single return, closed-circuit wind 
tunnel which can be operated at pressures from 0.1 to 10 atmospheres. It is capable 
of operating at Mach numbers from 0.05 to 0.50, and unit Reynolds numbers from 
300,000 to 49,000,000/m. The turbulence level is very low because of the nine 
turbulence reduction screens and because of the 17.6:1 contraction ratio. The test 
section is 0.91 m (3 ft) wide, 2.29 m (7.5 ft) high, and 2.29 m (7.5 ft) long. 

Figure 1 gives the details of the experimental setup. 



Figure 1 
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THREE -COMPONENT LDV SYSTEM USED IN LTPT 


OF 




The single-axis, five beam optical configuration uses the standard two-color, 
two -component beam pattern with the two green beams arranged in the horizontal plane 
and the blue beams arranged in the vertical plane (ref. 12). The green beams are 
used to measure the U component and the blue beams to measure the V component . A 
third green beam is placed along the optical axis bisecting the angles between the 
original two green beams. The addition of this beam creates two additional fringe 
patterns, from which the W component is obtained. Bragg cells are used in the two 
outer green beams to separate the three signals obtained from the three green fringe 
patterns. (Fig. 2.) 



Figure 2 
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OPTICAL ACCESS THROUGH SIDEWALL WINDOW 


The optical access to the model was limited to a single window of diameter 0.75 
m. This window was surrounded by a drum 0.978 m in diameter by 0.711 m in depth. 

This further complicated the design of the laser velocimeter (LV) system, resulting 
in a reduced cross-beam angle. Also, an off-axis backscatter collection mode had to 
be used in order to reduce the effect of background reflections from the splitter 
plate. (Fig. 3.) 



Figure 3 
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MODEL MOUNTED IN TEST SECTION 




The LV system described here resulted in a sample volume of diameter 80 
micrometers and a length of 120 micrometers. In order to make measurements very 
close to the juncture, the LV beam system was oriented at an angle of 11 degrees with 
respect to the horizontal wing. Flare from the juncture surfaces restricted the 
closest measurement distance to 3 . 5 mm from the vertical splitter plate, and to 1.0 
mm from the horizontal wing. Tridecane, a kerosene derivative, was successfully used 
as the seeding material in this experiment (ref. 14). (Fig, 4.) 



Figure 4 
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COMPARISON OF LV DATA WITH HOT-WIRE RESULTS 


In order to determine the accuracy of the LV data, measurements made at 
stations located farthest from the juncture have been compared here with the two- 
dimensional turbulent boundary layer results from Klebanoff*s hot-wire measurements 
over a flat plate with zero pressure gradient (ref, 34) . Even though these stations 
were chosen to minimize the influence of the juncture, some three-dimensional effects 
are expected to be present in the LV data. The mean velocity component U shown in 
figure 5 compares very well with the hot-wire results. 


75% chord, z = 1.18 In. (3 cm) 
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Figure 5 
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MEAN VERTICAL VELOCITY 


In a flat plate boundary layer, the mean velocity component V is non-zero 
because of the displacement effect of the boundary -layer growth; it is directed away 
from the plate and has small magnitude. The values of mean velocity component V 
measured by the LV system satisfy this criterion* (Fig* 6.) 


75% chord, z- 1.18 in. (3 cm) 
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Figure 6 
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MEAN SPANWISE VELOCITY 

The mean velocity component W is expected to be zero in a two-dimensional flow. 
But the flow in the juncture can still influence the measurements at the location 
where the data are being compared (fig. 7). This can partially explain the non-zero 
distribution of the W component. 
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Figure 7 
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CHORDWISE VELOCITY FLUCTUATIONS 


In the case of the turbulence intensity, u/u„ (where u is the RMS fluctuating 
component of velocity in the x direction) , the LV data obtained inside the boundary 
layer compare well with Klebanoff's hot-wire results (ref. 34). The value outside of 
the boundary layer is higher than the freestream turbulence intensity of 0.04 
percent, obtained from hot-wire measurements made in the same facility. The value is 
also higher than the LV resolution limit of 1 percent. This can be due to the 
following reason: hot-wire measurements normally filter out the very low frequency 

tunnel oscillations (< 2 Hz) whereas in the case of LV measurements, no such 
filtering is carried out. This will result in the higher- than-normal distribution of 
turbulence intensity u' as measured by the LV system. (Fig. 8.) 

75% chord, z = 1.18 in. (3 cm) 
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Figure 8 
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VERTICAL VELOCITY FLUCTUATIONS 

Compairison of the turbulence intensity, v/\x^ (where v is the RMS fluctuating 
component of velocity in the y direction), with Klebanoff*s results (ref* 34) within 
the boundary layer is very good, but outside of the boundary layer the turbulence 
level is high. This is due to the fact that the lowest turbulence level that could 
be measured with the present LV system (or any other LV systems) is about 1 percent, 
(Fig. 9.) 

75% Chord, z = 1.18 in. (3cm) 
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SPANWISE VELOCITY FLUCTUATIONS 


The overall distribution of the turbulence intensity, w/u^ (where w is the RMS 
fluctuating component of velocity in the z direction) , is much higher than 
Klebanoff’s reference profile (Ref. 34). Part of it can be attributed to the three- 
dimensional effects of the juncture. Also, the narrow cross-beam angle between the 
laser beams that were used to obtain w-related components resulted in lower signal- 
to-noise ratio; this could have introduced some inaccuracies in this and other w- 
related components (fig. 10). 
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TURBULENT SHEAR STRESSES 


There is excellent agreement between LV measurements of the Reynolds stress uv 
and Klebanoff’s hot-wire results (Ref. 34). The accuracy of w-related Reynolds 
stresses and vw was not as good because of the reasons detailed earlier. 

In conclusion, a single-axis, five-beam, three-component laser veloclmeter 
system has been used In a major experiment. Satisfactory results have been obtaine 
with the LV system in the juncture flow. Limited optical access to the tunnel proved 
to be a problem for a three component LV system in determining the third component 
(W). (Fig. 11.) 


75% chord, z = 1.18 In. (3 cm) 
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PURPOSE OF INVESTIGATION 



The development of flow diagnostic techniques has gained renewed momentum in 
recent years due to the rising need to validate Computational Fluid Dynamics (CFD) 
results. While CFD has made great strides in the understanding of certain basic flow 
fields, there are still voids in the understanding of flow physics dealing with 
boundary- layer transition, turbulence, and separation. 

For instance, the requirements of CFD for three component data of mean , 

velocity, turbulence levels, and Reynolds stresses are limited to measurement 
techniques that can handle both mean and fluctuating flow quantities. While laser 
velocimetry (LV) systems have an advantage of being nonintrusive and very good at 
measuring mean velocity in unseparated and separated flow fields, they are weak in 
measuring low disturbance levels. In contrast, hot wire techniques are considered 
good in low turbulence flow fields, yet poor at measuring mean flow quantities. 

These strengths and weaknesses are expected to be compounded in the transonic flow 
regime where compressibility influences the outcome of each measurement. The flow 
regimes where the different instruments have general agreement is also expected to 
change with increased Mach number. These areas of mutual agreement will therefore be 
an area of interest in this test series. Thus, it is the basic purpose of this 
investigation to compare 3-D measurements obtained with a hot wire system to those 
obtained with an orthogonal LV system (fig. 1). j 

Current flow diagnostic research efforts are focusing on the higher order flow j 

field data bases, such as those generated by laser velocimetry, hot wire anemometry, j 

and multi-hole pressure probes. These instruments are being used in studies that j 

range from simple 2-D flow fields, such as a flat plate boundary layer, to complex j 

3-D efforts, which include unsteady vortex flows generated by a delta wing. * 

Recent low- speed comparisons of results obtained with LV and hot wires (refs. i 

35 and 36) have revealed strengths and weaknesses of each instrument. To extend this j 

comparative process to transonic speeds, the Basic Aerodynamic Research Facility was i 

modified for a customized orthogonal 3-D LV system, " 

• Comparison of three-component measurements from hot-wire 
anemometry and orthogonal LV system 

• Hot-wire anemometer suitable for low-level fluctuations 

• Laser velocimeter suitable for higher levels of fluctuations ^ 

and is nonintrusive = 

• Area of mutual agreement of two-methods to be determined - 

Figure 1 
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TEST FACILITY AND FLOW CONDITIONS 


A test plan (fig. 2) has been derived that will focus on the differences in the 

two measurement techniques at both sub- and transonic speeds. The first phase is 

expected to show areas of mutual agreement at low speeds (Mach <0.4), To eliminate 
mean velocity gradient and surface problems the initial test series will be performed 
behind a uniform turbulence grid. The classic turbulence decay format will represent 
magnitude and length scale variations. The second test series will be 

performed in a sub- and transonic 2-D shear flow field in the wake of a flat plate. 

Utilizing the velocity gradient, without a surface, will focus on the control volume 
size of both the LV and the hot wire. Moving upstream to the turbulent boundary 
layer of the flat plate will then introduce a surface effect, again concentrating on 
control volume size and laser flair. The final test series of this comparative study 
will address a 3-D vortical flow field at transonic speeds (0.5 < Mach < 1.3), 


• Atmospheric — continuous tunnel 

• Subsonic tests at Mach numbers from 0.1 to 0.4 for u, v, w from 
both LV and hot wire measurements 

• Transonic tests at JVl^ch numbers from 0.4 to 1.2 for u, v, w from 
LV system and u,p, Tj, from hot-wire system 

• Top and side walls of test section made of glass 

Figure 2 
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TEST SECTION OF BASIC AERODYNAMIC RESEARCH FACILITY 
FOR USE WITH ORTHOGONAL LV SYSTEM 

The test section of the Basic Aerodynamics Research Facility (fig. 3) was 
designed for optical access for an orthogonal 3 component LV system (ref. 15). This 
constraint required the slots to be moved to the corner of the test section. 
Comparison of the performance of the new LV test section to that of the uniform 6 
slot configuration showed minimal deviation in the Mach number distribution through a 
Mach number range up to 1.3. The primary effect was a degradation of the re-entry 

flap performance and diffuser efficiency. 

As part of the initial test plan, a seeding study will be initiated to 
determine particulate tracking ability. This will be performed over the entire Mach 
number range of the tunnel. The stagnation line on a sphere will be used as the test 
case. Spheres of different sizes will be used as computations of the stagnation flow 
may be influenced by model blockage . 



Figure 3 
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OPTICS OF THREE -COMPONENT ORTHOGONAL LV SYSTEM 

The 3-D orthogonal capability of the LV optics (fig, 4) makes the Basic 
Aerodynamics Research Facility a one-of-a-kind transonic facility. The flexibility of 
the paneled test section is enhanced by the flexibility of the LV optical system 
itself. Not only is the system capable of 3-D orthogonal measurements but also can 
be arranged in off-axis backscatter or forward scatter off-axis 3-D orientations. 

This allows the researcher to optimize the optical access to any variety of 2 -D or 3- 
D models. It also gives the researcher the ability to compare different 3-D LV 
optical schemes directly. 



Figure 4 
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LV TRAVERSING SYSTEM TO OFFER FLEXIBILITY IN STUDYING 
DIFFERENT 3-D OPTICAL CONFIGURATIONS 

The LV traversing system is designed to offer flexibility in orienting the 
transmitting and/or receiving optics to either a forward scatter or backscatter 
system (fig. 5). This allows the researcher to optimize the optical access to any 
variety of 2-D or 3-D models. It also gives the researcher the ability to compare 
different 3-D and LV optical schemes directly. 
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SUMMARY OF TEST EXPERIENCE ON TEST SECTION SIDEWALLS 


Figure 2 summarizes the experience at LaRC with UTSI (University of Tennessee 
Space Institute) skin friction balances. The results shown were measured on the test 
section sidewalls of the 0.3 -me ter Transonic Cryogenic Tunnel, or 0.3-m TCT, and the 
Unitary Plan Wind Tunnel. The comparison is presented in the incompressible plane 
since the Mach number ranges from low subsonic to transonic in the 0.3-m TCT, and 
through the supersonic range in the Unitary tunnel. The Karman-Schoenherr flat plate 
skin friction formula is included for comparison. The shaded area represents data 
taken over the history of the Unitary tunnel with floating element skin friction 
balances (references 37, 38, and 39.) 

The present results of testing on tunnel sidewalls should not be used to judge 
the accuracy of the skin friction balances. Rather, the data shown here should 
simply be taken as evidence of operational experience. Unless extraordinary 
precautions are taken, tunnel sidewall boundary layers are not classical flat-plate 
turbulent boundary layers. At a minimum, they are non-adiabatic and affected by wall 
roughness. Figure 3 shows that the 0.3-m TCT data level shown here can be 
represented by an equivalent wall roughness of only .02 mm. Since the ratio of 
boundary layer length to roughness height is the scaling parameter, and the boundary 
layer length increases with tunnel size, the large Unitary Plan Wind Tunnel would be 
only one -fourth as sensitive for the same absolute roughness height. 

The appropriate conclusion to be drawn here is that the balance is capable of 
operation in environments as diverse as the cryogenic, transonic, high- shear rate of 
the 0.3-m TCT, and the high- temperature supersonic environment of the Unitary Plan 
Wind Tunnel. 



Figure 2 
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COMPARISON OF 0.3-M TCT DATA AND WALL ROUGHNESS EFFECTS 

The data shown here (fig. 3) for the 0.3-m TCT are the same as shown on the previous 
plot, but have not been transformed to incompressible coordinates. Curves are shown 
for constant values of distributed roughness. The actual roughness was not measured. 
The approximate formula for rough flat plate flow used to generate the curves was 

Cf - (2.87 + 1.581og(x/€))'^-^ 

where e is the roughness height and x is distance form the leading edge, reference 
40. References 17 and 18 discuss the details of applying this formula to a test 
section wall boundary layer, by calculating an equivalent flat plate length. The 
curve labeled smooth in figure 3 was calculated from the relation, 

Cf - 0.027 / 


also from reference 40. 

The intent of this figure is not to promote a rough wall prediction method or 
fuXly explain the data trends . It does serve to demonstrate the severe effects of 
small roughness heights in a turbulent boundary layer, and that an equivalent 
roughness height of only .02mm is sufficient to match the data. Further discussion 
may be found in references 17 and 18. 

Current plans call for further testing of UTSI balances on a large flat plate 
in the NTF. Floating element balances as well as other types of skin friction 
measuring devices will also be tested for comparison. The surface finish will be 
carefully controlled and extensive boundary -layer profile surveys will be conducted. 
The result will be a boundary layer much better understood than the test section 
sidewall cases . 

The basic feasibility of the UTSI balance to operate in cryogenic conditions 
has been demonstrated. Carefully controlled testing will be required to establish 
limits on accuracy. 

M = 0.7 M = 0.4 

oo OO 

T = 300 K o • 

T = 220 K A 



Figure 3 
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PROPOSED NEW RESEARCH 



Work to be done under a new grant from LaRC to UTSI will Include the use of 
fiber optics to read the movement of the belt rather than strain gages. This 
technique is expected to be less sensitive to temperature changes, thus simplifying 
balance calibration and use at cryogenic temperatures. Also electrical noise and 
error due to gage heat will be eliminated. Finally, eliminating the strain gage 
removes the primary barrier to miniaturization. 

A concurrent effort is also under way to fabricate these balances using the 
relatively new wire -cut method (electron discharge machining using a wire for an 
electrode). The use of the wire-cut technique has the potential to lower the cost of 
the balances by reducing the part count and simplifying the assembly procedure. The 
combination of the wire -cut technique and the use of fiber optics may reduce the cost 
per balance sufficiently to allow them to be tailored to a specific test. It is 
conceivable that as many as 10 balances could be dedicated to a single airfoil 
configuration. A typical layout is shown in figure 4 on a 14 percent thick, 10-inch 
chord airfoil. Each balance would be contoured to match the airfoil surface. This 
would simplify measurement of the location of transition, separation and transonic 
shocks , 



41 1 


Figure 4 


N90-12528 


RECENT FLOW VISUALIZATION STUDIES 
IN THE 0.3-m TCT 


Walter L. Snow, Alpheus W. Burner, and William K. Goad 
Langley Research Center 
Hampton, Virginia 


412 


SCHEMATIC OF TYPICAL FLOW VISUALIZATION SETUP IN 0.3-M TCT 


Light beams are altered by refractive index changes and indeed flow induced 
refractive index changes provides the impetus for conventional visualization 
techniques such as schlieren and shadowgraph. Unfortunately effects related to the 
flow can be masked by refractive index inhomogeneities external to the test section. 
For high-pressure cryogenic facilities this is especially true in the interface 
between the region of test and the more benign environment where equipment is housed. 
Since the refractive index is directly proportional to density, it depends upon 
pressure and temperature through the ideal gas law. Even mild temperature 
inhomogeneities can be detrimental when coupled with high pressures as occurs in the 
plenum. The simple shadowgraph scheme depicted in figure 1 was used to assess the 
flow quality of the Langley 0.3 -Me ter Transonic Cryogenic Tunnel. The light source 
was a pulsed (15 microsecond) 75-watt xenon arc lamp focused onto a 0.51 mm diameter 
pinhole. Light was collimated with an f/8 mirror with a 1.22-m focal length. Data 
were recorded on film with a 70mm format reflex camera. 
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REPRESENTATIVE SHADOWGRAPH RESULTS PRIOR TO IMPROVEMENTS 

The simple shadowgraph scheme depicted in figure 1 was used to characterize the 
flow quality of the Langley 0.3 -Me ter TCT with results shown in figure 2. This 
collage of photos was taken at Mach 0.65 conditions with stagnation pressure and 
temperature as indicated. The completely uniform light field which characterizes the 
pretest condition becomes increasingly mottled and jumbled as temperature decreases 
and pressure increases. This behavior is typical of earlier optical experiments 
conducted in this facility (see refs. 41, 42, and 43). 
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Figure 2 
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MODIFIED SHADOWGRAPH SYSTEM TO EXCLUDE PLENUM EFFECTS 


Several of the tests devised to determine the source of the image degradation 
are described in reference 19. The crucial test involved installation of tubes to 
penetrate the plenum. The tubes could be purged or evacuated to allow direct viewing 
of the test section. Figure 3 shows the same layout as figure 1 with the hatched 
area designating the region impacted by the isolation tubes, A bellows type of 
arrangement was required to accommodate the large differential expansion experienced 
in the structure of the facility. Details of the assembly are found in the 
reference . 


Top view — With evacuation in equipment pods (cross hatched) 



Figure 3 
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EFFECT OF PLENUM EVACUATION ON SHADOWGRAPH RESULTS 

When the penetration tubes were evacuated the quality of the shadowgraphs 
improved dramatically. Results for the worst test condition of figure 2 (T,. - 180 R, 
p - 5 atm, Mach 0.65) are shown in figure 4. Residual "graininess" in the evacuated 
case is believed to be due to the flow field itself. Earlier difficulties were 
brought about by density inhomogeneities in the regions of strong temperature 
gradients especially in those areas under high pressure. Provisions for evacuated 
paths to and from the test section should be seriously considered in future flow 
visualization efforts for facilities of this type. In scaled up cryogenic facilities 
such as the National Transonic Facility equipment is completely housed within the 
plenum area (ref. 42). Much of the critical optical path would be in the controlled 
environment of the equipment enclosures. Nevertheless, regions near windows where 
extreme temperature gradients occur should be evacuated if possible to lessen the 
effects of convection induced mixing. 
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ABSTRACT 


Gortler vortices arise in laminar boundary layers along 
concave walls due to an imbalance between pressure and centri- 
fugal forces. In advanced laminar-flow control ( LFC) super- 
critical airfoil designs (reference 1), boundary-layer suction is 
primarily used to control Tol Imien-Schl icht ing instability and 
cross-flow vortices in the concave region near the leading edge 
of the airfoil lower surface. The concave region itself is com- 
prised of a number of linear segment s pos it ioned tolimitthe total 
growth of Gortler vortices. Such an approach is based on physi- 
cal reasonings but rigorous theoretical justification or experi- 
mental evidence to support such an approach does not exist. An 
experimental project was initiated at NASA Langley to verify this 
concept. In the first phase of the project an experiment was 
conducted on an airfoil whose concave region has a continuous 
curvature distribution- Some results of this experiment have 
been previously reported (references 2 and 3) and significant 
features are summarized in this paper. (Figures 1-3). 


MOTIVATION 


• LFC SC airfoil development 

• 8-FT experiment 

• Influence of curvature distribution 



• Influence of suction 

• Study of Gdrtler vortices in external flow 


Figure 1 


BACKGROUND 


• Previous experiments used 'bumps' on opposite walis to 

ensure attached flow in the concave region or 
were conducted In channeis 

• Results indicate invariance of waveiength with 

freestream parameters 

• Information on effects of Suction, Pressure gradient, 

and Curvature distribution not available 

Figure 2 

PRESENT EXPERIMENT 

• Airfoil model 

— Continuous, smooth curvature distribution 

- Concave region 0.175 <x/c <0.275 

- Perforated titanium panel 0.17 < x/c < 0.225 

- Convex region 0.275 <x/c <0.50 

• Flow visualization 

• Laser velocimetry 

Figure 3 
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AIRFOIL MODEL AND TEST SET-UP 


The 1. 83-m airfoil model was tested in the NASA Langley Low- 
Turbulence Pressure Tunnel ( LTPT) . The model consists of two 
parts - a structural element and a test insert (Figure 4). With 
this design several test region geometries can be examined. 
Attached laminar boundary-layer flow is maintained by means of 
suction through a perforated t itanium panel located in the compres- 
sion part of the concave region. The suction region is divided 
into three spanwise suction strips» each controlled independently 
by its own needle valve. The 10— percent chord flap is used to 
control the stagnation point location. The experiment was con- 
ducted at atmospheric pressure and the chord Reynolds number was 
varied from 1.0 million to 5.9 million yielding a Gortler number 
of up to 16. 




Figure 4 
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FLOW VISUALIZATION 


X>RIG1NAL 

OF 


A thin layer of solid white biphenyl sublimating chemical 
was sprayed over the black model surface. The flow pattern is 
made visible due to the differential surface shear stress distri- 
bution under the layer of Gortler vortices which leave behind a 
trace of alternating white and black bands on the model surface. 
Elapsed times of about 30 minutes to 1 hour were required for the 
pattern to emerge clearly. A pair of white and black bands con- 
stitutes a wavelength. The wavelength was determined by taking 
the average of the number of pairs of bands over a 15 cm to 45 cm 
span. Beyond the concave region, the bands decrease considerably 
in contrast but are visible back to the jagged transition line. 
The decrease in contrast is attributed to damping in the convex 
region. (Figure 5). 



Rc “ 2.24 Million 



Rc = 3.67 Million 


Figure 5 
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SPANUISE VARIATION OF NORMALIZED STREAMWISE VELOCITY 


A specialized, single— axis, three— component laser veloci- 
meter was used to study the flow field in the test region (refer- 
ence 2). The velocity field was measured by scanning the laser 
control volume spanwise in steps of 0.05 cm at a constant nominal 
height above the model surface. The control volume was then 
moved normal to the surface to the next measurement height and 
scanned back. This process normally took approximately five 
hours per chord location. The minimum and the maximum chord 
Reynolds numbers for laser velocimeter measurements were 1.0 and 
3.67 mil 1 ions, respect ively. The spanwise distributions of mean 
streamwise velocity, normalized with respect to the local 
boundary-layer edge velocity, are given in Reference 2 for 
different chord locations and chord Reynolds numbers. Typical 
results at x/c = 0.25 are shown in Figures 6(a>and6(b) for the 
lowest and the highest test chord Reynolds numbers. The pertur- 
bations are observed to vary across the boundary layer with an 
essentially constant periodicity along the span for a given chord 
Reynolds number. Appreciable difference in periodicity is clear- 
ly seen with a change in chord Reynolds number. Maximum spanwise 
fluctuations of as much as 45-percent in streamwise velocity can 
be observed in Figure 6( b) . 


Rc = 1.0 X 10®, G = 5.679, x/c = 0.25 


x/c-0.25, Rc - 3.67 million 


u/U 




O y - 0.22 cm 

□ y - 0.18 cm 

O y “ 0.13 cm 

A y • 0.08 cm 

y - 0.04 cm 


( a ) 


( b) 


Figure 6 
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NORMALIZED STREAMWISE VELOCITY PERTURBATION 


Normalized streamwise velocity perturbation functions at 
various chord locations are shown in Figure 7<a> for a chord 
Reynolds number of 1.0 million. The maximum amplitude of the 
streamwise perturbation grows to a maximum at the end of the 
concave region ( x/c = 0.275) and then falls rapidly in the convex 
region. Linear theory was used to compute theoretical eigen- 
functions corresponding to the measured wavelength. For the 
purposes of comparison, normalized streamwise velocity pertur- 
bation is shown in Figure 7(b) corresponding to chord Reynolds 
number of 3.67 million. The maximum amplitude reaches a value of 
nearly 15% of the boundary-layer edge velocity but it is instruc- 
tive to note that the flow was still laminar. 


□ Expertmenf 


t_| 

.02 .04 .06 .08 .10 


x/c = 0.25, G’ 5.679 x/c - 0.275, (3 1.6 18 




< a> 


x/c = 0.25, G = 9.563 



ORIGINAX PAGE IS 
OF P(X)R QUALITY, 


( b) 

Figure 7 
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EFFECTS OF SUCTION 


The influence of moderate amounts of additional suction is 
shown in Figures 8<a)and(b). Nominal suction corresponds to the 
minimum suction required to maintain an attached laminar boundary 
layer in the test region. The porous suction strip ends at x/c 
= 0.225 and at this location an increased suction level appears 
to have consistently reduced the streamwise velocity pertur- 
bation. However, at x/c = 0.25, where there is no local sue 
tion, the streamwise velocity perturbation is consistently larger 
than those obtained with nominal suction in the compression re- 
gion. The available results are not sufficient to make major 
conclusions on the effects of suction. A more detailed study is 
necessary. 


x/c-0.225. G-12.108 


— Theory 
Experimenl 
□ Nominal Suction 
A 50% Excess Suction 
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x/c-0.25. 0 = 6.023 
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Figure 8 
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VARIATION IN WAVELENGTH AND MAXIMUM PERTURBATION AMPLITUDE 


The variation in wavelength with Gortler number is shown in 
Figure 9<a>. The Gortler number G^ is based on the conditions 
at the location of the minimum radius of curvature in the concave 
region. Theoretical results were computed using Floryan's code 
(reference 4). A good agreement is observed. Unlike all pre- 
vious experiments, significant variation in wavelength was ob- 
served with changes in freestream parameters, as seen from this 
f igure. 

The variation in normalized maximum disturbance amplitude is 
shown in Figure 9(b). Local suction exists in the region 0. 17 < 
x/c 5 0.225. In this region the growth appears to be gradual but 
is followed by a sharp increase to the maximum at the end of the 
concave region. The perturbation amplitudes decrease rapidly in 
the convex region. Increased levels of suction decrease the 
amplitude locally but appear to trigger a more rapid growth 
downstream. 



0 .2 .4 .6 .8 

Wavelength, X, cm 


< a ) 
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PHASE RELATIONS 


Current theoretical models for velocity perturbations due to 
the presence of Gortler vortices assume the streamwise and normal 
velocity components to be in phase and the spanwise component to 
be out of phase by 90'-. However, cross-spectral analysis of 
experimental data consistently gave results at variance with 
these assumptions (Figure 10(a>>. The streamwise ( u) and the 
normal ( v) velocity components were out of phase with each 
other by 180'-' whereas the spanwise component ( w) was out of 
phase by 90'-' or 270'- with the streamwise component ( 270-' or 90'-' 
with the normal component) depending on the height above the 
surface. The w— component changes sign across the center of the 
vortex layer. Such a behav io r is consistent with the correct 
physical model of Gortler vortices shown in Figure 10(b). Based 
on these results the following correct theoretical model is 
suggested for velocity perturbations' 


u( X, 

Y' 

Z > = U< X, 

y> 

Cos 

az 

exp (/ ^^dx> 

V< X, 

Y' 

Z ) = v< X , 

y> 

Cos 

( az 

■t It) exp ( J P^dx) 

W< X, 

Y' 

Z > = w( X , 

y> 

Cos 

( az 

± l(/2) exp < f ^^dx). 


MEASURED PHASE RELATIONS PHASE RELATIONS FOR VELOCITY PERTURBATIONS 

FOR VELOCITY PERTURBATIONS (Theoretical model) 



Figure 10 
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SUMMARY OF RESULTS 


• Flow visualization and Laser Velocimetry results corroborate 

dependence of vortex wavelength on freestream parameters 

• Measured values agree well with theoretical wavelengths for 

maximum amplification 


• Velocity perturbation profiies determined in the concave and 

convex regions 

• Velocity perturbation profiles and vortex growth trends agree 

well with theory 

• Vortex strength increases in the concave region followed by 

a rapid damping in the convex region 

0 Influence of moderate amounts of suction obtained 

• Phase relationships between velocity components agree well with 


the correct physical model 



FUTURE EFFORTS 


• More detailed boundary layer measurements to better 

define stability parameters 

0 Extend the scope of the experiment to study the influence of 
suction and curvature distribution more thoroughly 

# Study vortex interaction between Tollmlen-Schllchting waves 

and Gdrtler vortices 

0 Transition In the presence of Gortler vortices 
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INTRODUCTION 


Laminar -flow control (LFC) has the potential _ for significant drag 
reduction in a wide range of aerodynamic vehicle applications. The primary 
objective of LFC is the achievement of extended areas of laminar flow by 
delaying boundary-layer transition where the primary drag reduction benefit is 
in the decrease in the skin friction. Two competing methods to apply wall 
suction are the use of discrete slots or a contiguous porous skin along the 
aerodynamic surface as discussed in (Refs. 1 and 2). The porous surface 
provides a nearly uniform continuous suction in contrast to a series of narrow 
slots which provide a discontinuous spatially localized suction distribution. 

The estimation of suction quantities for LFC applications can be 
calculated by assuming continuously applied suction such that the growth 
factors for instability waves are below a prescribed level. Ref. 3 used 
boundary-layer theory and showed that if the same amount of suction is 
distributed at wide porous strips instead of continuous distribution, the 
integrated amplitudes of Tollmlen-Schlichting waves are essentially unchanged. 
However it would be Incorrect to extrapolate this^ result to the case of 
localized suction of a slot because the upstream influence of the slot or 
"sink effect" is not properly accounted for in boundary -layer calculations, 
which are parabolic in nature. Ref. 4 utilized a triple deck theory to 
account for the upstream effect and Ref. 5 provided experimental confirmation 
of these results for wide porous strips. However, the limitations of low wall 
suction velocities, an Inherent assumption of the analytical development, 
raises a question of the accuracy of this procedure for the case of discrete 
suction through a series of multiple slots. 

The first part of the present study concentrates on a comparative 
examination of the growth of instability waves in a Blaslus boundary-layer 
flow by these two methods of wall suction. The second part of the paper 
Includes design criteria for suction hardware to minimize velocity 
perturbations at the entrance of the slot. The slot blockage^ study 
demonstrates that a boundary-layer Instability develops downstream which may 
lead to premature transition. The lower limit of this instability for a slot 
blockage width equivalent to the local boundary-layer displacement thickness 
represents the most restrictive criteria for the suction rate through a given 
slot. 


EXPERIMENTAL CONFIGURATION 


This study was conducted In the Lockheed-Georgia LTVT (Low Turbulence 
Wind Tunnel) which has a 0,6m X 0.9m rectangular test section by 6.5m length. 
The free stream turbulence level is of the order .02-. 04% of the free stream 
velocity. The flat plate depicted in Figure 1 has an elliptical leading edge 
and is constructed of 1.27cm thick aluminum plate with dimensions 0.9m X 3.0m 
length. Chordwise and spanwise static pressure measurements were made with a 
total of 32 static surface ports. The pressure gradient was set close to zero 
over the test region by the adjustment of the tunnel walls with a Cp variation 
of +0.25%. A variable deflection trailing edge flap attached to the 
downstream end of the plate controlled the leading edge stagnation line. 

Three suction locations were chosen on the flat plate at x=l. 22, 1.525 and 
1.83m from the leading edge. The suction power was provided by a radial vaned 
pump attached to a vacuum chamber, where the total mass flow was maintained at 
a constant value by a sonic nozzle at the entrance of the chamber. The flow 
rate from each panel was controlled by a ball socket valve and monitored by a 
calibrated pressure drop across a porous membrane installed in each pressure 
line. To accurately measure the mean flow and the dynamic fluctuations of the 
instability waves, the boundary-layer surveys were made with a computer driven 
three-axis probe survey mechanism of resolution .015mm which was mounted 
outside the test section of the tunnel. The hot wire anemometer used in this 
study was a Disa type 55D01- 
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GEOMETRY OF SLOT AND POROUS STRIP PANEL 


Figure 2 shows the details of the suction chamber design, where a 2-D 
slot panel or a wide porous-strip panel was alternately installed above each 
suction chamber. The porous panels were 2.3cm wide and lOmm thick and 
contained electron beamed holes of diameter 0.25mm spaced in a diamond pattern 
with a porosity value of 20% open. The slot panels were the same thickness 
material with a saw cut of 0.25mm width. A comparison of the slot vs. wide 
porous-strip data was made by replacing the slot panel with the porous panel 
without altering the external boundary conditions of the experiment. The slot 
dimension was 16% of the displacement thickness of the boundary layer as 
compared to the porous strip width of 15.0 displacement thicknesses. The 
0.38cm diameter metering holes were displaced from the centerline of the 
suction chamber by 0.64cm with a pitch of 75wg. The suction chamber had a 
depth of 21wg and a width of 105Wg. 


ELECTRON-BEAM PERFORATIONS OR 
SLOT WIDTH =0.025 cm 



Figure 2 
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COMPARISON OF MEAN BOUNDARY-LAYER PARAMETERS 


From dimensional 
Blasius profile, the 
discrete wall suction 


analysis, with the initial condition specified as a 
change in the mean boundary- layer parameters due to 
can be represented by the functional relationship: 


5*/e = g(yK/S*,w/5*,R5*) 

The first term represents the streamline height of the suction mass flow 
tar upstream of the influence of the wall suction normalized on the 
displacement thickness, the second term is the non-dimensional width of the 
wall suction, and the third term is the local Reynolds number. For low 
suction rates the following expression can be derived for the first term: 

yj^^/ 6* = 1 .88{m)^'^2. where m=wv/U5*. 


data presented in this section Rc* = 2100, which corresponds to a 
5* = 1.65mm and the unit Reynolds number per meter of R 
-1.336 X 10 -The suction flow was maintained at the same value at each 
streamwlse station during the acquisition of the mean and RMS velocity data. 

igure 3a shows the boundary -layer shape factor comparison for slots and 
porous strips for a suction value of m=.035 at successive downstream 
locations. The corresponding displacement thickness variation normalized on 
tne reference value is shown in Figure 3b* 
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COMPARISON OF LOCALIZED Cp AND SHEAR 


The variation of C^ in the near vicinity of the suction region at x-1.22m 
was measured by a small static probe of diameter 0.28*. The comparison of 
slot vs. wide porous strip is shown in Figure 4a for the suction parameter 
m=0.035. The main difference in the two suction methods is the higher peak 
value in C„ of the slot and larger adverse gradient across the slot-suction 
width. The corresponding skin friction data comparison shown in “8“ 
ure 4b was obtained from the mean velocity gradient of the anemometer data 
at a normal distance of 0.15mm from the surface. The significant difference 
is the larger magnitude of the shear near the downstream edge of the slot. 



Figure 4 
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GENERATION OF TOLLMIEN-SCHLICHTING WAVE 


The Introduction of a controlled disturbance in the boundary layer was 
accomplished using the vibrating wire technique. An electromagnet! cal ly 
excited stainless-steel wire of diameter 0. 025mm was positioned 0.76m from the 
leading edge and offset from the plate by 0.2ram spacers and allowed to freely 
vibrate over a length of 0.75m (see Figure 1 previously). A two dimensional 
array of permanent magnets located underneath the surface provided a constant 
magnetic field and a power amplifier with a sinusoidal wave generator provided 
the excitation current for the wire. A weight of 5 Newtons (1.1 lbs) attached 
to the wire provided the tension force, where the natural frequency of the 
taught wire was measured at f=210 Hz. This arrangement allowed a constant 


disturbance amplitude over a spanwlse extent of length 0.5m. The measured 
disturbance wave in the boundary layer is compared to linear stability theory 
in Figure 5a for a dimensionless frequency of F=64.4 X 10“®* Figure 5b shows 


the corresponding mean velocity profiles. 



Figure 5 
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SPECTRAL PLOT OF DELTA FUNCTION VERSUS DISCRETE FREQUENCY 

A cursory examination of the frequency which had the largest amplitude 
growth over the test region was ascertained by the generation of a delta 
function input to the wave generator at a frequency of 5 Hz. The anemometer 
was positioned downstream of the last slot at a position in the boundary layer 
of y/5* =0.6. The spectral plot of the anemometer fluctuations shown in 
Figure 6a demonstrates a wide band of instability waves were amplified. The 
suction test experiment was conducted at the relatively low non-dimensional 
frequency of F=28.3 X 10“®, which corresponds to the maximum growth of the 
discrete frequency of the sine wave input at 127 Hz and is shown in Figure 6b. 
The maximum amplitude of the wave was maintained in the linear range over the 
test section region for the duration of the experiment. 
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COMPARISON OF AMPLITUDE GROWTH OF INSTABILITY WAVE 


The amplitude growth of the instability wave shown in Figure 7 represents 
a comparison of the porous strip versus slot for two suction quantities. The 
data in this figure represent the maximum fluctuation amplitude normalized on 
the maximum amplitude at =1200. To demonstrate the effect of suction 
width on the growth of the instability wave, tape of 0.05mm thickness was used 
to cover 80^ of the width of the porous strip. This configuration had 8 rows 
of holes which were open on the porous strip panel. Figure 7 shows the 
amplitude growth for three cases, namely: 


a. ) wide porous strip of width w/5* =15.0 

b. ) narrow porous strip of width w/8* =2.5 

c. ) and the slot of width w/5* =0.16 


The effects on stability of the narrow porous strip as compared to the 
slot are nearly identical. 



1200 1250 1300 1350 1400 1450 1500 1550 1600 1650 1700 


Figure 7 
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SUCTION HARDWARE DESIGN TO MINIMIZE PERTURBATIONS 


Ref. 6 has demonstrated that premature transition of the boundary layer 
can occur if there is a breakdown of the internal flow in the suction chamber. 

Also, an analytical study (Ref. 7) and an experimental study (Ref. 8) have shown 
nonlinear interactions between three-dimensional and two-dimensional waves 
which excite exponential growth rates which may lead to premature transition. 

The amplitude of the three-dimensional waves in which the nonlinear 
interaction occurs has the low nominal value of u'/U=0.01X. The frequencies 
where resonance occurs are those which have significant linear growth rates 
and the corresponding second harmonics and subharmonics. Therefore, from an 
LFC standpoint it is desirable to minimize the perturbations from the ducting ! 

system, since the generation of Internal finite disturbances will propagate [ 

into the external boundary-layer flow. For this study a closed loop water 
tunnel was used for both flow visualization and the quantitative data reported 
in the subsequent sections. The measurements were obtained with a TSI* 
anemometer using a hot film probe with linearizer, and a Nlcholet FFTf spectrum 
analyzer was used to obtain the frequency data from the velocity fluctuations. 

Figure 8a shows the parameters of the slot flow which were investigated, and 
Figure 8b shows the relationship of the width of the separated region on the i 

parameters Rg and a which are discussed In Ref. 9. | 

*Thermal Systems Incorporated (TSI). 

■*‘fast Fourier transform. 
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LAMINAR BUBBLE TRANSITION IN SLOT 

Flow visualization using dye for a streamline marker demonstrated that 
the flow separates from the leading edge of the slot and reattaches within 4.0 
slot widths depending on the slot Reynolds number and the boundary- layer shear 
parameter (a) as outlined in Ref. 9. For this experiment the suction chamber 
was designed with large dimensions to isolate perturbations from the shear 
layer flow in the slot. Velocity spectra were obtained with a boundary-layer 
probe located 0.25S* downstream of the slot and 0.158* from the surface. The 
velocity spectra shown in Figure 9 demonstrate that transition and breakdown 
of the slot laminar bubble occur at a slot Reynolds number R„ =450. The peak 
spectral data collapses to a dimensionless Strouhal frequency fa/v„ 
=0.030, where the length scale a is the maximum bubble width in the slot. 



f (Hz) 
Figure 9 
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SUCTION CHAMBER OSCILLATIONS 

Flow visualization using dye markers for the streamlines was used to 
observe the flow within the chamber below the slot. The two dimensional flow 
from the slot remains as a coherent jet and impinges on the floor of the «uct. 
The jet then separates into two helical stuctures which are drawn into the 
metering holes. A study was undertaken to examine the effect of varying the 
geometry of the suction chamber on the RMS fluctuations at the mouth of the 
slot with no external boundary-layer flow above the slot. The hot film Pyobe 
was positioned at the entrance of the slot. The metering holes had a pitch 
p/w =50 with the hole Rh< 1200, to minimize pertubations of the flow through 
the^holes. Figure 10 shows the spectral data for three geometries under 
investigation. The maximum amplitude occurs at a Strouhal frequency S fw/Vg 
=•0135 where, at higher frequencies, the amplitudes have an exponential decay- 
The minimum RMS of the oscillations for a given flow rate occurs for a chamber 
depth = 17wg and chamber width = 40wg. 
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METERING HOLE INSTABILITY 


Flow visualization showed that laminar separation of the streamlines 
occurs at the entrance of sharp -edged holes with reattachment within one 
diameter. An idealized experiment was studied where a round tube of diameter 
Z.Bcm and length of two diameters was mounted flush in the center of a flat 
plate. The anemometer was positioned in the center of the exit hole and 
spectral data were taken at various suction rates. For this case the laminar 
separation bubble transitioned at Rj^ZOOO. The data in Figure 11 shows a well 
defined spectral peak at a Strouhal frequency fd/u=0.81. At larger Rj a wider 
spectrum of fluctuations occurs as breakdown of the separation bubble 
encompasses a wider physical dimension of the core flow. Since these data were 
taken in an environment with no turbulence at the entrance flow it is 
recommended that for LFC applications Rj<1500. 



f(Hz) 


ANEMOMETER 



Figure 11 
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SLOT BLOCKAGE INSTABILITY 


Ref. 10 has shown that a three-dimensional boundary-layer instability 
from a single row of suction holes will cause premature transition above a 
given suction rate. Blockage in the slot is analogous to the suction hole 
instability in that the mean boundary layer develops a strong spanwise 
gradient of vorticity adjacent to the obstruction in the slot. The three- 
dimensional Instability growth rate will be dependent on the slot blockage 
width Wk and the value which can be expressed in terms of the slot 
Reynolds number by the relationship Rj^=2.0Rg. For the blockage study, a thin 
rectangular shim was placed flush across the slot width of dimensions Wh 
/Wg=1.3 with an approaching boundary-layer displacement thickness of 8* = 1.5 
Wg and with a free stream Rs*=1300. 

Flow visualization confirmed that the boundary-layer vorticity coalesced 
behind the blocked region and was convected downstream with a discrete 
frequency when the slot Rg>200. To explore this instability further the hot 
film probe was positioned downstream of the blocked region at x/8*=20 at a 
vertical displacement of y/S* =1.0. The suction rate was varied and the 
velocity fluctuations were recorded and are shown in Figure 12. 

The spectral results show a large fundamental frequency with higher 
harmonics which demonstrate the pulse like nature of the Instability for Rg 
>200. For Rg>300 the amplitude of the instability wave is sufficient for 
imminent transition downstream. The fundamental disturbance has a well 
defined Strouhal frequency of fyj^/Uoa=0.036. In most cases this frequency is 
outside the Tollmien-Schlichting neutral stability curve and is therefore 
highly damped. However, at a slot Rg=110 the velocity perturbations are 
significant and may interact with the Tollmien-Schlichting instability. The 
spectral data show that for Rg<75 the slot blockage instability has a 
negligible growth rate. 



Figure 12 
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CONCLUSIONS 


Detailed mean velocity and disturbance amplitude measurements were 
conducted in a Blaslus boundary-layer flow with wall suction applied at three 
downstream locations. The main emphasis of this study was a direct comparison 
of the growth rate of the Instability wave with discrete spanwise slots versus 
wide porous strips. The results demonstrate that the local effects of suction 
through slots or very narrow porous strips have a greater beneficial effect on 
the stabllty of the boundary-layer flow relative to the suction influence of a 
wide porous strip. 

Codes which use continuous suction for the growth rates of the 
instability waves to determine the suction quantities for a multiple series of 
slots will be quite conservative in the estimation of the suction quantity. 
To accurately determine the mean boundary-layer changes from a narrow suction 
strip the "sink” effect of the local suction must be accounted for in the 
neighborhood of the slot. This has been demonstrated in the Navier- Stokes 

analysis of Ref. 11. 

The second part of the paper concerned an experimental study to provide 
guidelines for suction-chamber design and flow rates to minimize internal 
oscillations which propagate into the boundary-layer flow. Based on this 
study, the following observations are summarized: 

1. The flow inside the slot has a characteristic inflectional profile 
due to the presence of a separation bubble, which reattaches within 
four slot widths. The experimental data show that the resulting 
shear-layer transitions at Rg>450. 

2. Oscillations from the suction chamber below the slot exit are caused 
from the two-dimensional jet interactions at the boundaries of the 
suction chamber. The velocity perturbations at the slot entrance can 
be minimized by a suction chamber of depth=17wg and width=40wg. 

3. The flow into the metering hole has a separated region which 
reattaches within one diameter. The perturbations of the entrance 
velocity becomes significant for R^>2000. For LFC applications it is 
recommended that Rj<1500. 

4- Localized slot blockage from any external debris in the environment 
causes a mean spanwise gradient of vorticity in the boundary layer 
downstream of the blocked region. The local three-dimensional 
boundary-layer Instability which results for wu/5* = 1.0 has large 

amplitude growths for Rg>110 with imminent transition at Rg>300. The 
blockage instability has negligible growth rates for Rg<75. 
Therefore, the suction Reynolds number Rg<75 represents the upper 
limit to desensitize the three-dimensional instability for 
intermittent blockage of the slot during flight. 
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LIST OF SYMBOLS 


a 



U 

Uoo 

u' 


u 

V 
V' 

V 
X 

y 

yk 

a 

5* 

n 

0 

p 

T 


b 

d 

k 

o 
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laminar separation bubble length scale inside slot 
defined as (Pq - P„)/l/2p liar 
diameter of metering hole 
frequency 

dimensionless stability parameter 2nfv/Uoo^ 

suction parameter v w/U8* 

pitch of metering holes 

free steam static pressure 

surface static pressure 

unit Reynolds number Uo,/v 

defined as 

Reynolds number Uoox/v 
slot Reynolds number VgWg/y 
hole Reynolds number vd/v 

tfisplacement thickness Reynolds number Uoo8*/v 

Strouhal frequency 

mean velocity in boundary layer 

free stream velocity 

velocity fluctuation 


RMS of velocity fluctuation 

streamlThe velocity upstream of suction influence 

average suction velocity at surface 

velocity fluctuation 

width dimension of wall suction 

streamwise coordinate 

normal coordinate to surface 

streamline coordinate upstream of suction influence 
dimensionless velocity gradient parameter 8U/5y (Wg/Vg) 
displacement thickness 

y(U„yvx)l/2 

momentum thickness 
kinematic viscosity 
density 
shear stress 


SUBSCRIPTS 

blockage 
metering hole 

suction streamline parameter in the boundary layer 

reference value 

slot 

free stream 
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LAMINAR- FLOW CONTROL (LFC) EXPERIMENT OBJECTIVE 

A large chord swept supercritical laminar- flow control (LFC) airfoil has been 
designed, constructed, and tested in the NASA Langley 8-Foot Transonic Pressure 
Tunnel (8-Ft. TPT). The LFC airfoil experiment was established to provide basic 
information concerning the design and compatibility of high-performance supercritical 
airfoils with suction boundary-layer control achieved through discrete fine slots or 
porous surface concepts. It was aimed at validating prediction techniques and 
establishing a technology base for future transport designs and drag reduction. Good 
agreement was obtained between measured and theoretically designed shockless pressure i 
distributions. Suction laminarizatlon was maintained over an extensive supersonic 
zone up to high Reynolds numbers before transition gradually moved forward. Full- i 

chord laminar flow was maintained on the upper and lower surfaces of the slotted 
suction wing at the design Mach number of 0.82 up to a chord Reynolds number of 1 2 x 
10°. When accounting for both the suction and wake drag, the total drag could be 
reduced by at least one-half of that for an equivalent turbulent airfoil. Specific | 

objectives for the LFC experiment are given in figure 1. j 

I 

I 

i 


Conduct basic aerodynamic and fluid dynamics research program 
on a high-performance, swept supercritical, LFC airfoil to determine: 

• Ability to laminarize over extensive supercritical region 

• Ability of stability theories to predict transition and suction 
laminarization requirements 

• Relative merit of slotted and perforated suction surfaces for 
LFC and HLFC 

• Effects of surface conditions and boundary layer influences 
on laminarization 


Figure 1 
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LFC EXPERIMENT DESIGN CONSIDERATIONS 


The Langley LFC airfoil experiment was established to provide basic information 
concerning the design and compatibility of high-performance, swept supercritical 
airfoils with boundary-layer control. Several methods of laminarlzation were 
considered and more than one will be employed during the experiment. The integration 
of laminar-flow control with swept supercr itical wing technology required advanced 
materials and fabrication processes to construct a suitable airfoil. The tests 
required modifications to the wind tunnel to provide low levels of free-stream 
turbulence and acoustic disturbances to ensure that the background disturbance 
environment would not severely affect the maintenance of laminar flow. In addition, 
the necessity for a very large chord airfoil posed a problem that required extensive 
analysis and modification to the 8-Ft. TPT test section to provide a near 
disturbance- free representation of the flight environment. The test section walls 
were contoured with a foam liner to conform to the theoretical streamline shape 
(including wing and wall boundary-layer effects) of the flow about the wing at the 
design condition. The tunnel stagnation chamber was fitted with suitable honeycomb 
and screens to reduce turbulence at transonic speeds, and the test section downstream 
of the model was fitted with an adjustable sonic throat to block diffuser noise 
(figure 2), 


• Design advanced high-performance, large chord, 
swept supercritical LFC airfoil 

• Design and modify 8-FT TPT test section for simulation 

of free air about infinite-span yawed wing at transonic speeds 

- Contoured walls 

• Design and modify 8-FT TPT for flow quality improvement at 
transonic speeds 

- Honeycomb, screens, sonic choke 


Figure 2 
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TEST SETUP FOR LFC EXPERIMENT IN THE 8-FT. TPT 

A schematic of the overall LFC experiment in the Langley 8-Ft. TPT is shown in 
figure 3 along with tunnel modifications. The major component was a large cnord, 23 
swept supercritical LFC airfoil of aspect ratio near one, vhich spanned the full 
tunnel height. Laminar-flow control by boundary-layer removal was achieved by 
suction through closely spaced fine slots extending spanwise on the airfoil 
surface. After passing through the slots, the air passed through metering holes 
located in plenums beneath each slot and was collected by spanwise ducts with nozzles 
located at the ends. From the duct/nozzles . the air passed through airflow system 
evacuation lines, through airflow control boxes which controlled the 3 

suction to each individual duct nozzle, and through sonic nozzles to a 10,000 ft /min 
compressor which supplied the suction. All four walls of the tmnel ^^^oured 

in order to produce a transonic wind tunnel flow which simulated unbounded 
flow about an infinite yawed wing at model design conditions. The contoured liner 
was shaped to conform to computed streamlines around the wing and corrected for 
growth of the wall boundary layer. The success of the LFC experiment depended_ to a 
large extent on the environmental disturbance levels in the test section. Isolation 
of the test section from downstream disturbances was achieved by an adjustable two 
wall choke (sonic throat). Reduction of upstream disturbances such as pressure and 
vorticity fluctuations was achieved by the installation of a honeycomb and five 
screens in the settling chamber. 





I 
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MODEL INSTALLED IN TUNNEL - DOWNSTREAM VIEW 

Figure is a downstream view of the upper surface of the model taken 
immediately upstrecun of the model. This figure illustrates the smooth streamline 
contour of the liner and how it blended with the model. The dark areas at the top 
and bottcm liner model juncture regions are suction panels in the "collar" around the 
ends of the model to control the growth of the boundary layer in these regions. The 
dark area on the left vertical wall and downstream of the model is one of the 
flexible two-wall chokes (sonic throat). The choke plate on the opposite wall is 
hidden behind the model. The dark area immediately in back of the model is the 
tunnel test section access door followed by the downstream high-speed diffuser. 
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Figure 4 
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MODEL INSTALLED IN TUNNEL - UPSTREAM VIEW 


Figure 5 is an upstream view of the finished liner and wing trailing edge as 
seen from the test section diffuser entrance where the liner faired into the original 
tunnel lines. The LFC model extended from floor to ceiling and blended with the 
liner. The offset of the wing mean plane from the tunnel centerline may be seen as 
well as the development of the liner floor and ceiling step which resulted from the 
differential spanwise flow displacement in the tunnel channels "above" and "below" 
the wing surfaces. The dark vertical area on the left of the photograph and 
downstream of model trailing edge is the edge of the test section access door. The 
dark rectangular area ahead of the model is the tunnel contraction region. 



Figure 5 





ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH^ 



458 


MEASURED AND DESIGN PRESSURE DISTRIBUTIONS 


Measured and design chordwise pressure distributions on the upper and lower 
surfaces of the LFC model are shown in figure 6 for two chord Reynolds numbers at the 
design Mach number of 0.82. In general, these representative results indicate 
measured pressure distributions very close to design. Shockfree flow is shown for a 10 
million Reynolds number and nearly shockfree flow for 20 million. The slightly 
overall higher velocities on the upper surface and the chordwise deviation frcrni the 
design pressure distribution were attributed to classical problems associated with 
wind tunnel testing, wall interferences and model deformation under design air 
loads. The velocity field between the upper surface and tunnel wall (supersonic 
bubble zone) was slightly higher than predicted due to the liner contour and 
inability to COTipletely account for boundary-layer displacement effects in the design 
analysis. The irregularity in the upper surface pressure distributions appears to 
be associated with local surface contour deviations, caused by deformation of the 
model under load, and perhaps amplified by the proximity of the contoured wall. As 
chord Reynolds number Increased above 10 million, transition moved forward rapidly on 
the lower surface and the flow could no longer sustain the adverse pressure gradient 
leading into the trailing-edge cusp, and separation occurred at about 80-percent 
chord. This separated flow changed the local effective area distribution of the test 
section resulting in a slightly higher free-stream Mach number and increased upper 
surface shock strength at a chord Reynolds n\jmber of 20 million. 
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MEASURED AND THEORETICAL SUCTION DISTRIBUTIONS 


The measured chordwise suction coefficient (Cq) distribution required to maintain 
full-chord laminar flow over both surfaces at 0.82 Mach number and 10-million chord 
Reynolds number is shown in figure 7 compared to the theoretical suction 
distribution. The required suction level was higher than the theory over most of the 
upper and lower surfaces. Only about one-third of the total suction requirement is 
attributable to the upper surface and the remaining two-thirds to the lower 
surface. The higher suction requirements were due to the upper surface 
overvelocities and the surface pressure irregularities, as well as to higher suction 
required to overcome the cross-flow instabilities associated with the steep pressure 
gradients on the upper and lower surfaces and the minimization of centrifugal Taylor- 
Gortler-type boundary-layer instabilities and interactions in the concave regions of 
the lower surface. The overall higher suction levels are also influenced by tunnel 
disturbance levels which are inherently higher than free-air turbulence levels 
expected in flight. 



Figure 7 




EFFECT OF REYNOLD'S NUMBER ON UPPER-SURFACE TRANSITION PATTERN 


A planform view of the laminar region of the LFC wing upper surface is shown in 
figure 8 with flush-mounted surface thin-film sensors for transition detection 
indicated. The thin-film sensor oscillograph output trace for fully laminar flow is 
essentially flat, while that for fully turbulent flow resembles "white noise" and is 
easily distinguishable. Intermediate cases are assigned a measure based on the 
percentage of time in a local sample that is taken up with turbulent bursts. On 
these diagrams, lines have been drawn to indicate the estimated boundaries of the 
laminar, transitional, and turbulent regions. For a Reynolds number of 10 million 
(at design Mach number), the flow is laminar to the trailing edge over most of the 
span. As the Reynolds number increases to 20 million, the transition zone sweeps 
forward. The bottom end of the wing has more upstream transition at all Reynolds 
numbers, and this is thought to be an endplate effect. 


Moo=0.82 

Rc=10 X 10® Rc=20 X 10® 



Figure 8 



TRANSITION VARIATION WITH REYNOLDS NUllBER 

The data presented in figure 9 show the chordwise extent of laminar flow 
achieved on the upper surface for several Mach numbers up to the design Mach number 
of 0.82, as determined by a grid of flush-mounted surface thin film gages (fig. 8). 

At Rp = 10 million, full-chord laminar flow could be maintained over the upper and 
lower surfaces for all Mach numbers. As Reynolds number was increased for constant 
Mach number, transition moved gradually forward on the upper surface. The Reynolds 
number at which this forward movement began was dependent on Mach number and occurred 
at progressively lower Reynolds numbers as Mach number increased. For the design 
Mach number of 0.82, the forward movement began between 11 and 12 million and reached 
about 65-percent chord at Rq = 20 million. Transition on the lower surface moved 
more rapidly than on the upper surface and occurred near the leading edge for M = 

0.82 and Rp = 20 million. It was concluded that suction lamlnarization over a large 
supersonic zone is feasible to high-chord Reynolds numbers even under non-ideal 
surface conditions on a swept LFC airfoil at high lift. 


oo 



Figure 9 



UPPER- MD LOWER-SURFACE TRANSITION VARIATION WITH Rq 

At a chord Reynolds number of 10 million, the airfoil is fully laminar on both upper 
and lower surfaces (fig. 10). For chord Reynolds numbers above 11 million, 
transition moves forward on both the upper and lower surfaces. Between 11 and 13 
million, the upper surface transition moves Just upstream of 80 percent, and then more 
gradually to about 65 percent at 20 million. On the lower surface, the transition 
line moves from the trailing edge at a Reynolds number of 10 million to about 75 
percent at 1 3 million. At a Reynolds number of 1 5 million, the lower-surface 
transition has moved to 30-percent chord (just beyond the end of the favorable 
pressure gradient) and the flow in the aft cusp is separated. 


Mjo= 0.82 



Figure 1 0 
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EFFECT OF REYNOLDS NUMBER ON DRAG AT DESIGN MACH NUMBER 

The contributions to total section drag coefficient are shown in figure 11 for 
the Mach number of 0.82 over the chord Reynolds number range of 10 to 20 million. 

The assignment of suction drag contributions to the upper or lower surface is 
accountable since the suction drag is computed duct-by-duct. The wake drag is 
separated into upper and lower surface components on the basis of the assumption that 
the wake can be divided between the upper- and lower-surface at the point on the wake 
rake where the stagnation pressure loss is the greatest. The data indicate that the 
larger contribution to the total drag is from the lower surface. The sharp rise in 
wake drag on the lower surface between 14 and 15 million Reynolds number is 
associated with rapid forward movement of transition (fig. 10) and separation of the 
boundary layer in the pressure recovery region of the lower aft cusp. 



Figure 1 1 








MEASURED TOTAL DRAG ON AIRFOILS WITH/WITHOUT SUCTION CONTROL 


A comparison of the measured drag of the present LFC supercritical airfoil 
(identified as SCLFCd )-051 3F) with an equivalent turbulent supercritical airfoil is 
shown in figure 12, The total drag of the swept supercritical LFC airfoil with 
suction slots (solid symbols) includes the suction drag penalty required to maintain 
full-chord laminar flow and thus represents drag levels obtained with the maximum 
extent of laminar flow at the design lift conditions. The open symbols indicate drag 
levels obtained with an unswept supercritical airfoil with fully turbulent (tripped) 
attached flow in the same Langley 8-Ft. TPT. In general, the results indicate that 
about 60-percent drag reduction may be achieved with boundary-layer control over the 
speed range when compared with a turbulent drag level of about 80 counts. 


Data Facility Airfoil 



(L/D) 


max 


• 8-ftTPT SCIFC(1I-0513F 23° 183 

O 8-ftTPT SCTurb 0° 75 
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PRESSURE DISTRIBUTIONS ON SWEPT LAMINAR SUCTION MODELS 

Figure 13 shows the upper-surface pressure distributions at the design lift 
condition (C • .3) on two earlier low-speed LFC airfoils of different sweep as 

Li 

compared to the supercritical LFC design pressure distribution. These low-speed LFC 
designs, based on standard NACA airfoil profiles, have favorable pressure gradients 
over the first 50 percent or more of the chord, and much less severe adverse pressure 
gradients aft than the supercritical design. 


Upper surface, Rq=10 X 10® 
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MEASURED MINIMUM PROFILE DRAG ON SWEPT LAMINAR SUCTION MODELS 

Figure 14 shows only the upper- surface measured minimum profile drag 
coefficients corresponding to the upper-surface pressure distributions of 
figure 13* The breakdown shows that, as might be expected, the supercritical design 
has a somewhat larger suction drag penalty than the NORAIR model, which has the 
greatest extent of favorable pressure gradient (fig. 13 ), and only a slightly larger 
suction drag penalty than the University of Michigan 5~ by 7~Foot Tunnel model while 
the wake drag contribution is about the same in all cases. 


Upper surface, Rc = lO X 10® 
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VARIATION OF TRANSITION FOR SIMULATED HLFC 


A natural follow-on to full chord active laminar-flow control will be a hybrid 
configuration which will combine suction over forward regions of the upper surface 
with natural laminar-flow concepts over rearward regions. An attempt was made to 
simulate such hybrid laminar-flow conditions by turning off suction over rear 
portions of the LFC model while maintaining suction for 0 < x/ c < .25 on the lower 
surface. The simulation involved first establishing laminar flow to the most rearward 
thin film location (about 95 percent chord), then progressively turning off upper 
surface suction starting at the most rearward suction slots and measuring the extent 
of laminar flow downstream of the turn-off point. The variation of transition 
location on the upper surface with chordwise extent of suction is shown in figure 15 
for M = 0.82 and two chord Reynolds numbers. The data Indicate that extensive laminar 
flow could be maintained beyond the point where suction stopped. Clearly, this 
simulated HLFC approach must be considered somewhat "non-ideal istic” and the results 
conservative since no effort was made to seal or smooth the unsucked regions of the 
airfoil. Furthermore, the pressure distribution (fig. 6) is not like that which would 
be designed for a hybrid configuration. 

Accomplishments of the swept LFC airfoil experiment in the 8-foot TFT are 
listed in figure 16, and current and future tests are outlined in figure 17- 


Transition 
location, 
% chord 



L.S. suction, O^x/c^.25 
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SWEPT LFC AIRFOIL EXPERIMENT IN 8-FT TPT 

Accomplishments 


• Achieved shockless pressure distribution 

• Achieved full chord laminar flow on upper and lower surfaces 
for 0.4^M<o^0.82 

• Maintained full chord laminar flow at subcriticai speeds 
and over large supercritical zone at transonic speeds 
for Rc-22 X 10® 

• Achieved about 80% drag reduction for upper surface only 
and 60% for both surfaces for 0.4 ^Moq^0.82 

• Demonstrated feasibility of combined suction laminarization and 
supercritical airfoil technology 


Figure 1 6 


SWEPT LFC AIRFOIL EXPERIMENT IN 8-FT TPT 

Current and Future Tests 

• Complete evaluation of porous-suction upper surface model 
for comparison with geometrically identical slotted model 

• Evaluate hybrid LFC design with reshaped non-suction 
aft upper surface 


Figure 1 7 
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ABSTRACT 

An analytical study of 1 inear— ampl if ying instabilities of a 
laminar boundary layer as found in the experimental data of the 
LaRC/8— Ft. laminar-flow control < LFC) experiment was completed 
and the results are presented. The LFC airfoil used for this 
experiment was a sweptr supercritical design which removed suc- 
tion air through spanwise slots. The amplification of small dis- 
turbances by 1 inear processes on a swept surface such as this can 
be due to either Tol Imien-Schl icht ing ( TS) and/or crossflow ( CF> 
mechanisms. This study consists of the examination of these two 
instabilities by both the commonly used incompressible (SALLY and 
MARIA) analysis and the more involved compressible ( COSAL) analy- 
sis. A wide range of experimental test conditions with varia- 
tions in Mach number# Reynolds number# and suction distributions 
were available for this study. Experimentally determined tran- 
sition locations were found from thin— film techniques and were 
used to correlate the "n-f actors" at transition for the range of 
test cases. 



BACKGROUND/ MOTIVATION 


The design of laminar-flow airfoils is dependent upon the 
use of transition prediction techniques which accurately model 
the characteristics of the boundary layer. A technique which is 
widely used today is the e^ method that was developed sepa- 
rately by Sraith^ and Van Ingen^. This semi-empirical technique 
examines the tendency of a small disturbance wave within the 
laminar boundary layer to amplify over a given surface to the 
point where breakdown to turbulence occurs. The logarithmic 
amplification ratio, or n-f actor as defined in figure 1, is a 
measure of the spatial growth in amplitude of a disturbance wave 
from the initial point of instability. Stability theory assumes 
that the most amplified disturbance (for a critical combination 
of frequency, wave-angle, and wavelength) is the driving mecha- 
nism behind transition. Based on calibrations of a linear 
incompressible theory with low-disturbance wind tunnel data, 
stability analyses have shown fairly consistently that transition 
corresponds to n-f actors on the order of ten. However, since 
stability theory does not account for the influence of the free- 
stream disturbances, it seems fair to assume that n-factor 
calibrations are only valid between data sets with similar 
disturbance environments^. 

This research was conducted for several reasons: to examine 

the amplification of small disturbances by linear processes on a 
modern swept airfoil; to verify the use of a linear stability 
theory as a design tool through transonic speeds; and to compare 
and contrast the incompressible and compressible methods. 


• Accurate transition prediction techniques are essential for 
laminar-flow airfoil design 

• e'’ method developed separately by Smith and Van Ingen 

- Semi-empirical technique 

- Relates stability theory (O-S Eq.) to transition 

- Uses amplification ratio as indicator of transition 

• Solutions to O-S equation give local amplification rates 

n = ln(A/Ao) = (local amplification rate) dx 

Xo 

• items to note: 

- Solutions depend on proper selection of f, and X 

- Theory does not account for initial disturbance levels (Aq) 

• Objectives of analysis 

- Validate codes used in design (SALLY & MARIA) 

- Examine effect of compressibility (COSAL) 

Figure 1 
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BASIC CORRELATION METHOD 


Figure 2 depicts the basic approach employed for this 
research. Measured pressure and suction distributions were used 
to calculate streamwise and crossflow velocity profiles along the 
chord^, A stability analysis of these profiles revealed the 
amount of TS and CF disturbance amplification from the initial 
point of instability. Comparisons of the measured transition 
locations with these calculations determined the correlated 
n-factors at transition. 

Solutions to the governing stability equation, which is the 
celebrated Orr~Sommerf eld (OS) equation, provide the necessary j 

information for analyzing the growth of disturbances within the 
boundary layer. Two of three disturbance variables, frequency 
(f), wavelength (A), or wave— angle < ) , must be specified with 

the third given by the solution to the OS equation. Analysis of 
TS disturbances was accomplished using a f ixed frequehcy and 
wave— angle method. For incompressible flow it was assumed that 
the TS disturbances which amplify most are those moving in the 
direction of the free— stream ( ij) = 0). Thus, the incompressible 
TS (SALLY) analysis consisted of analysis of a range of frequen- 
cies with a wave— angle of zero. Since the two— dimensional dis- 
turbance assumption could not be used for compressible calcu- 
lations, the compressible TS (COSAL)^ analysis consisted of { 

analysis of, first, a range of frequencies with a wave-angle of i 

zero, then a range of wave angles for the most critical fre- 
quency. Suffice to say that the COSAL analysis can be quite time 1 

consuming. The CF analysis was accomplished with the MARIA^ code \ 

which gives approximate solutions using a fixed wavelength and 
frequency method. : 



tBoundary-layeri Stability i 

Experimental ' (CEBECI) '(SALLY, MARIA or COSAL)' Transition 

data I analysis j analysis | n-factors 
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Figure 2 
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EXPERIMENTAL DATA 


The 8-Foot LFC experiment had all the basic information 
necessary for a stability theory correlation - a pressure and 
suction distribution as well as transition measurements®. 

Because the upper surface was designed as ^he primary test 
surface^ only upper surface data were analyzed. A wide range of 
experimental test conditions with variations in free-stream Mach 
number <M^^>, chord Reynolds number <R ), and surface suction ( C„> 
were available for this study. As shown in figure 3 , the avail- 
able data encompassed both the subsonic and transonic regimes. 

For the subsonic cases, the pressure distributions all had a 
strong leading-edge peak as is shown in figure 3 by a represen- 
tative case at M^ = 7. The suction distributions available for 
these cases were all of the type designed for the airfoil with 
suction over the full chord. Only the overall suction levels 
were changed between cases without any attempt to optimize the 
distributions. At moderate to high overall suction levels, full- 
chord laminar flow was achieved for the entire subsonic range of 
M and R„. 

The transonic cases available were all at the design Mach 
number < M^ = .82). As is shown in figure 3, the pressure 
distributions for these cases closely followed the design distri- 
bution except for the slight waviness and slightly higher velo- 
cities towards the nose. The waviness in the pressure distri- 
butions was thought to be due in part to the model deflecting 
under loads. A sonic bubble encompassed most of the upper 
surface and culminated in a weak shock at the higher Reynolds 

numbers ( R„ > 12 million). Full-chord laminar flow was achieved 
c 

until the shock appeared and then transition occurred ahead of the 
shock. Because of this, analysis was done on the "hybr id"-type 
suction distributions in which suction was systematically shut 
off at the trailing edge, moving towards the leading edge. 


Upper surface 


SUBSONIC TRANSONIC 

. 4 . . 6. ,7 .82 




Change distribution 




X/C 


Figure 3 
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SAMPLE SUBSONIC CALCULATIONS 


Figure 4 shows the results of an incompressible TS analysis 
of a sample subsonic case < = .6 and = 10 million). For the 

pressure and suction distributions shown in the figure/ n~factors 
were calculated for a range of frequencies with a wave-angle of 
zero. The combination of the strong adverse pressure gradient 
near the leading edge and the full— chord suction distribution 
confined the TS amplification to within the first 15% chord 
followed by a stable region to the trailing edge. For this case, 
with its relatively low overall suction level, transition was 
measured to be around 10—15% chord. As is shown in figure 4 , 
the TS growths reached a logarithmic amplification ratio of 
n = 10 in this region, for the critical frequency of f = 10 kHz. 


-1.0p 
-.8 - 
-.6 
-.4 - 
-.2 - 

Pressure 
coefficient, 0 

.4 

.6 

.8 
1,0 


.6, Rc = 10 million 
(x/c)tr “• 1 0“* ** Sc 



I 


I 
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Figure 4 
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SAMPLE SUBSONIC CALCULATIONS 


Figure 5 shows the results of an incompressible CF analysis 
of the sample subsonic case. The analysis indicates that CF 
instability is not strong enough to cause transition in this 
case. In fact, over the entire subsonic range of M^ and R 
cases available, CF logarithmic amplification ratios did not ^ 
exceed n = 2.5 in the nose region where the cross-flow instability 
is considered to be most critical. 


Mflo = .6, Rc=10 million 
INCOMPRESSIBLE CROSS FLOW (MARIA) 



Chord station, x/c 


Figure 5 
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RESULTS OF SUBSONIC ANALYSIS 


The results of the subsonic incompressible analysis are 
summarized in figure 6 which shows the calculated transition 
n— factors versus the subsonic range of Mach numbers for the 
Reynolds numbers of 10 and 20 million. Since CF amplification 
was shown to be insignificant# transition was due entirely to TS 
amplification. The figure shows the TS logarithmic amplification 
ratios that correlated with transition. The TS transition n— factors 
varied from n = 8.5 to 10.5 over the range of M^ and R^, ana- 
lyzed. This variation in n-factor is attributed to the fact that 
the wind tunnel disturbance environment is changing with M^ and 


INCOMPRESSIBLE TS 



Figure 6 
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CONCLUSIONS OF SUBSONIC ANALYSIS 


The main conclusions that were drawn from the incompressible 
subsonic analysis^ are listed in figure 7. The transition pro- 
cess was found to be dominated by TS amplification. For these 
subsonic cases/ transition occurred when TS logarithmic amplifi- 
cation ratios in the range of n = 8. 5 to 10.5 was reached near the 
leading edge. 


• CF disturbance amplification at leading edge 

- Small 

- Unaffected by suction 

^Decreasing overall suction levels 

- Decreased amount of laminar flow 

- Allowed higher TS growths at leading edge 

• Therefore, transition process 

- Dominated by TS disturbances 

- Occurred when njs the range 8.5 to 10.5 
was reached 


Figure 7 
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SAMPLE TRANSONIC CALCULATIONS 


Figure 8 shows the results of an incompressible TS analysis 
of a sample transonic case < = .82 and = 20 million). For 

the pressure and suction distribution shown in the figure, 
n-factors were calculated for a range of frequencies with a wave- 
angle of zero. The combination of the wavy pressure distribution 
and suction being turned off after 8% chord allowed strong TS 
amplification to begin. For this "hybrid” case, transition was 
mgssured to occur between 20 and 28% chord. As is shown in the 
figure, TS disturbances grew to a logarithmic amplification ratio 
of n = 10 to 13 in this region, for a critical frequency of 
f = 9 kHz. 


Moo = 82 million 



Suction 

coefficient, 

Cq 


Figure 8 
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SAMPLE TRANSONIC CALCULATIONS 


Figure 9 shows the results of an incompressible CF analysis 
of the sample transonic case. As is shown in the figure^ CF 
disturbances grew quickly to a peak logarithmic amplification 
ratio of n = 4.5 at about 5% cord and then rapidly damped out. 
There was no calculated CF amplification in the region of 
measured transition. For this reason# it was concluded that CF 
disturbances did not contribute to transition. The case shown 
here illustrates the "worst" case calculations in terms of CF 

disturbances over the entire range of M^^ and R . 

^ c 


Mqq= .82 20 million 

INCOMPRESSIBLE CROSS FLOW (MARIA) 
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VARIATION OF TRANSITION FOR SIMULATED HLFC 


During the experiment/ an attempt was made to simulate 
possible "hybrid" laminar-flow control (HLFC) conditions by 
turning off the suction level at various chord locations. The 
HLFC cases consisted of measuring the extent of laminar flow as 
the upper surface suction was systematically shut off# at the 
same time maintaining lower surface suction on only the first 25% 
chord. This simulated HLFC approach is considered somewhat "non- 
ideal istic" since no effort was made to seal or smooth the un- 
sucked regions of the airfoil. Furthermore/ it is clear that the 
current pressure distribution (measured or theory) is not like 
that which would be designed for HLFC. However, these cases 
provide data for which transition occurs far forward of shock 
location thereby being ideal for a stability theory correlation 
at transonic speeds. 

Figure 10 shows the measured variation of transition 
location on the LFC airfoil upper surface with extent of suction 
for M^ = .82 and two chord Reynolds numbers. For Rc = 10 
million/ transition occurred roughly 20% chord past the point 
where suction was shut off, until full-chord laminar flow was 
achieved. For R^ = 20 million, transition occurred roughly 10% 
chord past the point where suction was shut off, until the shock 
present on the rearward end of the model tripped transition. The 
data clearTy indicate the extent of laminar flow that could be 
maintained beyond the point where suction stopped for transonic 
airfoils with weak— shock and shockless conditions. 


SWEPT LFC AIRFOIL SLOTTED UPPER SURFACE 

L.S. suction, O^x/c^.25 



Figure 10 
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TRANSONIC RESULTS 


Figure 11 indicates the results of the incompressible 
stability analysis of the HLFC data that has shockless flow on 
the upper surface ( R^ = lo million). The figure shows the range 
of TS amplification at transition versus the extent of suction 
applied to the airfoil. For each case in which the suction was 
turned off at a different chord location, symbols were plotted 
for the calculated incompressible TS logarithmic amplification 
ratio that corresponded with the measured transition zone. The 
open symbols indicate the n-f actor calculation that corresponded 
to the chord location of the last fully laminar thin-film gage, 
while the closed symbols indicate the calculation that corre- 
sponded to the first fully turbulent gage. The transition zone 
shown on the figure helps illustrate the correlated n-factors 
for the extent of suction used. For this model at test condi- 
tions = .82 and Rc = 10 million, the TS transition 

^~f®ctors for cases with suction past 30% chord are slightly 
lower than those found with little or no suction. The cases with 
the extent of suction less than 30% chord have incompressible TS 
amplification ratios at transition of n = 10. The cases with 
suction past 30% chord have correlated n— factors closer to n = 7. 
Perhaps this indicates that the transition process on swept 
LFC/HLFC airfoils becomes increasingly non-linear as the extent 
of desired laminar izat ion increases. 


Calculated n-factors for variable chordwise suction 
Moo = -82 Rc - 10 X 10® 



Figure 11 
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TRANSONIC RESULTS 


Figure 12 indicates the results of the incompressible 
stability analysis of the HLFC data that has a weak shock formed 
on the rear of the upper surface < F<c =20 million). As is shown 
in figure 11 the range of TS amplification at transition is 
plotted against the extent of suction applied to the airfoil. 

For the cases ranging from no suction to suction over the first 
30% chords the same trend of correlated n— factors of n = 10 was 
found. Suction beyond 30% chord, however, shows that transition 
is increasingly being dominated by non-linear processes as the 
sxtent of suction approaches the shock location. These results 
demonstrate that for a wind tunnel environment, transition can be 
very sensitive to the formation of a shock. 


Calculated n-factors for variable chordwise suction 
Mco = -82 Rq= 20 X 10® 

INCOMPRESSIBLE TS 



Figure 12 
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CONCLUSIONS OF TRANSONIC ANALYSIS 


cone 1 us i ons that were drawn from the incompressible 
analysis of the transonic data are listed in figure 13. For 
cases in which full -chord suction distributions were used, 

Tirahs it ion was sensitive to the rearward attachment point of the 
sonic bubble. If shockless flow was attained < < 12 million) 

then transition could be maintained very close to the trailing 
r^rnidg^r However, if a weak shock appeared < > 12 million) then 
^ trans it ion occurred ahead of the shock for a wide range of over- 
all suction levels which could not be explained with a linear 
stability theory. Turning off suction ahead of shock, which 
s imulated an HLFC concept, allowed transition to occur as a conse- 
quence of disturbance wave amplification and good theory correla- 
tions were found for the available range of R^. The CF disturbances 
were generally small and unaffected by suction, while TS disturb- 
ances grew rapidly to transition once suction was turned off. 

The incompressible n-factor correlations at transition for this 
transonic data were n = 10 when suction was applied up to 30% 
chord and n = 7 when suction was applied further than 30% chord. 


• Full-chord suction cases 

- Transition occurs ahead of shock location for Rc-12 million 
which could not be explained by linear stability theory 

- Transition occurs at trailing edge for shockless flow (Rc'^12 million) 
•'Hybrid* suction cases 

- Turning off suction ahead of shock allowed natural transition and 
good theory correlations for range of Rq 

- CF disturbances small and unaffected by suction 

- TS disturbances grew rapidly to transition after suction turned off 

• Incompressible correlations at transition 

- 10 near leading edge 

- 7 in the mid-chord region 


Figure 13 



SAMPLE TRANSONIC CALCULATIONS 


Figure 14 shows the results of a compressible TS analysis of 
the sample transonic case < = .82 and = 20 million). For 

the pressure and suction distribution shown^ compressible 
n-f actors were calculated first for a range of frequencies with a 
wave angle of zero, then for a range of wave angles for the most 
critical frequency. Since transition was measured between 20 and 
28% chord, the compressible TS n— factor at transitions corre- 
sponds to n = 5 - 7.5 for the most critical disturbance 
<f = 7 kHz and ^ = -50^). 


Mjq-,82, Rc =20 million 
(x/c)tr^ ■20-.28C -»njs = 5-7.5 
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COMPARISONS OF SALLY AND COSAL RESULTS AT TRANSITION 


Comparisons of incompressible and compressible stability 
analysis methods were conducted for a range of M and the 
results are tabulated in figure 15. As shown in the figure, 
compressibility effects reduced the calculated n-factors at 
transition. The difference accounts for nearly 25% reduction of 
the calculated ampl if icat ion for the ~ 0. 4 to over 50% 

reduction at = .82. Both have fairly consistent correlations 

with the incompressible n— factors of 9 to 10 and the compressible 
n-factors of 5 to 6. However, as noted earlier, the COSAL anal- 
ysis involves solutions based on a second unknown ('P) which 
increases the time and expense of each compressible case. 


Moo 

Rc 

njs-SALLY 

n-|-3-COS AL 

.82 

20 

10.0 

4.8 

.70 

14 

10.5 

6.3 

.40 

12 

9.0 

6.8 


• Both have fairly consistent correlation 

• COSAL introduces another unknown - U 

• COSAL is more expensive and time consuming 


Figure 15 



SUMMARY 


A wide range of low-disturbance wind tunnel data with 
variations in M^^, and suction were ana 1 yz edT b^ed 6ti a • - 

linear boundary-layer stability theory. The data# which w^r 
a large chord/ swept, supercritical LFC experiment, ehcomjpTssed 
both the subsonic and transonic flow regimes. The analysis 
showed that as far as 1 aminar izat ion through the control of dis- 
turbance amplification, the model performed as designed. The CF 
amplification was successfully controlled through design by using 
a proper combination of Reynolds number, pressure gradients and 
sweep angle. The TS ampl if icat ion could be successfully controlled 
by suction as long as strong non-linearities (e.g. shocks) were 
not present on the airfoil. As is shown in figure 16, the TS 
logarithmic amplification ratios that were calculated at transi- 
tion were n = 9 to 11 for incompressible theory and n = 4 to 6 
for compressible theory. Based on these results, it seems fair 
to conclude that linear stability theory can be used as a guide 
in the design of laminar-flow airfoils through transonic speeds. 
However, a word of caution is needed. Increasing conservatism in 
the choice of a design n— factor is required with increasing 
lengths of LFC, especially at transonic speeds. 



( 


I 

• For the range of subsonic and transonic cases 
TS disturbances dominated the transition process 

• Transition corresponded to 

- Incompressible values of 9=^njg=^11 

- Compressible values of 5=^n-pg=^6 

• Linear stability theory can be used as a guide in 
the design of laminar flow airfoils at transonic speeds 

• Increasing conservatism in choice of n-factor is required 
with Increasing length of LFC 


Figure 16 
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GROWING BOUNDARY LAYRR 



The boundary layer on a flat plate is illustrated in Figure !• The 
streamwise, normal, and spanwise directions are denoted by x, y, and z, re- 
spectively. The displacement thickness 6* increases in the streamwise 

direction. At any distance xq from the leading edge, one can define a 
Reynolds number Re based on the velocity u and kinematic viscosity v 

in the free stream and the local displacement thickness. 

The transition process of the growing boundary layer is influenced by 
significant non-linear, three-dimensional and non-parallel effects. 
Unfortunately, existing computer resources are only adequate for treating two 
out of these three effects. 
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PARALLEL BOUNDARY LAYER 


A common compromise is to study the parallel boundary layer (see Figure 
2) instead of the true, growing one. Here the focus is on the vicinity of 
some point xq (see Figure 1) and the approximation is that the displacement 
thickness remains constant (in x) at the value 6 • As a consequence, 

only the non-linear and the three-dimensional efPects are taken into 
account* The neglect of the non-parallel effects should be serious only if 
there is appreciable growth on the scale (X^^) of the Tollmlen-Schlichting 

(TS) waves* Since the mean flow in the parallel boundary layer is uniform 
in X, a Fourier approximation in x is highly accurate; moreover, only one 
spatial wavelength needs to be resolved for the temporal transition problem* 
Thus, highly resolved computations can be performed, well into the strongly 
non-linear regime* 








Non-linear 

Three-dimensional 


Figure 2 
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BOUNDARY-LAYER CONTROLS 


Figure 3 depicts three types of laminar-flow control (LFC) techniques. 
The dimensionless parameters describing self-similar solution of the Falkner- 
Skan boundary-layer equations are defined in the figure. These are 3 for 
pressure gradient, F^ for suction, and x for heating. In the last case, 
the free stream and wall temperatures, denoted by T^ and T^, 
respectively, (and given in degrees Kelvin), differ and the kinematic 
viscosity depends upon the temperature T. 


Pressure gradient 


Uoo a 
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Suction 


Heating 
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CONTROLLED, PARALLEL BOUNDARY LAYER 


The essential features of the controlled, parallel water boundary layer 
are summarized in Figure 4, The mean flow is taken to be the solution of the 
appropriate self-similar boundary-layer equations. The temperature dependence 
of the viscosity and conductivity are taken from Ref. 1. (However, the 
specific heat and density are presumed to be independent of temperature; their 
values in the free stream are used throughout the flow field.) The parallel 
flow assumption permits the use of a periodic approximation in x and z. 
Shown below is the form of a single Fourier component. This alone is 
appropriate for the linearized stability problem, but for the full non-linear 
simulations reported in this paper a finite Fourier series is employed in 
both X and z. 


• Parallel flow assumption 

U(X, Y, Z, T) = 0(Y) e' - a>T) 

• Mean flow described by Falkner-Skan equation with pressure gradient, 
suction and/or heating controls 

• Viscosity and conductivity based on empirical formulas for water 

• Reynolds numbers based on displacement thickness and free 
stream conditions 


Figure 4 
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LINEAR THEORY RESULTS FOR HEATED BOUNDARY LAYERS 

The linear stability of heated water boundary layers has been examined by 
Wazzan, Okamura, and Smith [2] and by Lowell and Reshotko [1]. Wazzan et al. 
included the effects of temperature only Insofar as they modified the mean 
flow; their stability analysis admitted only velocity and pressure 
fluctuations. The linear results of Lowell and Reshotko, on the other hand, 
Included the effects of temperature fluctuations. Figure 5, taken from Ref. 
I, compares the neutral stability curves that ensue under heating control 
for T = 60°F when the temperature is allowed to fluctuate (solid lines) 
and when it is held fixed at its mean value (dashed lines). Note that the 
boundary layer is actually more stable than is suggested by linear theory 
calculations which neglect temperature fluctuations. Note also that although 
the flow Initially becomes more stable as the wall is heated, its stability 
eventually degrades with additional wall heating. Lowell and Reshotko have 
performed further calculations in which the effects of the temperature upon 
the density are included. They found that this effect moves the neutral 
stability curves to the left, e.g., for the T„ = 90°F case, the solid curve 
moves one-third of the distance to the dashed curve when density fluctuations 
are admitted. 


HEATING CONTROL 
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INCOMPRESSIBLE HEATED BOUNDARY LAYER 


This paper examines LFC techniques in the non-linear regime via direct 
numerical solution of the time-dependent, Incompressible Navier-Stokes equa- 
tions, with temperature dependent viscosity and conductivity. These equations 
(Figure 6) are augmented with an equation for the time evolution of the tem- 
perature. In order to conform to the parallel flow assumption, forcing terras 
need to be added to the momentum and temperature equations. These are denoted 
by Fy and F0 respectively. See Zang and Hussalni [3, 4] for more details. 


V • u = 0 

Su 1 

■^ + w X u = -VP + —V • (/<Vu) + Fu 

It ^ 

(*^ = V X u 


P = p+(l/2)|up 
ff = (T- T^)/(r„ - T^) 


Figure 6 
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NUMERICAL METHODS 


In this work, numerical methods are needed for three problems. The mean 
flow is calculated from the Falkner-Skan boundary-layer equations via a 
fourth-order finite-difference scheme [5]. The linear eigenmodes (and linear 
stability properties) are computed by a Chebyshev Tau method [6] . The time 
dependent Navler— Stokes equations are solved by a Fourier— Chebyshev colloca- 
tion method [7], The mean flow and the linear eigenfunctions are used for the 
Initial conditions of the direct simulation (Figure 7). 


• Non-linear mean flow 

4th-order compact finite difference scheme 

• Linear modes 

Chebyshev Tau method 

• Navier -Stokes solution 

• Fourier -Chebyshev coilocation in space 

• 3rd-order Adams-Bashforth on expiicit terms 

• Crank-Nicholson on implicit terms 
(verticai diffusion, pressure and continuity) 


Figure 7 


SPLITTING METHOD 


The Navler-Stokes algorithm used a splitting method which is outlined in 
Figure 8. The first (velocity) step accounts for the advection and diffusion 
terras. (Although not indicated below, the temperature equation is also 
integrated in this step.) The second (pressure) step enforces the incompress- 
ibility constraint. The boundary conditions in the velocity step are chosen 
to minimize the slip velocity which is present after the pressure step. More 
details are given by Zang and Hussalnl [7]. The book by Canute, Hussainl, 
Quarteronl, and Zang [8] contains an exhaustive discussion of spectral methods 
for simulations of incompressible flow. 

The results reported in this paper were obtained on the NAS Cray 2. The 
collocation grids ranged from 16 x 48 to 64 x 64 (in x and y) for the 2- 
D cases and from 16 x 48 x 8 to 64 x 64 x 16 (in x, y, and z) for the 3- 
D problems. Typical simulations took several thousand time-steps and covered 
from two to five periods of the primary 2-D TS wave. 


Velocity step (tn t*) 

= u* X w* + V • (a* V M*) 

<n) = «(!, tn) 

«*lbndy =£* 


Pressure step (t* — > f„+i) 

mT* = - V P** 

V-y** = 0 

u** • h = g • h 


Final result tn+l) = ^n+l) 


Figure 8 
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MEAN FLOW AT Re = 8950 


The specific cases chosen for this Investigation of non-linear stability 
are Illustrated in Figure 9, The Reynolds number was 8950 and was chosen in 
part because some linear theory results were reported for this Reynolds number 
in Ref. 1. The free-stream temperature was 293 K. The amounts of pressure 
gradient, suction, and heating were chosen so that the flow was neutrally 
stable. (The linear theory results for the heated case did not Include the 
temperature fluctuations.) The mean flow of both the controlled and 
uncontrolled cases are given in the figure. In the heated case, the viscosity 
varies by 50% and the conductivity by 7%. The specific heat and density 
(whose variations are not accounted for in the non-linear simulation) vary by 
0.1% and 1%, respectively. 



0 .5 1.0 1.5 0 .5 1.0 1.5 

U T 


Figure 9 
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MEAN FLOW CHARACTERISTICS FOR Re = 8950 


The use of the LFC techniques has a dramatic effect upon the actual dis- 
placement thickness of the boundary layer. This is quantified in Figure 10. 
The Reynolds number based on the displacement thickness of the corresponding 
uncontrolled boundary is also listed there. Note that in terms of the growing 
boundary layer, the controlled cases correspond to different positions xq, 
with the pressure gradient case having the greatest distance from the leading 
edge and the uncontrolled case having the least distance. 


CONTROL 

^ Uoo 

*^^Blasius 

NONE 


8,950 

P » 0-55 

0.W8 

16,330 

T = I-IO 

1-3986 

11,090 


Figure 10 
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BOUNDARY-LAYER MODES FOR Re = 8950 


The accuracy of the fully non-linear, time-dependent Navier-Stokes code 
is documented in Figure 11. For each of the linear modes specified in the 
first four columns, a highly accurate temporal eigenvalue u and eigenfunc- 
tion were generated using the Chebyshev Tau code. This eigenfunction was then 
input at very low amplitude into the non-linear code as the inital condi- 
tion. The code was run for two TS periods and the growth rate of the eigen- 
function was measured. This is listed in the last column. Listed next to it 
is the growth rate (Imaginary part of u) produced by the linear stability 
code. This growth rate is effectively zero — the real part of o) is 


roughly 0.04. 

Keeping 

in mind the 

size of 

the real part of 

0 ), it is deal 

from Figure 11 

that the 

non-linear 

code is 

accurate to four or 

five digits. 





Linear 

Computed 

Control 

Mode 

a 

P 

growth rate 

growth rate 

Pressure 

TS2-D 

0. 168 

0.000 

0.000095 

0. 000096 


TS 3-D 

0. 168 

0. 168 

-0. 001012 

-0. 001028 

Suction 

TS2-D 

0. 162 

0.000 

0.000093 

0.000093 


TS 3-D 

0,162 

0. 162 

-0.000968 

-0. 000993 

Heating 

TS2-D 

0. 150 

0.000 

0.000093 

0.000097 


TS 3-D 

0. 150 

0. 150 

-0. 000798 

-0. 000793 


Figure 11 
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INITIAL CONDITIONS 


The initial conditions for the direct numerical simulations are depicted 
in Figure 12* The initial velocity field consists of the appropriate mean 
flow plus a 2-D TS wave and (possibly) a symmetric pair of oblique 3-D waves 
(only one of which is shown). The arrow denotes the streamwise direction and 
the indices and label the streamwise and spanwise Fourier components 
relative to the wavenumf)ers a and 3 of the 3-D TS waves given in Figure 
11* The initial amplitude of the waves is measured by their maximum stream- 
wise velocity relative to the free-stream mean flow velocity. 


Mean flow 2-D TS wave 3-D TS wave 



0 .5 1.0 1.5 

U 


Figure 12 
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GROWTH OF FINITE AMPLITUDE 2-D WAVES 


The first set of results pertains to finite amplitude effects upon 2-D TS 
waves. Results of simulations for which the Initial amplitude of the 3-D wave 
was zero and the initial amplitude of the 2-D was 1/2, 1, 2, and 4% are 
summarized in Figure 13. Suction control yields behavior similar to that for 
pressure gradient control [4] and is therefore not given. Two types of 
heating control simulations were performed: in one case the temperature was 
held fixed at its initial mean value (corresponding to the linear theory study 
of Wazzan et al. [2]), and in the other the temperature fluctuations were 
properly accounted for (as in the work of Lowell and Reshotko [1]). The 
initial conditions for both types of heating simulations, however, were 
identical. 

The strength of the perturbation at any Instant is measured here by the 
kinetic energy E of the (k^,k 2 ) = 0.0) Fourier component of the velocity 
field. This is plotted on a semi-log scale in Figure 13, with the time 
measured in units of the period of the 2-D TS wave; the energy at t = 0 is 
denoted by Eq. Thus, one indication of the Impact of non-linear effects is 
the departure of the curves from a straight (and nearly horizontal) line. 
Another is the failure of the curves to lie on top of each other. 




LD 

3 heating 


0 I 2 3 “1 5 

TIME 


I 

I 


Figure 13 
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NON-LINEAR EFFECTS ON 2-B WAVES 


In all cases, the initial growth of the 2-D TS wave is faster than 
exponential. Note that in the non-linear regime the heated boundary layer 
appears to be less stable, l.e., subjected to faster growth of the 2-D wave, 
when temperature fluctuations are included. This is precisely the opposite of 
how temperature fluctuations affect the linear stability results. It is also 
apparent that non-linear effects lead to a more rapid destabilization of a 
water boundary layer controlled by heat than one controlled by pressure 
gradient . 

The tables below (Figure 14) summarize the departure from linear growth 
for 2-D waves with initial amplitudes of 1/2% and 2%. The numbers give the 
ratios (after 1, 2, and 3 TS periods) of the actual amplitude of the 2-D wave 
to that given by linear theory. These numbers, of course, just reinforce the 
data presented in Figure 13, 


NON-LINEAR/LIHEAR AMPLITUDE RATIO FOR l/2t 2-D HAVES 


TIME 

PRESSURE GRADIENT 

HEATING (FIXED T) 

HEATING 

1 

1-OOA 

1-005 

1-001 

2 

1-013 

1-017 

1-016 

3 

1-022 

NON-LINEAR/LI HEAR AMPLITUDE 

1-030 

RATIO FOR 2X 2-D HAVES 

1-080 

TIME 

PRESSURE GRADIENT 

HEATING (FIXED T) 

HEATING 

1 

1-050 

1-065 

1-189 

2 

1-189 

1-257 

1-605 

3 

1-378 

1-567 

2-387 


Figure 14 
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FUNDAMENTAL VERSUS SUBHARMONIC SECONDARY INSTABILITIES 


The transition to turbulence is a strongly three-dimensional process. 
The next part of this investigation focuses on finite amplitude effects upon 
the secondary Instability of the primary 2-D wave to small, 3-D perturba- 
tions. The secondary instabilities may be categorized as fundamental or 
subharmonic. These are identified in flow visualization of the early three- 
dimensional stage of transition as ordered or staggered arrays of lambda 
vortices [9], and have been explained by weakly nonlinear theory [10]. Figure 
15 sketches the periodic array of vortices associated with these secondary 
instabilities. The distance Lj^ is the length of the primary 2-D TS wave and 
it is given by 2Tr/a. 

The 3-D waves that lead to the fundamental instability are TS waves, 
i.e., solutions to the Orr-Sommerfeld equation, whereas the subharmonic 
Instability arises from the interaction of the 2-D wave with a solution of the 
vertical vortlcity (or Squire) equation with streamwise wavenumber a/2. 
These so-called Squire modes are all linearly stable. 



I 

I 


Figure 15 


506 


'FimRilillilllillimiimilllli ' II I mil mm 1 1 li III"I|I|H(||I 



SECONDARY INSTABILITY FOR PRESSURE GRADIENT CONTROL 


Figures 16, 17, and 18 summarize the results of numerical simulations of 
the secondary instability of controlled boundary layers. In all cases, the 
initial 3-D amplitude was 0.01% and the initial 2-D amplitude varied between 
1/2% and 4%. In these plots, the kinetic energy is shown for the 2-D Fourier 
component (k^,k ) = (1,0) and for the appropriate 3-D component— ( 1 , 1 ) for 
the fundamental instability and (1/2,1) for the subharmonic. The 3-D curves 
are labeled by the amplitude of the 2-D wave for the simulation. 


SECONDARY INSTABILITY at Re = 8950 
Fundamental Subharmonic 



Figure 16 





LOG ( ENERGY] 


SECONDARY INSTABILITY FOR HEATING CONTROL (FIXED T) 

The secondary instabilities have the same general character here that 
they do In uncontrolled boundary layers: they are triggered by 2-D amplitudes 
on the order of 1% or more; their growth rate increases with the 2-D 
amplitude; they grow much faster than the primary wave; and the fundamental 
and subharmonic instabilities have comparable growth rates (Figure 17). 


SECONDARY INSTABILITY at Re = 8950 
Fundamental Subharmonic 



TIME TIME 

Figure 17 
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SECONDARY INSTABILITY FOR HEATING CONTROL 


Among these three cases, the secondary instability is strongest for the 
heated boundary-layer simulation which includes temperature fluctuations* But 
the more rapid growth of the 3-D waves in this case is clearly tied to the 
more rapid growth of the 2-D wave. Hence, the principal effect of the 
temperature fluctuations is upon the 2-D waves (Figure 18). 


SECONDARY INSTABILITY at Re = 8950 
Fundamental Subharmonic 



TIME TIME 

Figure 18 
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FLOW FIELD STRUCTURE 


Figures 19 to 23, taken from the work of Zang and Hussalni [4], illus- 
trate the flow field which develops for a fundamental (K-type) secondary 
instability in a heated boundary layer. The calculations were performed at 
a Reynolds number of 1100 and for wall heating of 2.75%* The Initial ampli- 
tude was 2.7% for the 2-D wave and 0.8% for the 3-D wave# This simulation was 
run for nearly five periods of the 2-D TS wave; a grid of 96 x 162 x 216 (In 
X, y, and z) was employed in the later stages. Figures 20 - 23 display the 
key features at a very advanced stage of the transition process. The flat j 

plate is indicated by the solid surface and the mean flow is from the lower 
right to the upper left. I 

The vortex lines Indicate the presence of a hairpin vortex. The spanwlse 
vorticity displays the strong detached shear layer which forms on top of the 
vortex. The normal velocity contours Indicate the regions in which low speed 
fluid from the wall region is Injected upward towards the free stream. The 
plane in the middle of the spanwlse direction is known as the peak plane. It 
is located at the center of the hairpin vortex. In this plane, the detached j 

shear layer is strongest and the upward normal velocity is most intense. Of | 

special interest here are the contours of the temperature perturbation. They ! 

represent the local departure from the mean temperature. Note that the ^ 

temperature fluctuations are strongly correlated with the normal velocity. i 

There is a temperature increase in the peak plane and an even stronger 
temperature decrease near the wall on the outer regions of the legs of the 
hairpin vortex* The line drawings in figure 19 are from a slightly earlier 
stage and are marked with specific contour levels. 


PEAK PLANE QUANTITIES OF HBL AT RE = 1100 

T = 1.10, c^ = 2.7% £j= 0,W 



Figure 19 


510 


I ( lliUttik IIHWHIIIIII II I i I I IlliilttI mil IlllilllllJItilili IM.JUI.; 



CONCLUSIONS 


ORIGIT^;tlr PAGE" IS 

D£ TOOK quality: 


The principal conclusions from this work are summarized below. For all 
of the LFC techniques examined here, finite-amplitude effects are 
destabilizing, i.e., finite-amplitude 2-D TS waves grow faster than predicted 
by linear theory. We also find, in direct contrast to the results from linear 
theory for low-amplitude waves, that temperature fluctuations exert a further 
de-stabilizing influence on finite amplitude 2-D TS waves. The controlled 
boundary layers are, of course, subject to intense 3-D secondary 
instabilities. The instantaneous growth rates of both the fundamental and 
subharmonic instabilities are strongly tied to the amplitude of the primary 
2-D wave. The principal finite-amplitude effects upon the 3-D secondary 
instabilities occur through the faster growth of the 2-D wave. 



VORTEX LINES AT T = 4 7/8 


Heated boundary layer K-type transition 


Figure 20 
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Figure 21 
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Figure 22 


ORIGINAL PAGE 
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TEMPERATURE PERTURBATION AT T = 4 7/8 


Heated boundary layer K-type transition 


Re = 1100 T /T = 1.0275 e, = 2.7% e, = 0.8% 





Figure 23 
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CONCLUSIONS 


• Finite amplitude effects lead to enhanced 
growth of 2-D waves 

• Temperature effects axe further destabilizing 
for finite amplitude 2-D waves 

• More rapid growth of the primary wave leads 
to more rapid development of 3-D secondary 
instabilities 

• The fundamental and subharmonic instabilities 
have comparable growth rates 

Figure 24 
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THREE-DIMENSIONAL BOUNDARY-LAYER PROGRAM 
LAMINAR FLOW CX)NTROL RESEARCH 
NASA LANGLEY RESEARCH CENTER 


The application of stability theory in Laminar Flow Control 
(LFC) research requires that density and velocity profiles be 
specified throughout the viscous flow field of interest. These 
profile values must be as numerically accurate as possible and 
free of any numerically induced oscillations. To date, a high 
percentage of the boundary-layer solutions for use in three- 
dimensional (3-D) stability analysis programs have been obtained 
using quasi-3-D procedures due to the general unavailability of 
3-D boundary-layer programs. 

Several mature 3-D boundary-layer programs based on finite- 
difference techniques (note: integral approaches are not 

considered in the present paper) currently exist, for example, 
references 1 through 4. Unfortunately, these programs are not 
available for general use by the LFC research community due to 
company proprietary controls. In addition to proprietary 
controls, most of the programs are not structured or adequately 
documented for ease of use; that is, their use would generally 
require direct interaction with and/or assistance from the group 
or individual that developed the program. 

Guidelines for the present research project are presented in 
figure 1 and are as follows: (1) develop an efficient and 

accurate procedure for solving the 3-D boundary-layer equations 
for aerospace configurations; (2) develop an interface program to 
couple selected 3-D inviscid programs that span the subsonic to 
hypersonic Mach number range; and (3) document and release 
software to the LFC research community. 


• Reasons for development 


# No general program exists that Is readily available to 
Ql! groups involved in laminar flow control research 


• No general program exists thot Is coupled direct I 
the most often used inviscid programs over the su 
to hypersonic Mach number range 


y with 
Dsonic 


• Most existing programs are at best second-order accurate 
ond often degenerate to first-order accuracy 


• Guidelines for development 

• User friendly 

• High-order accurocyj O(aZ^), 0(aX, aY)^ 

• Flexible interface 

• Document and release to LFC Research community 


(Finite difference nisnerlcs) 


FIGURE 1 


ORDER OF PRESENTATION 


As a result of time constraints for the present symposium, 
it was necessary to condense the material that would have been 
included in three presentations into a simple co-authored 
paper. Consequently, no attempt will be made to either (1) give 
a review of current boundary-layer literature or (2) compare 
numerical results from the present solution procedure with 
previously published numerical and/or experimental results. The 
material to be presented is outlined in figure 2 and is divided 
into two basic categories. The interaction between these 
^^begories will not be discussed due to time limitations. 

Within the inviscid area the focus of the presentation will 
be on the processing of the inviscid flow field data reguired as 
input to the boundary-layer program. The surface Euler equations 
he defined together with the numerical procedure for solving 
the nonlinear system. Several test cases will be presented where 
results obtained from standard interpolation procedures will be 
compared with results from the surface Euler equations. 

Within the viscous area the focus of the presentation will 
be on the numerical procedure and the establishment of its 
accuracy and numerical efficiency. The boundary-layer equations 
not be presented in equation form; however, their 
mathematical character and the required boundary and initial 
conditions will be discussed. The reader interested in the full 
equations can find them in several publications (see refs. 1 
through 6, for example). 


• Inviscid 

• Interface program 

• Inviscid data treatment 

• Surface Euler equations 

• Boundary conditions 

• Numerical procedure 

• Test cases 

• viscous 

• Governing equations 

• Boundary conditions 

• Initio! conditions 

• Tronsforirotion 

• Solution technique 

• Numerical scheme 

• Matrix structure 

• Test cases 


FIGURE 2 


519 


INTERFACE PROGRAM 


The interface program provides the required linkage between 
the most appropriate inviscid code for the configuration and 
flight Mach number regime and the three-dimensional boundary- 
layer program. The functions of the interface program are 
presented in figure 3. The interface program must be capable of 
(1) generating a sufficiently refined boundary-layer grid; (2) 
calculating the metric coefficients and related geometric 
parameters for this grid; (3) interpolating the pressure field 
from the relatively coarse grid used for the inviscid flow field 
solution to the relatively fine boundary-layer grid with j 

automatic smoothing capability as required; (4) solving the 
surface Euler equations on the wall boundary to obtain the 
inviscid velocity components and their gradxents in the 5' 
coordinate system; and (5) generating all required input files 
for the boundary-layer program. In addition to these primary ! 

functions, the interface program prepares diagnostic plot files , 

and will act as the interactive link for viscous-inviscid | 

interaction calculations. 

Success (measure of acceptance and use of software by LFC I 

research groups) of the boundary-layer program will depend, to a I 

large extent, on the structure of the interface software. 


•General function 

• Provide the required linkage between selected inviscid 
programs that span the subsonic to hypersonic Mach 
number range and the boundary- layer program under 
development 

• Specific functions 

• Generation of boundary- layer grid and related metric 
coefficients 

• Interpolation of inviscid pressure distribution. from 
the relatively coarse inviscid grid to the sufficiently 
fine boundary- layer grid 

• Solution of the surface Euler equations to obtain 
the inviscid velocity distribution over the 
boundary-layer grid 

• Generation of input and output files for the boundary- 
layer program 

• Interface for viscous-inviscid interaction 


FIGURE 3 
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TREATMENT OF INVISCID DATA 


The sensitivity of stability theory to oscillations in the 
viscous flow field requires that the inviscid data and body 
geometry (metric coefficients, etc.) be carefully treated in 
order to avoid the introduction of nonphysical oscillations . 
Alternate approaches to inviscid flow-field data treatment are 
presented in figure 4. A conventional approach has been to use 
spline interpolation to transfer data from the relatively coarse 
inviscid grid to the sufficiently fine boundary-layer grid. 

These interpolated results are then numerically smoothed. The 
interpolation approach is the most straightforward; however, 
large errors can occur in the interpolation and arbitrary 
smoothing of vector quantities. This approach yields results 
that are not consistent; that is, they do not satisfy the 
governing inviscid equations within an acceptable error bound. 

The present approach is to interpolate the pressure field 
from the coarse inviscid grid to the fine boundary-layer grid. 
This pressure distribution P^(5,n)» together with known initial 

and boundary values, is used in the numerical solution of the 
surface Euler equations . The advantage of the approach is that 
exact values of u and v can be enforced on certain boundaries, 
the velocity components (u,v) and their gradients in 
the 5,n plane are smooth, and the numerical results are 
consistent . 


• Conventlonol approach 

• Interpolate u, v, w, P from 
inviscid solution grid to 
boundary- loyer grid 


• Present opprooch 

• Interpolate pressure field 
from Inviscid solution grid 
to boundary- layer grid 


Transform u, v, w to u. v 
on boundary-layer grid 


Surface element 


• Dlsodyontoges 



Numerically solve the surface 
Euler system to generate 
u, V on boundary- layer grid 


Inviscid grid 

Viscous grid 


Advantages 


• Errors from Interpolotlon 
of vector quantities can 
be very large 


• "Exact" values of u, v con 
be enforced on known 
boundaries 


• Difficult to obtain 
smooth distributions 

• Numerical results are 
not consistent 


• Smoothness of solution of 
u, V and their gradients 

• Numerical results are 
consistent 


FIGURE 4 
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SURFACE EULER EQUATIONS 


The surface Euler equations are obtained from the three- 
dimensional boundary-layer equations in the limit as ? The 

system obtained from the ^ and n momentum equations are presented 
in figure 5. The system consists of two nonlinear hyperbolic 
equations for the C and n momentum equations and an algebraic 
relationship for the energy equation (perfect gas). The bj^ 
values are known functions of the grid system. The pressure is 
known from the interpolation procedure. The unknowns are 
p, u, and V, where p, u, and v are the density, ^-velocity 
component, and n -velocity components, respectively. The 
equations are a first order, nonlinear hyperbolic system where 
the local streamlines are characteristic. The system can be 
solved with an explicit march in C while remaining fully implicit 
in n provided the initial conditions and boundary conditions can 
be specified as indicated in figure 5. 


A method for calculating the inviscid velocity vectors from 
a specified pressure field was first reported by Cebeci and Meier 
(ref. 6) for incompressible flow around an ellipsoid at angle of 
attack. Vollmers (ref. 7) and Gleyzes and Cousteix (ref. 8) 
improved the approach by integrating the Euler equations. 

Anderson (ref. 3) described a method using experimental pressure 
distributions on rotating turbine blades using the Euler 
equations. The current procedure extends these concepts to 
complex geometry associated with wings and fuselages using a 
second-order numerical scheme. 


•Equation system obtained from boundary- layer equations 
In limit as «=« 

• Equations 

^-momentum 


p[biiu4^ + + (b 


-momentum 

,3 V 


12 ' 


,u^ f bji^uv + bic;V^)] = 


^13 


16ae 


^17an 


22 " “’23 ^ ‘>24“''/ “25''^“ = ">26fl " ‘>27l^> 


Kharacter Of system ■ ’ 

•First order, nonlinear hyperbolic system 
•Local streamlines are chorocterlstlcs 
•Explicit march In Implicit In n 



FIGURE 5 
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BOUNDARY CONDITIONS FOR FUSELAGE "TYPE" GEOMETRY 


As previously noted in figure 5, initial and boundary 
conditions must be specified for physical flows corresponding to 
aerospace configurations. These conditions are presented in 
figure 6(a) for a fuselage type body. The fuselage type body 
with a plane of symmetry relative to the flow field has well 
defined initial and boundary conditions. The initial conditions 
at the stagnation point in the physical plane are u = v = 0. 

These conditions from the required initial values along 

the 5=5^ line in transformed space. Along both the leeward and 

windward lines of the symmetry plane in physical space one can 

specify the conditions v = 0 and 8u/3ti = 0. These form the 

boundary conditions in the transformed plane for the 

lines n = n and n = n 

o max 


• Physical space (X.Y,Z) • Transformed space (c, n) 



FIGURE 6 (A) 
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BOUNDARY CONDITIONS FDR WING GEOMETRY 


Boundary and initial conditions for wing type geometry are 
somewhat more complex than those for fuselages. A schematic of a 
typical wing element is presented in figure 6(b) where a and A 
denote the angle of attack and leading-edge sweep/ 
respectively. At angle of attack, the leading-edge attachment 
line will be displaced from the geometrical leading edge as shown 
in the sketch. The location of the attachment line is not known 
a priori and must be obtained by an iterative solution procedure 
that positions the line such that u = 0 to within a specified 
error bound. This becomes the 5 = 5 ^ line in the transformed 

plane. The remaining initial condition along the 5 = 5 ^ line is 

obtained from the inviscid flow field results (see ref. 9). The 
boundary conditions along either the root chord or symmetry plane 
chord and a chord near, but not located at the wing tip, must be 
specified. The extrapolation conditions shown in the C/H plane 
in figure 6(b) for the wing tip and root chord regions have been 
used for the results contained in the present material. A more 
detailed discussion of the initial conditions for the leading- 
edge attachment line and boundary condition studies for the wing 
tip and root regions is presented in reference 9. 


• Physical space (X,Y,Z) • Transformed space (c,n ) 



FIGURE 6 (B) 


NUMERICAL PROCEDURE 


As previously mentioned, the surface Euler equations are 
hyperbolic. In order to satisfy the stability requirement, the 
difference stencil must be rotated in relation to the cross flow 
velocity as shown in figure 7 (a). 


• Difference stencil 


n 



FIGURE 7(A) 



DIFFERENCE APPROXIMATIONS 


The point upwind differences are used to approximate 
derivatives with respect to ^ (see fig. 7(b)). For derivatives 
\^xth respect to ri a 5— point stencil is used and rotated to 
satisfy stability where e = 1 for v < 0; e = — 1 for v > 0. 

Nonlinear terms are treated in an iterative sense by lagging 
the nonlinear quantity by one iterative step. The final set of 
discretized equations are "locally bloclc diagonal." 


a f 

g 



f 


i-i.j 


(i-g) 


+ f 


i-2,j 



' ^i-lj ^ 


at 1 

^i-l.j+2^ \ 

A(l+A) f 


e = +1 V < 0 

e = -1 V > 0 


FIGURE 7(B) 



PRESSURE COEFFICIENT CONTOURS 


C miOTNAf^A G S .. T 5 

OF. 'fOeft-QUAfeiTY 


A test case of interest to the Langley Research Center's LFC 
program involves experimental and analytical studies of the 
stability characteristics of a fuselage forebody shown in figure 
8. Figure 8(a) presents contours of the pressure coefficient 
distribution interpolated from a relatively coarse inviscid 
solution onto a fine boundary-layer grid. The initial inviscid 
distribution was obtained using the method of reference 10. 



FIGURE 8(A) 
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INVISCID PLOW-FIELD CONTOURS 
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Tlie pressure contours in figure 8(a) indicate a smooth and 
"well behaved" flow field; however, close examination of the 
cross flow velocity contours in figure 8(b) indicates a physical 
problem not easily solved by traditional boundary-layer 
procedures. The solid line on figure 8(b) indicates the locus of 
a V = 0 contour line. In the region above this line v > 0 
(positive cross flow) while below the line v < (negative cross 
flow) , Consequently, both the leeward and windward lines of the 
plane of symmetry are inflow lines. This presents a severe 
problem for traditional boundary-layer approaches where one 
normally obtains a numerical solution of the reduced set of 
boundary-layer equations along the windward line (generally an 
outflow line). Using this solution as one of the two required 
initial data planes, the solution can be marched in the direction 
of increasing n (direction of positive cross flow velocity) to 
obtain the solution for 5 = ^ ^ 

_ 1 O JllCl X 

surface Euler equations could be numerically integrated 

for E = £-»n <n<Ti because the entire n-line was treated 
^ o max 

implicitly. 



FIGURE 8(B) 
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SURFACE EULER VERSUS STANDARD INTERPOLATION 


Comparisons of interpolated u and v values with those 
obtained from solving the surface Euler equations are presented 
in figure 9. The interpolated values were carefully smoothed to 
get the results presented in figure 9(a) as opposed to no 
smoothing for the velocity field obtained from the surface Euler 
equations. In addition, the interpolated values are not 
consistent. 


Cessna fuselage forebody 
= 0.3; cx = 0" 




FIGURE 9 (A) 
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CROSS FLOW VELOCITY CONTOURS 


The carpet plots of cross flow velocity presented in figure 
9(b) are indicative of oscillations resulting from interpolation 
of the velocity components. Stability theory is critically 
sensitive to oscillations of cross flow velocity; consequently, 
the surface Euler approach should be used to obtain smooth and 
consistent input data for the boundary-layer solution. 


Cessna fuselage forebody 

Moo = 0.3j ot = 0“ 



FIGURE 9(B) 
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COMPARISON OF SURFACE EULER WITH INTERPOLATION 


Numerical results are presented in figure 10 for a wing 
having an NACA-0012 airfoil section. The inviscid flow field was 
obtained using the method of reference 10. Comparisons of the 
interpolated and smoothed streamwise velocity are compared with 
surface Euler results in figure 10(a). The agreement in 
magnitude appears to be satisfactory with the exception of a 
region just downstream of the leading-edge attachment line. The 
agreement appears to be independent of span location. 


NACA-0012; yv = 30°; cx = 8°; = 0.22; AR = 2.1 


(upper surface) 



X 


FIGURE 10(A) 
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COMPARISON OF SURFACE EULER WITH INTERPOLATION 


The stability equations will be integrated forward in space 
beginning at the leading-edge attachment line. A closer 
examination of the results presented in figure 10(a) is presented 
as a derivative with respect to x in figure 10(b). The 
oscillation in 9U/3^ presented in figure 10(b) for the 
interpolated values is not acceptable. It can be seen that the 
surface Euler equations yield smooth values of 30/3 C without 
requiring smoothing. These values are consistent with the 
governing equations. 


NACA-0012; ^ = 30“; cx = 8"j = 0.22; AR = 2.1 

(upper surface) 

A 



0 .025 .050 .075 

X 


FIGURE 10(B) 
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THREE-DIMENSIONAL BOUNDARY-LAYER SYSTEM 


Researchers involved in the numerical aspects of LFC are 
familiar with the three-dimensional boundary-layer equations; 
consequently, the equations are not presented. Readers 
interested in the equations are refered to references 1 through 
5. Certain characteristics of the system are presented in figure 
11 and warrant further discussion. 

The system consists of four nonlinear partial-differential 
equations (PDEs) of mixed type. The system is hyperbolic in 
planes parallel to the boundary surface (C,n planes) and 
parabolic in the direction normal to these planes. 

For perfect gas flow, the system is closed with algebraic 
equations for state and viscosity. The unknowns are u, v, w, T 
and p representing the velocity components in the C,ti,5 
directions, static temperature and density, respectively. 

Pressure is constant normal to the wall boundary. 

The present work focuses on laminar flow; however, the 
software has been structured to model the turbulent terms using 
eddy viscosity and eddy conductivity closure. Careful studies 
will be made in this area at some future date. This would allow 
a uniform approach for boundary-layer calculations to be 
developed for laminar, transitional, and turbulent flows. 


• Equations 

• Continuity 

• Momentum (2; X, Y) 

• Energy 

• State; viscosity 


Nonlinear PDE 
Aigebraic 


• Turbulence 

• Eddy viscosity/eddy conductivity | Not LFC requirement 


FIGURE 11(A) 



THREE-DIMENSIONAL BOUNDARY-LAYER SYSTEM 

Boundary conditions have been treated as generally as 
possible. Wall mass transfer can be specified together with 
either wall temperature or heat transfer distributions. Edge 
values are obtained from the interface software solutions of the f 

surface Euler equations. Initial values can currently be j 

generated for fuselage and wing geometries. The reader 1 

interested in the special equation sets for the initial value I 

planes is referred to reference 5. 

I 

i 


• Edge values 

• Surface Euler system 


• Wing "type" geometry ! 

• Leading edge attachment line | 

• Symmetry plane | 

•Root chord approximation j 

i 

FIGURE 11(B) f 


i 

I 


• Boundary conditions 

• Wall boundary 

• U; V (no slip) 

•W^ = W^ (X,Y) 

•T„ = (X.Y) 

or (X.Y) 

• Initial values 

• Fuselage "type" geometry 

• 3-D stagnation point 

• Symmetry plane 
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NUMERICAL SCHEME 


The difference stencil together with a sketch depicting the 
Raetz influence principle (ref. 5), and the difference 
approximations are presented in figure 12(a). Two separate 
stencils are used for convective derivatives in the 5 ,n-plane. 
The particular stencil utilized at each 5 point depends upon the 
sign of the cross flow velocity. For negative cross flow the 
Krause et al. (ref. 11) zig-zag scheme is used to satisfy the 
zone of influence principle. 


• Difference stencil 
• Positive cross flow 


Negative cross flow 



l.j,k-l/2 
O l.J.k-1 


1 










l.j-1 


Raetz influence principle 
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DIFFERENCE APPROXIMATIONS 


The difference approximations are presented in figure 
12(b). The approximations are second-order accurate in planes 
parallel to the wall and fourth-order accurate in the direction 
normal to the wall boundary. In the planes parallel to the wall, 
a weighted three -point upwind scheme is used to approximate the 
convective partial— derivatives (see ref. 12). 

In the direction normal to the wall boundary, a 2-point 
compact scheme with fourth— order accuracy is used (see refs. 13 
and 14 ) . 


- derivatives 

1 


NUMERICAL SCHEME 

Difference relations 


k 4 


II- 2 = (a^4>^ + + a34>i_2) 

l.j ■' 


►n - derivatives 


an , 


i J 
1 


= (b2^4>j + 
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= (C.ct,,, + CoC|>i)'^'2 +(C3<f>j + 

an ^ i ^ J ^ 1 ■' ^ 1 

1 J J 

ic = derivatives 


(qK - Q - 


V > 0 


V<0 


.k-1 


j-Ac(g^Q + 92Q ). Aj^(g 30 +94^1 ). . 

1 / J 


= 0(ac: 


where 


= Solution vector 
IJ 


FIGURE 12 (B) 
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SOLUTION TECHNIQUE 


Substitution of the difference relations into the 
transformed governing equations yields a set of seven nonlinear, 
coupled difference equations in seven unknowns (see fig. 13). 
Newton's method is used to linearize the system. The linearized 
system is solved using an efficient block-tridiagonal matrix 
inversion based on LU factorization. 

The two momentum equations yield a 4x4 block tridiagonal 
system in F, F', G, G', where F = u/ue, G = v/ve, F' = DF/a^ and 
G' = 3G/3^. The energy equation yields a 2x2 block for H and 
H'. The normal component of velocity, w is updated from the 
finite-difference form of the continuity equation. The system is 
solved in an iterative loop until convergence is obtained within 
a fixed error bound. 


• Difference equations are coupled and nonlinear 

• Newton's method used for linearization 

qP = qP‘1 + 

• The linearized system of equations in <sQ forms a block 
tri-diagonal system of the form 

Where 

|A} Ic} are block matrices of system 

Q is solution vector 

4x4 block tri-diagonal for momentum eqs 

2x2 block tri-diagonal for energy eq 

w updated from difference form of continuity eq 

FIGURE 13 


537 


ACCURACY OF METHOD 


The results of a study of the accuracy of the numerical 
scheme together with the effect of the number of grid points 
normal to the wall boundary for a fixed ^re presented in 

figure 14. The truncation error study for Blasius flow clearly 
indicates a slope of 4 to 1 as compared with 2 to 1 for the 
second-order scheme of reference 12. 

The numerical advantage of the fourth-order method is shown 
for 3-D stagnation point flow where streamwise velocity profiles 
(u/Ug) are presented for a fixed ^or 40, 16, and 8 points 

normal to the wall boundary (note that only a sample of the 
computed points is presented for clarity). Comparisons of the 
results for NPZ = 8 with NPZ = 40 indicate that the correct 
profile values can be predicted with as few as 8 points . This 
comparison represents a speed/storage advantage of 5 to 1. 


I 


Log^Q, 
(error 
of Cj) 


FIGURE 14 


Blasius flow 3-D stagnation point flow 



Log^o^ (no. of points) U/U^ 
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CYLINDER ON FLAT PLATE 


Numerical results for incompressible flow past a cylinder 
mounted normal to a flat plate are presented in figure 15. Total 
and cross flow velocity profiles are presented for 50, 25, and 12 
points normal to the wall as well as compared with the second- 
order theory of reference 12. The present fourth-order method 
P^o*3^ces excellent results for as few as 12 grid points normal to 
the wall boundary for this test problem. 


Total velocity Cross flow velocity 



UT/UT^ V/V^ 


FIGURE 15 
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PROLATE SPHEROID AT ANGLE OF ATTACK 


Numerical results for F'„ and G' for incompressible flow 

Vv 

past a prolate spheroid at 8® angle of attack is presented in 
figure 16. The results are presented for three x locations where 
n varies from 0° on the most windward line of symmetry to 180® on 
the leeward line. Of particular importance is the reverse region 
of cross flow: G'^ < 0; x = 1.52; 60® < n < 180® (note that 

G'w = (9G/95)^; G = v/v^ ) . 


Z 



n degrees n degrees 


FIGURE 16 
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SWEPT WING FLOW 


The final swept wing flow test case presented is for a wing 
having an NACA-0012 airfoil section. A schematic of the geometry 
and test conditions is presented in figure 17(a) together with 
calculated values of the chordwise ( 5 ) and spanwise ( ri ) skin 
friction coefficients. 


Schematic of geometry Skin friction coefficient 



FIGURE 17(A) 
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SWEPT WING flow: VELOCITY PROFTLES 


Calculated total and cross flow velocity profiles are 
presented in figure 17(b). Tbe numerical results were obtained 
for grid point distributions normal to the wall boundary of 12, 
20, and 50 points. The profile calculations are presented at 
chord stations of 0, 50, and 90 percent. The present fourth- 

order method obtained accurate results with as few as 12 points 
normal to the wall boundary. 


Total velocity 
X/C 



40X10'Ii 
r Present 

Tc L scheme 

NPZ = 50 

1-0 20 

□ 12 



Cross flow velocity 

,-4 

X/C = .92 


= .50 


.2, .4 .6 .8 1.0 1.2 
UT/UTE 



FIGURE 17(B) 
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CONCLUDING REMARKS 


In concluding, the current software package (interface 
software; boundary-layer software) is operational and has been 
tested for several 3-D flows (see figure 18). 

The interface program has been found to be a dependable 
approach for developing a user friendly procedure for generating 
the boundary-layer grid and transforming an inviscid solution 
from a relatively coarse grid to a sufficiently fine boundary- 
layer grid. The surface Euler equations used for this procedure 
yield smooth inviscid velocity components and gradients along the 
boundary-layer coordinate lines. These data are consistent with 
the governing equations in satisfying the boundary-layer 
equations in the limit as the distance normal to the wall 
boundary becomes very large. The interface program will 
eventually function as the iterative link between the selected 
inviscid software and the boundary-layer software for 
inviscid/viscous interactions studies. 

The boundary-layer program has been shown to be fourth -order 
accurate in the direction normal to the wall boundary and second- 
order accurate in planes parallel to the boundary. The fourth- 
order accuracy allows accurate calculations with as few as one- 
fifth the number of grid points required for conventional second- 
order schemes . 


• Basic programs have been developed and verified for 
several test cases 

•Surface Euler equation approach yields smooth 
and consistent Inviscid edge data 

• Boundary- layer procedure has been verified as 
computationally efficient with fourth order 
accuracy 


FIGURE 18 
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INTRODUCTION 



Lamlnarlzatlon of the boundary layer on the surface of aircraft wings can be 
accomplished by the use of concepts such as Natural Laminar Flow (NLF), Laminar-Plow 
Control (LFC) which uses suction over the entire surface of the wing, and Hybrid 
Laminar-Flow Control (HLFC) which uses suction only over 15 to 20 percent of the wing 
upper surface near the leading edge. The experimental data obtained at NASA Langley 
for the LFC wing have demonstrated that substantial reductions in wing profile drag 
can be obtained at transonic Mach numbers by the use of suction over the entire sur 
face of the airfoil.^ This drag reduction is primarily the result of the presence of 

laminar boundary— layer flow over a large portion of the airfoil surface. ^ 

Extension regions of laminar flow can also be maintained on wing surfaces with 
the NLF concept, which Involves the appropriate choice of pressure distribution to 
limit the amplification of the disturbances that trigger the transition of the lam- 
inar boundary layer to a turbulent boundary layer. Flight tests at the Dryden Flight 
Research Center on a variable sweep TACT F-111 fighter aircraft with an NLF wing 
glove have shown that laminar flow can be maintained over a large portion of the wing ; 

surface at transonic speeds.^ However, NLF is restricted to flight applications at 

low Reynolds number conditions and to wings with relatively low sweep angles. 

HLFC which combines the features of both LFC and NLF applies suction forward of ; 

the front wing-box spar to prevent the transition of the laminar boundary layer due 
to cross-flow and attachment-line instabilities that occur on swept wings. Laminar ^ 

boundary— layer flow is maintained aft of the front spar by the selection of a pres 
sure distribution which surpresses the growth of the disturbances due to Tollmlen- ? 

Schlichting waves or due to cross flow. Although the extent of laminar flow is less 
for HLFC than for LFC, the conventional wing box structure can be retained and the 
mechanical complexities are not as great. ; 

i 

The design of the HLFC wing sections to minimize the wake drag requires the E 

optimization of the drag due to laminar skin-friction and the turbulent boundary- 
layer separation. The traillng-edge recovery pressure is fixed due to the Kutta ^ 

condition for subsonic and transonic flow and the aft pressure gradient, free-stream s 

Reynolds number, and forward suction levels determine whether or not the turbulent I 

flow will remain attached close to the trailing edge. Examination of existing exper 
imental data on the LFC wing has shown that if turbulent separation occurs upstream 
of the 95-percent chord position on either the upper or lower wing surfaces, a pres- 
sure distribution necessary to maintain laminar flow cannot be realized. Therefore, 
the design of HLFC wings sections necessitates the use of reliable theoretical 
methods which accurately predict the locations of laminar boundary-layer transition 
and turbulent boundary-layer separation. 

Several finite-difference boundary layer, stability, and full Navier-Stokes 
equation solvers are available and have produced very encouraging results. However, 
these methods are not well suited for the routine optimization studies that were 
performed during the design of the HLFC wing section. Several new Integral 
boundary-layer methods, which are applicable to swept wings with varying amounts 
of surface suction, have been developed for the prediction of laminar, transition, 
and separating turbulent boundary layers. These methods have been developed for use 
at either subsonic or supersonic speeds, have small computer execution times, and 
are simple to use. The purpose of this presentation is to briefly outline the 
theoretical equations and assumptions which form the basis of these boundary-layer 
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methods and to present the results of several correlation cases with existing experi- 
mental data* The results of the application of these methods to the design of the 
HLFC wing scheduled to be tested in the Langley 8-Foot Transonic Pressure Tunnel will 
also be presented* 
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THEORETICAL METHODS USED DURING DESIGN OF HLFC WING 


The theoretical methods used during the design and analysis of the HLFC wing are 
listed in figure 1 and are classified as 1) potential/viscous design and analysis, 2) 
boundary-layer analysis, 3) wake and suction drag computations, and 4) a modified 
strip method for finite wings. 

Design of the upper surface of the HLFC wing in the local supersonic region at 
the design transonic free-streara Mach number was performed using the Perturbation 
Method of Characteristics technique.^ Modifications to the HLFC airfoi^ section were 
analyzed at transonic conditions using the Korn/Garabedian alrfojLl code^ and at sub- 
sonic conditions using the Multi- Component (MCARF) airfoil code. 

i 

I 

Analysis of design modifications on the viscous characteristics of the HLFC wing 
were performed using several recently developed Integral boundary^layer methods* 

These methods consist of 1) an integral compressible laminar boundary-layer method j 

for swept wings in the presence of suction at subsonic and supersonic speeds, ^ I 

2) criteria for prediction of laminar boundary-layer separation and reattachraent , ^ I 

3) criteria for prediction of the location transition due to either the amplification 

of Tollmien-Schllchting waves, cross flow, or leading-edge contamination, and 4) new 
Integral separating turbulent boundary layer* j 

Theoretical methods were also developed to determine the effect of design modi- 
fications on the drag characteristics of the HLFC wing. These methods account for | 

both the changes in the wake drag and suction drag as a result of applying suction on 
the upper surface leading-edge region of the HLFC airfoil* s 

A modified strip method was developed during this design study to account for 
the finite and swept wing properties of the HLFC wing* This strip method also j 

accounts for taper and both spanwise and chordwise pressure gradients* i 


# Potential/viscous design and analysis 

• Perturbation method of characteristics for Inverse design at transonic 
speeds 

• BGK IBauer-Garabedian-Korn) for transonic analysis 

• MCARF for subsonic analysis 

# Boundary-layer analysis 

• Integral compressible laminar boundary layer with sweep and suction at 
subsonic thru supersonic speeds 

• Short bubble and reattachment criteria 

• T.S. and C.F. transition criteria 

• Separating turbulent boundary layer method (AIAA 86-1832-CP) 

• Cn , and Cn 

□, wake U, suction 

• Modified strip method 


Figure 1 
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GOVERNING EQUATIONS 


The usual governing equations for compressible hydrodynamic laminar boundary- 
layer flow which consist of the continuity, streamwise-momentura, normal-momentum, and 
cross-flow-momentum equations are presented in figure 2, These equations contain the 
terms consisting of variable physical properties, such as, density, p, and dynamic 
viscosity, y. The values of these physical properties vary across the boundary layer 
as well as along the flow direction, and these variations are non-negliglble for 
boundary-layer flow at transonic speeds* At supersonic speeds, these variations in 
the physical properties of fluid within the boundary layer are quite appreciable* 

This means that there is a strong coupling between solutions of hydrodynamic and 
thermal boundary-layer equations at high transonic and supersonic speeds* 

In order to simplify the governing equations for solution by integral techniques 
while maintaining realistic, computational results for the hydrodynamic and thermal 
boundary layers, Stewartson’s transformations are used. These transformations refor- 
mulate the boundary-layer equations of motions into a transformed plane which is 
independent of the varying physical properties of fluid* The relations between the 
velocities in the transformed and physical planes are also shown in figure 2* 


GOVERNING EQUATIONS AND TRANSFORHATIONS 

GOVERNING EQUATIONS 
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TRANSFORMED BOUNDARY-LAYER EQUATIONS AND BOUNDARY CONDITIONS 


The transformed boundary-'layer equations using Stewartson’s tranf ormatlons are 
presented In figure 3* These tranforraed boundary-layer equations are applicable to 
infinite swept wings. The relation between the pressure gradient in the transformed 
and physical planes and the transformed boundary conditions needed to derive the In- 
tegral equations is also shown. The several groups of physical dimensionless param- 
eters used in the development of the integral compressible boundary-layer method are 
also presented. The subscript ”w" Indicates that the parameter is defined at the 
wall. 
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COMPATIBILITY CONDITIONS AND STREAMWISE MOMENTUM-INTEGRAL EQUATION 

WITH SUCTION 


The compatibility conditions used to derive the integral boundary-layer equa- 
tions are presented in figure 4. The streamwise momentum-integral equation with 
suction included was derived by integrating the transformed boundary-layer equations 
from the wall to the edge of the boundary layer and by making use of Leibnitz ^s rule. 
Furthermore, by making use of the dimensionless parameters defined on figure 3, a set 
of simultaneous equations is derived for the solution for the boundary-layer momentum 
thickness and form factor in the transformed plane in the presence of suction. The 
form of function F(S,Kg) is derived by the curve fit of the exact solution results 
for the Falkner-Skan type flow in the presence of suction. 
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BOUNDARY-LAYER PARAMETERS AND CROSS FLOW MOMENTUM-INTEGRAL EQUATION 


The application of reverse Stewartson’s transformations yields the relationships 
between the parameters in the physical and transformed planes as presented in figure 
5* The transformed cross-flow momentum-integral equation is derived in a manner sim- 
ilar to that for the streamwise momentum-integral equation. The numerical solution 
of this equation gives the computational results in the transformed plane in the form 
of the ratio of the boundary-layer thickness in the cross-flow direction to that in 
the streamwise direction as a function of the dimensionless parameter, S, which is 
the ratio of the distance along the normal section in the streamwise direction to the 
chord of normal section. The cross-flow boundary-layer thickness, 6^,, can be calcu- 
lated by multiplying the ratio 5 /6g, which is obtained by solving the cross-flow 
momentum-integral equation, by the streamwise boundary-layer thickness 6^. The 
dimensionless shape of the velocity profile w/W^ can immediately be calculated from [ 

the knowledge of momentum thickness, and displacement thickness. 
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THERMAL BOUNDARY-LAYER EQUATIONS 


Equations for thermal boundary layers in transformed and physical planes are 
presented in figure 6* As a first approximation, assuming that the external pressure 
gradient is negligible (which is approximately the case for the supersonic wing at 
low lift cruise conditions) and that the Prandtl number is close to unity, an analyt- 
ical expression for the temperature profile can be derived* This expression is then 
modified as shown to account for real pressure gradients and the exact value of 
Prandtl number* The real pressure gradient is accounted for through the polynominal 
expression for the velocity ratio u/U^ where the A^’s are functions of the local 
pressure gradient* Once the temperature profiles are computed, the variations of 
dynamic viscosity, kinematic viscosity, and density across the boundary layer can be 
calculated and their average values determined* These averaged values are then 
utilized to compute dimensionless parameters which are used to determine the transi- 
tion location of the laminar boundary layer* 
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CRITERIA FOR DETERMINING LOCATION OF BOUNDARY-LAYER NEUTRAL INSTABILITY 


The relationship between the dimensionless pressure-gradient parameter and an 
"equivalent” Reynolds number based on momentum thickness for neutral instability is 
presented in figure 7. The curve presented is used for determining the location of 
neutral instability for either the streamwise or cross-flow laminar boundary layer. 
The effects of pressure gradients, Mach number less than 1.3, and suction were 
accounted for during the derivation of this curve from the solution of Orr-Sommerfeld 
equations in conjunction with Stewartson’s transformations. 



Figure 7 
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EFFECT OF FREE-STREAM TURBULENCE ON DIFFERENCE BETWEEN TRANSITION- 
AND NEUTRAL-INSTABILITY REYNOLDS NUMBER FOR FLOW ON A PLAT PLATE 


The effect of free-stream turbulence on the transition location on a flat-plate 
laminar boundary-layer subjected to zero pressure gradient is shown in figure 8, The 
transition data were obtained by several investigators in the range of f ree-stream 
Mach numbers from low subsonic to a supersonic value of 3* The temperature profile 
data used to construct this curve were obtained either from experimental data or cal- 
culated from the previously stated theoretical expression for the temperature profile 
presented in figure 6. The curve shows the difference in "equivalent” Reynolds num- 
ber at the transition location and at the point of neutral stability. The equations 
for the averaged value of kinematic viscosity v and "equivalent" integral thickness 
6* which were used to derive this curve from the experimental data are also 
presented. 


*^tran ~ *^inst 


U 0 & n ? I I T 

= "n* = / 



Figure 8 



BOUNDARY-LAYER TRANSITION CRITERIA 


The criteria used in the compressible laminar boundary-layer method to determine 
the transition location due either to the amplification of Tollmien-Schlichting waves 
or to cross flow is presented in figure 9. The dimensionless expression on the ab- 
scissa contains several implicit and explicit physical parameters that are signifi- 
cant during the transition of laminar boundary layer. In deriving the curves shown 
in this figure for several values of free-streara turbulence intensities, use has been 
made of information presented in figure 8. It was assumed that the effect of free- 
stream turbulence on the transition of laminar boundary layer without a pressure gra 
dlent behaves in a manner similar to a laminar boundary layer with a pressure j 

gradient. 




Figure 9 
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EQUATIONS FROM APPLICATION OF SW,E? THEORY 

A schematic illustrating the application of "Strip Theory" for the prediction of 
laminar boundary-layer transition on a finite swept, tapered, and twisted wing is 
presented in figure 10. The wing is divided into a large number of finite strips 
oriented normal to the leading edge (more accurately normal to the local aerodynamic 
center-of-pressure line). The pressure distribution Cp on the edge as well as on the 
centerline of each strip can be determined either by the available two- or three- 
dimensional theoretical flow codes. The local Mach number distributions in the 
direction parallel and normal to the strip are also presented in this figure. These 
local Mach number distributions, Mj^ ^ and are input into the compressible lam- 

inar boundary-layer method to determine whether or not transition will occur due to 
either Tollmien-Schlichting waves or cross-flow instabilities. 
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CORRELATIVE STUDIES 


The present integral boundary-layer method for predicting transition was used 
extensively during the design of the Hybrid Laminar— Flow Control (HLFC) wing to be 
tested in the Langley 8-Foot Transonic Pressure Wind Tunnel. In order to establish 
the validity of the method prior to design of the HLFC wing, the correlative studies 
listed in figure 11 were performed to compare the computational results by the pre- 
sent theory with available experimental data. 


• Variable sweep wing with NACA 642AOI6 section 


• Phoenix wing 


• NASA LFC wing 


Figure 11 
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EXPERIMENTAL SETUP FOR TRANSITION MEASUREMENTS ON NASA AMES VARIABLE SWEEP WING 


O 

A schematic of the experimental setup used by Boltz in the NASA Ames Research 
Center’s 12-Foot Pressure Tunnel to determine the transition location on an infi- 
nitely swept wing at various sweep angles is shown in figure 12, The sweep angle was 
varied from 0° to 50® by insertion of wedges, and the wing tips were kept parallel to 
the free stream with appropriate wing tip extension. The wing section was an NACA 
642A015 airfoil. The chord of the wing was four feet, and the wing was mounted ver- 
tically on the turntable in a semispan manner. In the unswept position the wing had 
an effective aspect ratio of 5.0. The transition locations were experimentally de- 
termined at the various sweep angles and Reynolds numbers using both flow visualiza- 
tion techniques and small microphones located in the model. 





COMPARISON OF PRESSURE DISTRIBUTION FOR NACA 642A015 WING WITH 30° SWEEP 

The pressure distribution computed by the NASA Multi-Component Airfoil Program 
(MCARF) for the NACA bAoAOlS section at an angle of attack of 1° and sweep angle of 
30° is shown in figure 13. The experimental data shown in this figure were obtained 
by Boltz at the same free-stream conditions and sweep. The comparison shows excel- 
lent agreement between theory and experiment. 


A=30°, a = 1°. M^ = 0.27 



Figure 13 
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COMPARISON OF PRESSURE DISTRIBUTION FOR NACA eAjAOlS WING WITH 50® SWEEP 

The comparison between the theoretically computed pressure distribution and 
experimental data by Boltz for an angle of attack of 0° and a sweep angle of 50® is 
shown in figure 14. The computed pressure distribution agrees quite well with the 
experimental data except in the vicinity of the trailing edge. This discrepancy at 
the trailing edge can be attributed to turbulent separation due to outboard washout 
phenomenon which usually occurs at this relatively high sweep angle* 


A = 50°, a=0°, Mfl0 = 0.27 



Figure lA 
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TRANSITION CORRELATION FOR NACA 642A015 WING WITH 20® SWEEP 

The comparison between predicted and experimentally measured transition loca- 
tions for the NACA 642A015 wing with 20° of sweep is presented in figure 15. The 
theoretical results are shown as a function of free-stream Reynolds number for tran- 
sition locations due to 1) Tollmien-Schlichting waves or laminar separation with 
short bubble turbulent reattachraent and 2) transition due to cross flow. Also shown 
in this figure are the computed x/c locations for Tollmien-Schlichting (T.S.) and 
cross-flow (C.F.) neutral instabilities as a function of free-stream Reynolds number. 
The correlation results suggest the following: 

(1) Predicted transition due to T.S. by the present method agrees quite well 

with experimental data for Reynolds numbers less than 20 million. 

(2) For Reynolds numbers greater than 20 million, the predicted transition due 

to C.F. occurs upstream of that due to T.S., and the predicted locations 

due to C.F. agree well with the experimental data. 

(3) For a sweep of 20° and Reynolds numbers greater than 20 million, the 

transition is due completely to C.F. 


A = 20°.a = 0°, Mco = 0.27 


x/c 

Transition 
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TRANSITION CORRELATION FOR NACA 642A015 WITH THE 40° SWEEP 

The correlation of results between present method and Boltz's experimental data 
for the higher sweep angle of 40° is presented in figure 16* The theoretical data 
presented were computed using the pressure distribution generated by the MCARF pro~ 
gram. The following conclusions can be drawn from the correlations presented in this 
figure: 

(1) Transition locations due to T.S. and C.F. are predicted fairly well. 

(2) For Reynolds numbers greater than 8 million, theoretical computations indi“ 

cate that C.F. triggers the transition more abruptly and earlier than at 
the lower sweep angles. 

(3) There is a danger of leading^edge contamination due to C.F. for Reynolds 

numbers larger than 18 million at a sweep angle of 40°. 


A-40°,a=0°, M«)=0.27 
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transition 
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PLANFORM OF PHOENIX SAILPLANE AND AIRFOIL GEOMETRY 

The planform of the Phoenix sailplane and the geometry of the airfoil at the 
test location on the wing are shown in figure 17, The pressure distributions and 
velocity profiles were measured by Raspet^ at several chordwise locations at the test 
location for several values of lift coefficient. The maximum thickness-to-chord ra- 
tio for the test location airfoil is 0,15 and is located at x/c “ 0,35, The maximum 
camber is located at x/c = 0.70, The aspect ratio of the wing is approximately 18. 


Max thickness 14% Max camber 4.5% 
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PRESSURE DISTRIBUTION CORRELATION ON PHOENIX WING 
Cj^ = 0.76 AND Rj^, => 1.464 x 10® 

The comparison between the theoretical pressure distribution computed by the 
MCARF program and the experimental data of Raspet for the Phoenix wing at a lift 
coefficient of 0.76 and a chord Reynolds number of 1.464 million is shown in figure 
18. The results plotted in this figure show that the computed pressure distribution 
agrees quite well with the measured experimental data. These computed and experi- 
mental pressure distributions were input as boundary conditions to the integral 
boundary-layer method used to compute the transition locations presented 
subsequently. 


Cl=0.76, R|m = 1.464 X 10® 



Figure 18 
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PRESSURE DISTRIBUTION CORRELATION ON PHOENIX WING 
= 1.14 AND Rjj = 1.2 X 10^ 

The comparison of pressure distribution computed by MCARF and the experimental 
data for Phoenix wing at lift coefficient of 1.14 and a chord Reynolds number of 
1.2 million is presented in figure 19. The computed pressure distribution agrees 
fairly well with experimental data except near the trailing edge where turbulent 
boundary-layer separation is present. 


Cl=1.14, Rn=1.2 X 10® 



Figure 19 
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TRANSITION CORRELATION FOR PHOENIX SAILPLANE WING 

The correlation for transition and turbulent separation locations on the test 
location airfoil of the Phoenix wing is presented in figure 20. The results are 
plotted as a function of lift coefficient. These results indicate that the chordwise 
location of turbulent separation is predicted fairly well. However, some discrepancy 
exists between the computed and experimental transition location for lift coeffi- 
cients in the range of 0.6 to 1.1 when the theoretical pressure distributions were 
used as input boundary conditions. In order to examine the effect of using different 
boundary conditions, transition locations were also computed using experimental pres- 
sure distributions as input boundary conditions. As seen in figure 20, the differ- 
ence in the theoretical transition locations using the two different boundary condi- 
tions is not significant. It should be emphasized that the author Raspet did not 
measure the location of transition directly but. Instead, Inferred the location from 
the measured velocity profiles; therefore, there may be some discrepancy in the 
actual measured transition locations. 


x/c 

Turbulent 

separation 


x/c 

transition 
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(1) Computed transition with B.C. of computed pressure distribution 

(2) Computed transition with B.C. of experiment pressure distribution 

(3) Computed turbuient separation with B.C. of computed 



Figure 20 
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CORRELATION OF TRANSITION DATA FOR LFC WING WITH PARTIAL SUCTION 

A » 23® AND == 0*826 

The correlation of the transition locations for the NASA Langley LFC wing with 
partial suction on the upper surface are presented in figure 21* The chordwise 
location of suction was varied from an x/c of 0*0 to 0,5* The transition measure- 
ments were made at a free-stream Mach number of 0.826 and chord Reynolds numbers of 
10 to 20 million. Theoretical computations were performed using boundary conditions 
of experimentally measured suction coefficients as input boundary conditions. The 
predicted transition locations are in good agreement with the measured transition 
locations. The results of the three correlation cases presented have demonstrated 
the validity of the new integral boundary-layer method. 

I 
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x/c 

transition 
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Figure 21 
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DESIGN CONSTRAINTS OF HYBRID LAMINAR-FLOW CONTROL WING 


The design constraints and objectives for the Hybrid Laminar Flow Control (HLFC) 
wing to be tested in the Langley 8-Foot Transonic Pressure Tunnel are illustrated on 
figure 22* The design constraints and objectives that were laid out by Mr. P. J. 
Bobbitt, who also directed the HLFC design and optimization studies, are summarized 
as follows: 


Design Constraints 


(1) The geometry of the lower surface of the HLFC wing must be the same as that of 

the LFC wing. 

(2) The geometry of the upper surface panel #1 of the HLFC wing must be Identical 

to that on the present LFC wing and suction must be applied only through 
panel #1. 

(3) The sonic bubble height for the HLFC wing must be no greater, and, if possi- 

ble, smaller than that for the LFC wing. 

Design Objectives 

(1) The geometric shape of the upper surface panels and must be derived by 

the inverse perturbation method of characteristics so that laminar boundary- 
layer flow is maintained up to x/c = 0.6 on the upper surface of HLFC wing 
for C^ = 0.45 and = 0.82, 

(2) The computed turbulent separation location on the upper surface of the HLFC 

wing must be aft of x/c =• 0.95, 


M 




# Suction applied to panel 1 

• Geometry of panels 2 and 3 altered for: 

• Sonic line height constraint 

• Laminar B.L. flow for x/c ^ 0.6 

• Turbulent separation aft of x/c % 0.95 

Figure 22 
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COMPARISON BETWEEN THEORETICAL AND EXPERIMENTAL PRESSURE DISTRIBUTIONS 

FOR LFC AND HLFC WINGS 

A comparison of pressure distributions between the design HLFC wing and the pre- 
sent LFC wing is presented in figure 23* The application of suction is limited to 
0*025 < x/c < 0*26 on the upper surface of the HLFC wing. The pressure distribution 
for the HLFC wing (shown by the dotted line) constitutes the boundary condition 
necessary to achieve the design constraints and objectives outlined on figure 22* 



x/c 


Figure 23 
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PREDICTED TRANSITION AND TURBULENT SEPARATION LOCATIONS FOR HLFC 

The computational results for the predicted transition and turbulent boundary- 
layer separation locations on the upper surface of the HLFC wing are presented in 
figure 24, Predicted results are shown for transition due to Tollmlen-Schlichtlng 
(T.S*) wave amplification, transition due to cross flow (C,F,), and instability due 
to T.S# and C.F, The chordwlse location of computed turbulent separation is also 
presented* These results are plotted as a function of free-stream Reynolds number 
and for a wing C|^ = 0*44 and free-stream « 0*811, 

As shown in figure 24, transition due to T*S* with short bubble reattachraent 
occurs at x/c = 0,58 for the range of Reynolds numbers shown* However, transition 
due to C*F* takes over at a Reynolds number of 15 million* It is assumed with the 
present theory that the transition occurs on the wing due to whichever phenomena 
appears first* Thus, cross flow essentially determines the transition location at 
the higher Reynolds number. In addition, there is a danger of leading-edge contami- 
nation due to C*F* instability at Reynolds numbers larger than approximately 20 mil- 
lion* The computed turbulent separation location is downstream of x/c = 0*95 for the 
entire range of Reynolds numbers* 


A= 23°, Cl = 0.438. = 0.81 1 

1 


x/c 


Free stream Reynolds no. X 10~® 


Figure 24 
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CONCLUSIONS AND RECOMMENDATIONS 


The results of the design studies of the HLFC wing using the integral boundary-layer 
methods of Goradia lead to the following conclusions and recommendations: 

(1) The perturbation method of characteristics was found to be extremely useful for 
determining the geometric changes needed in the baseline LFC airfoil shape to obtain 
the desired HLFC pressure distribution. 

(2) The compressible laminar boundary-layer method with suction was found to be 
quite accurate in predicting the extent of laminar flow for swept wings in the pre- 
sence of suction. Also with this method the suction requirements can be determined 

for different local velocity profiles such as those for Blaslus, separating laminar, 1 

or asymptotic suction velocity distributions. ! 

(3) The turbulent boundary-layer separation method Includes terms to account for the j 

rapid Increase in the turbulent fluctuations in the flow near separation; these terms j 

greatly improved the accuracy of the prediction of the location of turbulent I 

separation. 

(4) The Integral boundary-layer methods execute very rapidly on the computer making | 

it possible to analyze several hundred configurations in a relatively short period of I 

time. i 

(5) The results of the correlative studies generally showed excellent agreement be- 
tween the theoretical predictions and the experimental data. The results also showed j 

that, for values of Reynolds number and wing sweep of practical interest for commer- ; 

clal and fighter aircraft, cross-flow instabilities were predominant in triggering \ 

transition. * 

I 

(6) The analysis of the final HLFC design showed that less than one count of suction = 

drag coefficient was required in conjunction with the appropriate pressure dlstribu- I 

tlon to achieve laminar flow on the upper surface of the wing to the 60-percent chord 

location. This one count of suction drag results In a corresponding 30 to 40 count 

reduction in the wake drag coefficient which suggests that HLFC is a very lucrative _ 

and promising concept for viscous drag reduction at both transonic and supersonic | 

speeds. - 

■ 

(7) In order to prevent laminar transition due to cross flow at large values of I 

sweep and Reynolds number, it is recommended that both the chordwise and the spanwise ; 
pressure gradients be tailored to minimize the growth of the boundary— layer distur— _ 

bances. With the use of a modified strip theory, arbitrary spanwise pressure gradl- ^ 

ents can be accounted for when using the present Integral boundary-layer methods. | 

The present integral boundary-layer methods have been extended to compute both lami- 
nar and turbulent boundary-layer flow and to predict transition locations at super- | 

sonic Mach numbers. The effects of vacations in temperature profile and physical = 

flow properties across the boundary layer, suction, wing sweep, wing taper, and wing | 

twist are also accounted for in these methods. These methods are not limited to , | 

wings alone, but can also be used for the analysis of fighter aircraft fuselage with = 

suction at supersonic flight conditions. Correlation studies are currently under way ^ 

to determine the validity of these extended Integral boundary-layer methods. 
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Results are presented from a numerical simulation of transition control in plane 
channel and boundary layer flows. Details of the channel flow control are available 
in reference 1. The analysis is based on a pseudo-spectral/finite-difference semi- 
implicit solution procedure (ref. 2) employed to numerically integrate the time- 
dependent, three-dimensional, incompressible Navier-Stokes equations in a doubly 
periodic domain. In the channel flow, we find the active periodic suction/blowing 
method to be effective in controlling strongly three-dimensional disturbances. In 
the boundary layer, our preliminary analysis indicated that in the early stages, 
passive control by suction is as effective as active control to suppress instabili- 
ties. Our current work is focused on a detailed comparison of active and passive 
control by suction/blowing in the boundary layer. 


1 


! 

I 

■ 
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GOVERNING EQUATIONS 


CONTINUITY EQUATION 


MOMENTUM EQUATION 



u . 


3U. 


J 3X. 




1 ^ ^ ^ 
p ax. ax. ax. 

I J J 


THERMAL -ENERGY EQUATION 

il + IL = A- 

at ^ ax. PC ax .ax . 

1 r vr J 

CONSTANT PROPERTIES; NO VISCOSITYj TEMPERATURE FIELD UNCOUPLED 
EQUATIONS NOMDIMENSIONALIZED BY U^, h 

FLOW DRIVEN BY A CONSTANT MEAN PRESSURE GRADIENT 2/Re, Re = U^h/v 
CONVECTIVE TERMS PUT INTO A FORM THAT CONSERVES ENERGY AND MOMENTUM 

!!i + u !!!i) = _ ^ + -2 , + 1_ 

at j 'ax. ax.' ax. Re il Re ax. ax. 

J 1 1 J J 

p* . p/e t I ujuj 
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SOLUTION PROCEDURE 


SAME TECHNIQUE AS THE VELOCITY FIELD 

i 

Non-dimensional ized with (T-T )/(T -T ); h and U 

u w U U ! 

Adams-Bashforth 2-step method for the advective terms 
Crank-Nicholson implicit scheme on the diffusive terms. 

Periodicity along allows and i 

I 

Two-D Fourier transform in the x,-x^ plane j 

i 

• The Pseudo spectral method in the x^j^x^ directions I 

! 

Finite differences with variable mesh along the x^ direction j 

Solution in Fourier space as a tri diagonal system ] 

Back transformed into physical space to obtain temperature | 

fi el d at (n+1 ) * 


IMPLEMENTATION OF BOUNDARY CONDITIONS 

No-slip B.C. 

ui = U2 = U3 = 0 

p from x^-momentum equation 
Suction B.C. 

Flow homogeneous along x^ , 

Incoming mass flow rate must equal to the outgoing mass flow rate 

3<u„> 

-- - = 0 
3X2 

<U 2 > = CONST, or g(x^ , x^) 

Hence velocity magnitude and direction at one wall must be preserved 
throughout the flow field to satisfy continuity 

Physically plausible condition is suction-blowing or periodic b. cond. 
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TRANSITION IN WALL-BOUNDE D FLQW? 


2-D Tollmien-Schl ichting waves. 

Formation of streamwise vortices. 

Formation of shear layers away from the wall due to 
vortici ty- induced velocity. 

Secondary instability 
(kinks and spikes). 

Breakdown into smaller scales, formation of wall shear, 
hairpin eddies. 

Turbulent spot - horseshoe vortex-turbulence. 


MODEL PROBLEMS 


Periodic plane channel flow 


Periodic boundary layer. 






MAXIMUM PLANE-AVERAGED RMS VELOCITIES 

The temporal development of plane- averaged maximum velocities is presented. 
These velocities provide comparisons between the controlled flow (three-dimensional 
control) and no-control cases for the channel flow. 




No control 


With control 

a . 

RMS 

d. 

RMS 

b. 

RMS 

e. 

U3 RMS 

c. 

Ug RMS 

f. 

U2 RMS 


c 
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TEMPORAL DEVELOPMENT OF U-COMPONENT FLUCTUATING VELOCITY AND ITS HARMONICS 

The evolution of the various Fourier modes indicated that all amplitudes are 
significantly reduced, and after the third control wave they all decay rapidly. 


No control 


With control 






NORMALIZED ONE-DIMENSIONALIZED WAVE SPECTRUM 


In the uncontrolled case we observe that energy transfer to the high wave num- 
bers is indicated by a full spectrum. In the controlled case, this does not occur 
and energy is concentrated in the low wave numbers preventing the development of 
higher harmonics. 



c|» = <>{ 1 ) 

a. Without control b. With control 
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MAXIMUM U-RMS ALONG X3 (SPANWISE DISTRIBUTION) 

In the uncontrolled flow, peak- valley splitting develops. The control wave 
does not prevent peak- valley splitting, but reduces the amplitudes. The uncontrolled 
and controlled distributions remain in phase. 


c nAxinun ui* Rns alomg X3 at t=20 


G fiAXimjn ui’ Rns along x3 at t=30. 


MAxinun ur rhs along x3 at t=40. 


o NAXINUN Ul* RNS ALONG X3 AT 1=60 




With control 


b. With control 



















U^-FLUCTUATIONS ALONG 

In the uncontrolled case a strong negative spike develops between T = 40 and 
T = 50* No evidence of spike formation and nonlinear distortions is observed in the 
uncontrolled case. As T = 60, the controlled distribution is nearly sinusoidal, 
whereas the uncontrolled case shows a broad frequency content. 


o Ul‘ FLUCTUATIOMS ALOIHG XI AT T = 20 . 


o U1 FLUCTUATIOMS ALONG XI AT T=50 


FLUCTUATION ALONG XI T=40 


With control 


With control 











SURFACES OF CONSTANT TEMPERATURE 


In these figures, three-dimensional representations of the temperature field 
(treated as a passive scalar) are displayed. The uncontrolled flow displays evidence 
of strong mixing and a highly convoluted temperature surface, while the controlled 
flow is relatively uniform and indicates local laminar ization. 



(-)47T (-)4tt 


589 








CONCLUDING REMARKS 


* 2D and 3D wave interaction in channel flow. 

* Preliminary calculations indicate comparable effects of passive and active 
control in the periodic boundary layer. 

* Passive temperature effective to tag flow dynamics. 

* Need for space-evolving numerical experiments. 
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ABSTRACT 


The growth and decay of a wave packet convecting in a boundary layer over a 
concave -convex surface and its active control by localized surface heating are 
studied numerically using direct computations of the Navier-Stokes equations. The 
resulting sound radiations are computed using linearized Euler equations with the 
pressure from the Navier-Stokes solution as a time-dependent boundary condition. It 
is shown that on the concave portion the amplitude of the wave packet increases and 
its bandwidth broadens while on the convex portion some of the components in the 
packet are stabilized. The pressure field decays exponentially away from the surface 
and then algebraically, exhibiting a decay characteristic of acoustic waves in two 
dimensions. The far-field acoustic behavior exhibits a super-directivity type of 
behavior with a beaming downstream. Active control by surface heating is shown to 
reduce the growth of the wave packet but have little effect on acoustic far field 
behavior for the cases considered. 

Active control by sound emanating from the surface of an airfoil in the vicinity 
of the leading edge is experimentally investigated. The purpose is to control the 
separated region at high angles of attack. The results show that injection of sound 
at shedding frequency of the flow is effective in an increase of lift and reduction 
of drag. 
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OBJECTIVES 


The objective is to study the growth and decay of a wave packet developing over 
a concave-convex surface and the resulting acoustic radiations. An engineering 
application is the flow through a wind tunnel contraction. The concave portion is 
known to be potentially unstable to upstream disturbances while the convex portion is 
stabilizing. This instability can be reduced by active surface heating. The 
intensity of resulting sound may still be sufficient to trigger instability 
downstream, for example, on a model leading edge. 

The experimental part is concentrated on separation control by sound emanating 
from the surface of an airfoil at high angle of attack. This reduces the extent of 
the separated region, thus improving performance. 


OBJECTIVES 

• Drag reduction by active surface heating 

• Separation controi by sound emanating from the surface 


TOPICS OF DISCUSSION 


■ Numerical 

• Study of instability over a concave-convex surface 

• Active control by surface heating 

• Resulting sound radiation 

■ Experimental 

• Control of flow separation over an airfoil 
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APPLICATION 


The geometry considered is typical of a wind tunnel contraction. The 
disturbances entering a wall boundary layer in a contraction may grow due to the 
unfavorable pressure gradient, a source of sound radiation, resulting in both 
standing and downstream propagating waves. In addition, a patch of turbulence (or 
vorticity) convecting through the pressure gradient created by the contraction 
geometry undergoes a variable rate of strain, thereby generating sound. The 
perturbed field from the boundary layer as well as from the sources embedded in the 
inviscid flow interact with the leading edge of the model is a possible source of 
transition. 

This presentation includes a study over only one wall of the contraction with 
disturbances externally imposed in the boundary layer. 


WIND TUNNEL CONTRACTION 
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COMPUTATIONAL MODEL 


The computations are carried out using two-dimensional laminar, compressible 
Navier-Stokes equations written in conservation form. The finite-difference scheme 
used is an explicit predictor-corrector method which is second-order accurate in time 
and fourth order on convective terms in space. Such an accuracy is needed for 
reducing numerical dispersion in wave propogatlon problems. 

A steady state is first obtained starting with a boundary-layer solution. An 
unsteady disturbance in the form of the Orr-Sommerfeld solution is then specified at 
inflow and an unsteady solution is obtained. 


Q = Qo+Q' 


• 2-D Navier-Stokes equations: Wt + Fx+Gy=0 

W = (p, pu, pv, E)' 

• Predictor-Corrector F & G are standard flux functions 

• Explicit 

• 4th Order accurate on convective terms in space 

• 2nd Order accurate in time 

• Operator splitting 

• Qq is steady state 

Q' is Orr-Sommerfeld solution 
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BOUNDARY CONDITIONS 


For a complete description, certain boundary conditions must be imposed. At 
inflow, characteristic boundary conditions are used. Three incoming characteristic 
variables are obtained from known solutions. In addition, one outgoing 
characteristic is extrapolated from the interior. 

At the top and downstream, radiation boundary conditions are used. This 
accelerates convergence to steady state and facilitates the upper boundary to be 
brought closer to the surface. 

At the surface, velocities and temperature are specified while the pressure is 
obtained from normal momentum equation. 


• Inflow 

■ Specify three incoming characteristics 

■ p+pcu 

■ V 

2 

■p-pc^ 

■ Extrapolate outgoing characteristic 

■p-pcu 

•At top and outer boundaries 

■ Radiation condition 

■ Pt - pcut = 0 downstream 

■ Pj - pcv^ = 0 top 

• At the wall specify u,v and T 
Pressure from normal momentum equation 
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MEAN FLOW 


The surface has portions of both concave and convex curvatures. The algebraic 
equation describing the curve is a seventh-degree polynomial with first three 
derivatives zero at both ends. This curved surface is accompanied by flat portions 
both up and downstream. The curve rises to 0.02 ft. in a 1 .2-ft. distance. The 
total length in the stneamwise direction is 3.6 ft. 

The results presented are for an inflow Mach number of 0.4 and a unit Reynolds 
number, based on freestream values, of 3.0 x 10^, The mean Mach number and pressure 
distribution are shown. The pressure gradient induced by the surface curvature 
decelerates the flow at first and then accelerates to a maximum Mach number of about 
0.41 . 


MEAN MACH NUMBER AND PRESSURE DISTRIBUTION 


Me 
.41 

.40 

.39 

.38 
Cp 
-.2 

-.1 

0 

.1 

y. 

.02 

0 
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UNSTEADY FLOW - GROWTH OF DISTURBANCES 


The unsteady flow is obtained by introducing disturbances (in the form of a wave 
packet obtained from the Orr-Sommerfeld solution) on the mean flow. The growth of 
the disturbances is shown in terras of the root mean square (rms) of the massflux 
integrated across the boundary layer. The value at each station is normalized by the 
inflow value. The disturbance is Initially amplified on the flat portion. 

Additional amplification comes on the concave curvature due to the unfavorable 
pressure gradient. Farther downstream, the favorable pressure gradient results in a 
reduction of growth, followed by an exponential-like growth on the flat portion. 



I 

E 

I 
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UNSTEADY FLOW - INSTANTANEOUS TOTAL VORTICITY 


The instantaneous total vorticity is shown at successive times as the wave 
packet passes over the surface. The same range of vorticity is plotted for all the 
frames. The figure distinctively shows that as the disturbance passes over the 
concave portion, it changes in amplitude and broadens in bandwidth. Further, on the 
convex portion, drastic changes occur in the bandwidth showing that the favorable 
pressure gradient stabilizes certain components, as if the flow has passed through a 
bandpass filter. The phenomenon is strongly dependent upon the geometry of the 
curvature inducing the pressure gradient. 


y, ft 



0 .9 1.8 2.7 3.6 


X, ft 
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UNSTEADY FLOW - MASS FLUX ALONG THE SURFACE 

Further details of the unsteady flow are shown in terms of the mass flux 
variation with time at several locations along the surface for a fixed distance away 
from the wall in the boundary layer. It is clear that the initial packet grows both 
in amplitude and bandwidth over the initial portion and reaches a maximum at x = 1.0 
feet, just downstream of the concave portion, thus showing a strong non-linear 
behavior. This is followed by a significant decay in the favorable pressure gradient 
portion of the curvature and a growth again on the flat portion downstream. 


y/faj = .062, L = 1 


pu(t) pu(t) pu(t) pu(t) 



01 234 01 234 1 2345 34567 

tUoo/L t Uoo /L tUoo/L t U eo /L 
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SPACE-TIME GROWTH OF PERTURBATION PRESSURE 


This figure shows the space-time variation or the perturbation pressure at a 
fixed distance away from the wall. The additional information provided by this 
picture is that the crests of different waves reach their maximum at different 
positions on the surface indicating stability for some while continuous growth for 
others. The time rate of change of pressure is responsible for acoustic 
radiations. A larger height of the surface may further reduce the growth of the 
disturbances, but a resulting sound field could be larger. In this context, the 
resulting sound can be viewed as a source of instability. 


AT y/bj-10 



y. ft 
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ACTIVE CONTROL - EFFECT ON DISTURBANCE GROWTH 


Active control by surface heating is modeled by modifying the temperature 
boundary condition over a small portion of the surface. The control temperature is 
applied for a finite time interval coinciding with the passage of the wave packet 
over the control strip location. The control strip is located in the region of 
unfavorable pressure gradient at x = 0.9 feet. This location was selected because it 
is a region of large growth and also because the coupling of the control with the 
flow is favorable. 

The control is not optimized. Larger reductions in growth can be achieved by 
optimizing amplitude and phase as well as using multiple control strips. 


(]l^f 



.0 2 p ^ J >)>}>>} > I I' } 

o V > ^ ^7) ^^ Control strip 

I ^ 1 . 1 . I I I i 

0 .8 1.6 2.4 3.2 

X, ft 



i 
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DETERMINATION OF ACOUSTIC FIELD 

Due to the vastly different length scales between the acoustic far field and the 
boundary layer, it is not feasible to compute the acoustic far field concurrently 
with the Navier-Stokes simulation. The acoustic analogy introduced by Lighthill 
(ref. 1) provides a formalism to couple wave equations in order to determine the 
acoustic far field with the near field behavior of the viscous flow. The unsteady 
near field (generally including turbulent fluctuations) is used as a source term for 
a convective wave equation describing the far-field acoustic behavior. Our approach 
is to solve the linearized Euler equations using the pressure computed from the 
Navier-Stokes solution as a time-dependent boundary condition. We ignore the 
curvature of the surface because it is considered negligible On the acoustic length 
scale. 


• Due to vastly different length scales between the 
acoustic far-field and the boundary layer, it is not feasible to 
compute acoustics from the Navier-Stokes solution 

• The acoustic field is evaluated by solving linearized 
Euler equations using pressure from the N-S as time 
dependent boundary condition 

• This approach is an exact instance of the Lighthill's 
acoustic analogy 

(Ux + '^y) 

Pt + Px + 

Uj -h Ux + Px - 0 

Vt -I- Vx -h Py - 0 

B.C. pCx.o.t)-^ known from the N-S solution 
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WALL PRESSURE PERTURBATION ALONG THE SURFACE 


The acoustic radiations, which depend on the time rate of changes within the 
flow, are present even when the flow is being stabilized. In addition, acoustic 
radiations are known to produce instability leading to transition. Therefore, it is 
important to study not only the flow instability but the acoustic radiations as well. 

The pressure perturbation versus time is shown at different stations along the 
wall. This is used as an input to the linearized Euler equations to calculate the 
acoustic field. It may be noticed that the maximum amplitude and bandwidth are 
around x = 1.0 ft., which is on a location corresponding to maximum growth. 
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POWER SPECTRA OF WALL PRESSURE PERTURBATIONS 

The power spectral density of the wall pressure shows that in addition to the 
two input frequencies comprising the wave packet, two other frequencies appear 
between x = 0.8 ft. and 1 .2 ft. carrying most of the energy. This is a result of the 
instability produced by the concave surface. In the convex region and the flat 
region that follows, the energy content of these two additional modes is reduced 
drastically and the two original frequencies reappear at x = 1.6 ft. 


p(f)* 



P(f)* p(f)* 




p(f)* 
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DECAY OF TOTAL POWER WITH DISTANCE FROM THE SURFACE 


In order to establish the near- and the far-field behavior, total acoustic power 
is plotted against the distance from the wall. Two distinct regions can be 
identified, one close to the boundary (y ~ 5 ft.) where the decay is exponential and 
the other farther away with an algebraic decay, the acoustic far field. 
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ACOUSTIC PRESSURE AT 1 1 FT. RADIUS 


The acoustic pressure at a radius of 11.0 ft. from the surface is shown at 
various angular positions between 30° and 150°. Two distinct features are 
noticeable. One is the sound resulting from the input instability. This reaches the 
far field at different times due to convection. The other is the sound radiated due 
to the growth and decay of the wave packet. This component has strong beaming 
characteristics in the downstream direction. In addition, the modulation resulting 
from the passage of the wave packet on the curved surface can be clearly seen. Such 
strong amplitude and directivity have special significance when originating from a 
wind tunnel contraction since they may trip instability at the model's leading edge. 
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PSD OF ACOUSTIC PRESSURE AT 11 FT. RADIUS 


The power spectral density (PSD) of the radiated acoustic field is shown at 
11.0 ft. from the surface at various angular positions. Most of the sound is 
directed downstream with a widening in bandwidth near 0 = 150®. The predominant 
radiation frequencies are those resulting from the growth and decay of the wave 
packet. 





SEPARATION CONTROL BY SOUND OVER AN AIRFOIL 


This section concerns an experiment designed to study the effect of sound 
emanating from the surface of an airfoil on the separated flow at higher angles of 
attack. At high angles of attack the free shear layer, which sheds from the point of 
separation, displays the characteristic instability, whereby a sequence of discrete 
vortex forms. The aim is to study the coupling between the sound and the flow and 
its effect on separation. This coupling is utilized to Increase circulation, stall 
angle, as well as to reduce the drag of an airfoil. 

The experiment was conducted in the 2- by 6-inch wind tunnel at low Reynolds 
number. Smoke visualization and a hot wire were used to study the flow. 


• Experimental study of the effect of acoustic perturbation 
on pressure distribution, iift and drag 

• Shear iayer is very sensitive to the excitation field 
emanating from the surface in the vicinity of the separation 

• This sensitivity between sound and the locai field is used 
to impart dramatic favorabie changes to the flow 



0 .2 .4 .6 .8 1.0 

x/c 
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RESULTS - BEHAVIOR OF THE SHEAR LAYER 


At hiRh angles, a = 15 and 20 degrees, the flow separates over most of the upper 
surface Typical velocity power spectra at different locations over the upper 
^raoe or tS^alrfoil are sho»n (left). At x/o - 0.02. only a low-frequency peak la 
Observed. This peak frequency is approximately 50 Hz and corresponds to the 
"shedding” frequency of the wake region. Further downstream ® 

x/c = 0.22, other peaks appear with frequencies of 125, 250, and 500 Hz. This 
observation suggests that the separated shear layer forms and develops with a 
fundamental frequency of about 500 Hz at a distance away from the leading edge. We 
ScraSi t^a? re^redcmlnant frequcnolea In the shear layer vary with the second 
power of freestream velocity (right figure). 



612 


3U.CK OOBAPH 


RESULTS - SMOKE 


VISUALIZATION 


A smoke visualization technique is used to see the effectiveness of control by 
sound. For the uncontrolled case, a large region of separation can be seen on the 
upper surface. Separation control is achieved by sound injected from a narrow slot 
located at x/c = 0.15 on the upper surface. When forced exitation at one of the 
subharmonic frequencies is applied, the separated region becomes drastically 
reduced. The figure shows the controlled and the uncontrolled pictures. 

It is remarkable that a single injection gap can induce such a dramatic change 
in reducing the area of separation in contrast to the relatively modest control that 
is achieved at small angles of attack. 

\ 


Rec == 3 5 X 1 f 100 Hz 

UNCONTROLLED CONTROLLED 




or tWRliUALrnf^— 
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CONTROL OF LIFT AND WAKE PROFILES 


The effectiveness of control is plotted in terms of the lift coefficient versus 
the angle of attack and a typical wake profile at one angle of attack. The gain in 
lift over the whole range of angle of attack is apparent. Similarly, the defect in 
the wake profile becomes much smaller, indicating a reduction in drag. 


Reg=3.5 X 10^ 

A Controlled 
o Uncontrolled 



I 
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CONCLUDING REMARKS 


• Flow control using active surface heating and sound is a powerful 
and effective technique 


BENEFITS 

•Development of low drag configurations by control for laminar and 
turbulent flows 

•Optimization of control technique for flows with periodic and 
random behavior 

FUTURE DIRECTIONS 

•Extension to turbulent and high speed flows 

•Application to realistic aircraft configurations 

•Design of wind tunnel contractions with minimum flow and 
acoustic interference 
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NEAR-FIELD NOISE PREDICTIONS 
OF AN AIRCRAFT IN CRUISE 


John W. Rawls, Jr. 
PRC Kentron, Inc. 
Hampton, Virginia 



INTRODUCTION 


The application of laminar-flow wing designs will offer significant drag 
reduction and fuel savings on large commercial transports. Traditionally 
surface contamination, surface defects, vibration and turbulence have 
prevented practical application on commercial aircraft. These concerns are 
being overcome through innovative research. 

Noise generated by the aircraft is now being investigated as a cause of 
transition of a laminar boundary layer. An acoustic disturbance 
with the correct frequency and intensity is known to cause transition of a 
laminar boundary layer under laboratory conditions. It is speculated that 
aerodynamic noise sources may have the correct characteristics to cause 
boundary-layer transition under flight conditions. 

One important element in determining whether a laminar boundary layer 

presence of an acoustic source) is an accurate 
prediction of the noise source and its characteristics. Predictions yield not 
only the total noise level, which can be verified by measurements, but also 
informatpn concerning the characteristics of the individual noise sources. 

The prediction of the near-field noise environment of an aircraft in cruise 
has become an important element in the research and development of wing 
designs that can maintain laminar flow at a cruise speed of Mach 0.8. 


CAUSES OF TRANSITION 

• Surface Contamination 

Insect Residue 
Ice Crystals 

• Surface Defects 

• Vibration 

• Turbulence 

• Noise 


LOCKHEED NEAR-FIELD NOISE PREDICTION PROGRAM 


A computer program was developed by Lockheed-Georgia Company for the 
purpose of predictinq the near-field acoustic environment of an aircraft in 
cruise (Refs. 1, 2, 3). This program is designed to compute the 1/3 octave 
band spectra and Overall Sound Pressure Level (OASPL) at a point on an 
aircraft operating at high altitude and high-speed conditions. 

Noise sources predicted by the code can be divided into two categories: 
propulsion sources and airframe sources. The propulsion system of the 
aircraft is assumed to be a turbojet or turbofan engine. The components of 
the engine which produce noise are fan/compressor, combustor, turbine and 
jet. Acoustic suppression material may be used in the engine and is accounted 
for in the prediction procedure. 

The airframe itself is also a producer of noise. The two airframe 
sources predicted by the code are trail ing-edge noise and turbulent boundary- 
layer noise. 

The noise source algorithms are based on empirical formulations derived 
from data at near sea level conditions with relatively low or no forward speed 
effects. Since the noise predictions are for high speed and high altitude, 
corrections for forward speed and altitude are required. 

A forward speed correction is required to determine the correct emission 
angle and propagation distance at the time the sound is emitted. The 
correction is a function of the observed angle and distance at the time the 
sound is received. Forward speed corrections are also made to the OASPL to 
account for convective and dynamic amplifications. (A detailed discussion of 
these forward speed corrections can be found in Ref. 1.) 

An altitude correction is made by adjusting the acoustic impedance to 
correspond to the ambient conditions at altitude. Other altitude corrections 
such as atmospheric attenuation and atmospheric absorption are assumed to be 
small in the near field. 
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NEAR-FIELD NOISE SOURCES 


For the purpose of this report, the propulsion system will be a 
high-bypass turbofan engine. This propulsion system is widely used in large 
commercial transports where laminar-flow wing designs would be most 
beneficial. Noise prediction procedures are provided for the four basic jet 
engine components i.e., the fan/compressor, combustor, turbine and jet. With 
these basic components the noise generated by a high bypass turbofan can be 
modeled. 

Fan noise is generated by the fan blade interacting with the incoming 
flow. Noise from the fan can radiate forward through the inlet and aft 
through the secondary jet nozzle. In addition to the broadband component of 
the noise spectra, discrete tones occur at the blade passing frequency and its 
harmonics. If the tip speed of the fan blade is supersonic, subharmonics of 
the blade passing frequency are also generated which only propagate forward. 
Compressor noise is predicted by the same method used to predict fan noise. 
However fan noise usually dominates compressor noise and is not considered in 

this report. 

Combustor noise and turbine noise both propagate through the primary 
nozzle. Combustor noise is a low-frequency source while turbine noise is a 
high-frequency source. Like fan noise, turbine noise is composed of a 
broadband spectra and discrete tones occurring at the blade passing frequency 
and higher harmonics. Typically, the blade passing frequency of the turbine 
is so high that only the fundamental tone is included in the spectra. 

Jet noise, which is made up of three components, propagates from both the 
primary nozzle and the secondary nozzle. The first component of jet noise is 
mixing noise. Mixing noise is always present in a jet and is the result of 
jet interacting with the ambient air. The second and third components of 
jet noise are broadband shock and screech shock noise. Both broadband shock 
and screech shock occur only when the jet velocity exceeds the local speed of 
sound in the jet. The onset of broadband shock noise is predicted as soon as 
the jet flow conditions are supersonic. The onset of screech shock noise is 
not well defined and can usually be controlled with an appropriate nozzle 
design. Consequently screech shock noise is not considered in this report. 

The two airframe sources, turbulent boundary-layer noise and trailing- 
edge noise are also included for study in this report. 


620 


I I II \W¥ I . I lllllli III! I Mill mil hill ,11 II. I 111 ; II 11111111.11 illJ I ill 1.1 ..I lillJIIJII IlIUIUililiillUlllilltiiB 



NEAR-FIELD NOISE SOURCES 


TURBULENT BOUNDARY LAYER NOISE 



AFT-RADIATED FAN NOISE 
JET MIXING NOISE 
BROADBAND SHOCK NOISE 
SCREECH SHOCK NOISE 


SECONDARY JET 
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HIGH BYPASS TURBOFAN MODEL DESCRIPTION 


Input parameters required for the noise prediction of a high bypass 
turbofan engine are obtained using estimated performance parameters from a 
full-scale high bypass turbofan with a maximum continuous thrust rating of 
38,000 lb. Temperature and pressure ratios are known for the fan discharge, 
the combustor inlet and the low pressure turbine exhaust as a function of 
flight Mach number, altitude and thrust setting. The corrected rotational 
speeds of the fan and high pressure compressor, the corrected fan mass flow 
rate, the bypass ratio, and the fuel mass flow rate are also known. With this 
information, the necessary input data for a complete noise prediction of a 
high bypass turbofan are computed. 

Predictions are made for two observer locations on a line parallel to the 
engine axis as shown. 


NOISE STUDY OF A FULL SCALE 
HIGH BYPASS TURBOFAN ENGINE IN CRUISE 

Engine configured without acoustic suppression | 

I I 

Engine performance data given as a function of | 

I I 

* Aititude 

• % Thrust • 

•Fiight Mach No. 
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ALTITUDE VARIATION 


Noise levels are investigated by varying the altitude from 20,000 feet to 
40,000 feet while operating at maximum continuous thrust and a flight Mach 
number of 0.8. Parameters relevant to noise predictions and engine 
performance vary significantly with altitude. For example the mass flow rate 
required by the fan at 40,000 feet is 46% less than that required at 20,000 
feet. Since fan noise correlates directly with the mass rate a reduction in 
the OASPL as a function of altitude is realized. 

It would seem reasonable to expect total OASPL to also decrease with 
increasing altitude. However, this was not the case. Two offsetting factors 
contribute to this result. First, the acoustic impedance decreased with 
increasing altitude. Second, the intensity of the broadband shock noise, 
which is determined by the jet Mach number, increased as a function of 
altitude by about the same amount as the acoustic impedance decreased. This 
resulted in a constant OASPL as a function of altitude. 
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ALTITUDE VARIATION 

Mach 0.8 

100% Maximum Continuous Thrust 
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BOEING 757 GLOVE TEST DATA 


During the Boeing 757 glove test, data were obtained over an altitude 
range of 30,000 feet to 41,300 feet. The OASPL is plotted as a function of 
the secondary jet Mach number as shown. In this figure, the OASPL maintains a 
level of 126 dB for secondary jet Mach numbers between 1.0 and 1.28. When the 
jet Mach number exceeds 1.0, broadband shock noise is present and follows the 
trend of constant noise level as a function of altitude. When the jet Mach 
number is less than 1.0, the noise level is reduced as would be expected since 
broadband shock noise is not present in the secondary jet. 


Mach 0.8 

Altitude 30,000 / 41,300 Feet 
Microphone 10 



SECONDARY JET MACH NUMBER 


MICROPHONE LOCATION 
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THRUST VARIATION 


Noise levels are also investigated by varying the thrust of the jet while 
operating at a flight Mach number of 0.8 and an altitude of 30,000 feet. As 
expected the total OASPL increases with increasing thrust. However the 
increase is not as great as might be expected. A reduction of 64% in thrust 
results in only a 5 dB decrease at 30“ and 6 dB at 120°. The dominate noise 
source is again broadband shock noise from the secondary jet. In fact the 
remaining propulsion noise sources contribute very little to the total OASPL. 
Examination of the jet conditions indicates that additional thrust is 
generated by the primary jet with the jet Mach number increasing from .96 at 
64% thrust to 1.17 at 100% thrust. The secondary jet conditions on the other 
hand increase from 1.11 at 64% thrust to 1.18 at 100% thrust. It is worth 
noting that the secondary jet conditions are supersonic even at low thrust 
accounting for the presence of broadband shock noise. 
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TYPICAL JET CONDITIONS DURING CRUISE 


Broadband shock noise from the secondary jet has emerged as a very 
important noise source during cruise. Even though the secondary jet velocity 
at total temperatures is much lower than the primary jet, the Mach number in 
the secondary jet is typically supersonic when varying altitude and thrust. 

The supersonic conditions in the secondary jet are primarily the result of the 
low total temperature in the jet. Not only is broadband shock noise typically 
present in the secondary jet but the secondary flow acts as a shield to reduce 
the effect of the noise levels produced by the primary jet. 


1000 ft/sec < Vg < 1200 ft/sec 
500 °R <Tg < 575 °R 


1.00 < Mjg < 1.15 
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FLIGHT MACH NUMBER VARIATION 


Finally noise levels are examined by varying the flight Mach number from 
0.5 to 0.9 in 0.1 increments. Thrust and altitude are maintained at 100% and 
30,000 feet, respectively. At a flight Mach number of 0.5, broadband shock 
noise is absent from both the primary and secondary jets. An increase in the 
flight Mach number to 0.9 dramatically increases the total OASPL at both 
observer locations. An increase of 26 dB is observed at 30“ and an increase of 
12 dB is observed at 120°. This increase in noise is the result of several 
factors. One factor is that the engine performance parameters are at higher 
levels as the flight Mach number increases; another is that the forward speed 
correction for the noise sources is larger as the flight Mach number 
increases. 
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Flight Mach Number 
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AIRFRAME NOISE SOURCE SPECTRA 


Airframe sources are important in determining the total OASPL especially 
in cases where a surface is shielded from the propulsion system. Recent 
measurements on the wing of a Boeing 757 showed that noise levels on the upper 
surface were unaffected by increases in the engine RPM which is mounted by a 
pylon from the wing lower surface. This spectrum shows the contributions of 
turbulent boundary-layer noise and trailing-edge noise at an upper surface 
point near the leading edge of a Boeing 757 wing. The turbulent boundary- 
layer noise levels predicted by this code are typically 10 to 15 dB lower than 
the trailing-edge noise levels as shown. Therefore turbulent boundary-! ayer 
noise does not contribute to the total OASPL. 

Altitude 40,000 Feet; Mach 0.8; MIc 11 



.05 ,10 .20 .40 .80 1.60 3.15 6.30 12.50 


Frequency, 1/3 Octave KHz 
MICROPHONE LOCATION 
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TRAILING-EDGE NOISE 


Trailing-edge noise is the result of a sudden change in the pressure as a 
boundary layer is convected past a surface trailing edge. As the boundary 
layer approaches the trailing edge, the pressure differential between the 
upper and lower surface is forced to zero to satisfy the Kutta condition. 

This, in turn, results in an induced pressure field which propagates away from 
the trailing-edge region (Ref. 4). 

The algorithm for the prediction of trailing-edge noise is based on a 
flat plate analysis. The observer is restricted to be on the centerline of 
the plate and the flow must leave perpendicular to the trailing edge. 


Flow Observer 



Source Location 



Source/Observer geometry for trailing edge noise 
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TRAI LING-EDGE NOISE 


Trai ling-edge noise is a distributed source and not conducive to point 
source analysis. The wing of a Boeing 757 therefore Is approximated as a 
series of adjacent flat plates. Contour levels are generated on a grid of 121 
chordwise and 201 spanwise stations. This corresponds to a point every 
3 inches in both the spanwise and chordwise directions and gives sufficient 
resolution for smooth contour levels. 



Flat plate approximation of a Boeing 757 wing 



TRAILING-EDGE NOISE 


Trailing-edge noise is a function of primarily two factors, the flight 
Mach number and the thickness of the turbulent boundary-layer at the trailing- 
edge. An increase in the flight Mach number significantly increases the noise 
levels. An increase in the flight Mach number from 0.5 to 0.8 for example 
will increase the peak Overall Sound Presure Level by 21 dB. 

The effect of varying the turbulent boundary- layer thickness is perhaps 
of greater interest to laminar-flow studies, since a laminar boundary layer on 
the leading-edge portion of the wing will reduce the thickness of the 
turbulent boundary-layer thickness at the trailing edge. A laminar boundary 
layer is simulated over the first 40^‘ of the wing by using a length scale of 
60% of the chord to calculate the turbulent boundary-1 ayer thickness. The 
resulting reduction in noise level is unfortunately not very significant. 

Peak noise levels are reduced by only 3 dB. 



100% Turbulent Flow 
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MODELING THE ACOUSTIC FIELD 


The total noise at a point on the surface of the aircraft is computed by 
summing the contribution of the propulsion and airframe sources. In the case 
of the propulsion sources, the predictions are free-field so that the presence 
of a wing surface is not accounted for. The effect of shielding by a surface 
is also not predicted but can be simulated by judiciously choosing the noise 
sources to be predicted. 
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COMPARISON OF MEASURED VERSUS PREDICTED NOISE 


A noise prediction is made and compared with test data obtained on the 
wing of a Boeing 757 aircraft. The microphone is located on the lower surface 
of the wing near the leading edge of a test section designed for laminar 
flow. Three sources contribute to the total OASPL. They are jet mixing noise 
and broadband shock noise from the primary and secondary jets. In the low 
frequency range, .05 to .4 KHz, jet mixing noise is the dominate source. At 
higher frequencies, the spectra are dominated by broadband shock noise from the 
primary and secondary jets. Fan, combustor, turbine, and trail ing-edge noise 
were included in the noise prediction but did not contribute to the total 
OASPL in this case. 

Comparison of the total predicted noise and the measured noise shows 
levels are within 10 dB except for very high frequencies. It is believed that 
with sufficient test data, improvements to the code can be made especially in 
the area of forward flight corrections. 

Altitude 40,000 Feet; Mach 0.8; Fan Speed 4100 RPM; MIc 13 
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CONCLUSIONS 


The physics of the coupling of sound waves with the boundary layer is not 
yet well understood. It is believed; however, that for effective coupling of 
the sound waves and instability waves in the boundary layer, a matching of 
both frequency and wave number must occur {Ref. 5). This requires that the 
sound field to be accurately defined in both space and time. Currently 
analytical prediction methods lack sufficient accuracy to predict the noise 
levels from components of a turbofan engine. Although empirical methods do 
not yield the detail required for an analysis of the receptivity of sound by a 
boundary layer, valuable insight can be gained as to the changes in noise 
levels that might be expected under various operating conditions and aircraft 
configurations. 


Predictions 

Noise levels remained unchanged with Increasing altitude when 
flight Mach number and thrust are held constant 
at O.B and 100% respectively 


Noise levels Increase significantly with Increasing 
flight Mach number 


Noise levels Increase moderately with Increasing thrust 


Flight Data 

Limited flight data tend to support these conclusions 

although noise prediction program tends to over predict 
the total noise level 


Significant noise sources 


• Jet Noise 


Primary Jet 

Broadband Shock 

Secondary Jet 
Primary Jet 

Jet Mixing 

Secondary Jet 


• Forward Radiated Fan Noise 


> 'R’alllng Edge Noise 

Dominant noise sources dependent on flight Mach number 

Under normal cruise conditions broadband shock noise 
can be expected to be present In the secondary |et 
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